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Instrumentation; Miniaturization; Moisture content; Porosity = | = | 4 


ABSTRACT: Miniaturization of instrumentation has opened up the use of borehole 
geophysical techniques to geotechnical consultants, but most are not yet acquainted 


with the capabilities of the procedures. Borehole geophysics has been used widely in 


studies of nuclear power plant sites, waste disposal problems, and local and _— by 7 
-geohydrological relationships. Most investigations have been oriented toward 


| Stratigraphic correlation and assessment of lithologies. However, there is a wealth of } 


semi-quantitative or quantitative data available which is often overlooked. Because of 
the environments in which the geotechnical specialist works, exploratory borings are _ 


most frequently cased, and the radiation group of geophysical devices must be used. In 
properly engineering borings and with appropriate calibration, these devices can — 


REFERENCE: Crosby, James w., Ill (Prof. and Head of Geological Engrg. Section, _ 


Dept. of Civ. & Environmental "Engrg. , Washington State Univ., Pullman, Wash. S 


99164), Konstantinidis, Byron, and Davis, Paul, “Geotechnical Applications of 


Borehole Geophysics,” Journal of the Geotechnical Engineering Division, re. Vol. 
107, No. GT10, Proc. Paper 16585, October, 1981, pp. 1255- ee 


16591 DETERMINATION OF OF SOIL MODULI Trl THAN? 


' KEY WORDS: Damping; Dynamic properties; Field tests; Foundation 
Stability analysis; Laboratory tests; Shear modulus; Soil stresses; Strain; meee i 4 


' techniques. A supplementary set of laboratory resonant column tests was performed 


ABSTRACT: A field geophysical investigation was s performed to determine the shear — x 
and Young’s moduli as a function of depth for a site in southern Ohio. = 
q 


with the Drnevich resonant column device. Laboratory undisturbed specimens were ‘ 
excited in both the longitudinal and torsional modes to obtain both shear and Young’s ¥ _” 
moduli as a function of strain. Both laboratory and field data are presented and 
REFERENCE: Marcuson, William F., II] (Research Civ. Engr., "Geotechnical 
Laboratory, U.S. Amy Engr. Waterways Experiment Station, Vicksburg, Miss), and 
Curro, Joseph R., Jr., “Field and Laboratory Determination of Soil Moduli,” Journal -4 
the Geotechnical Engineering Division, ASCE, Vol. 107, No. GT10, Proc. "ior. 


investigation included crosshole, downhole, surface refraction, and surface vibratory 


me WORDS: measurement; Alaska; Flow 
Geological surveys; Geophysical exploration; Geophysical 
Marine engineering; Ocean bottom; Offshore structures; Sand waves Ss 7 


- ABSTRACT: Acoustic (geophysical) methods were applied to investigate ‘the sea-floor oO 
geological conditions, particularly sand waves, which might place constraints on 
construction in an area of potential offshore development in lower Cook Inlet, Alaska. 

_ Multisystem, high resolution, marine acoustic survey data were collected using 9 a 
sparker, tuned transducer, acoustipulse, fathometer, and side-scan sonar. Based on 65 

- sea- floor samples obtained in the northern portion of the study area, the sea-floor soils 

- consisted of over 90 percent sand by volume. This portion of lower Cook Inlet 

- contained the most numerous sand waves and contained sand waves of 5 ft (1.5 m) or 7 
_ greater height over an 80 percent length of the survey lines. The possibility of sand * 
wave movement is examined by reviewing the two phase (water-sand) flow regime and 

by repeat geophysical profiling. The present flow regime appears to be conducive to io P| 
_ the formation of ripples and dunes, but tidal flow durations may be insufficient for the 

formation of large bedforms. £3 


"REFERENCE: Mahmood, Arshud A. (Engr. McClelland Engineers, 
Houston, Texas), Ehlers, Clarence J., and Cilweck, Blase A., “Sand Waves in Lower | “ft 
Cook Inlet, Alaska,” Journal of the Geotechnical Engineering Division, ASCE, Vol. 
107, No. GT10, Proc. Paper 16583, October, 1981, pp. 1293-1307 Ok 
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KEY WORDS: Bulk Compression waves; Geophysics: ‘Laboratory w 


tests; Rayleigh waves; Seismic studies; Seismic tests; Shear modulus; Shear seth 
waves; Soil dynamics; Velocity 

_ ABSTRACT: Applications of scaled-down seismic cross-hole and steady-state Rayleigh 
wave tests to characterize dynamic soil properties in laboratory facilities are described. _ 
The nondestructive nature of seismic tests make them highly advantageous in 7 
laboratory settings. The major changes in test procedures and equipment to scale-down é. 
from field seismic tests are associated with: (1) Size and energy content of seismic | 
source; (2) receiver size and sensitivity; (3) wave speed timing; (4) electronic triggering 7 L 
requirements; and (5) precision of distance or spacing measurements, or both. The :Y 
seismic cross-hole test proved to be most effective in determining the spacial and — 
temporal variation of elastic soil properties and proved to be a very sensitive indicator ( 
of small changes in confining pressure. In addition, a sand bin testing facility which is ot 


fitted with air bags to soines surface confining pressure and new torsional cross-hole 


16599 SEISMIC TECHNIQUES IN THE LABOR: 


REFERENCE: Woods, Richard D. (Prof. ‘of Civ. Engrg., Univ. of Michigan, Asa 
Arbor, Mich.), and Henke, Robert, “Seismic Techniques in the Laboratory,” Journal of © 
the Geotechnical Engineering Division, ASCE, Vol. 167, No. GT10, eral Paper 7 


tests; Piers; Shafts; Soil Structural analysis; 


ABSTRACT: Compression wave propagation methods. provide an economical and 


rapid nondestructive means of assessing the structural integrity of drilled shafts shortly 
after construction. Several source-receiver configurations can be employed. The most 
successful configuration used involved a mechanical hammer with which a vertical | 
impulse was applied at the top of the shaft and vertical velocity transducers which | 
were cast in place in the shaft at preselected depths. Successful application of this 
configuration required identification of reflected compression wave arrivals with — 
transducers allowing travel times of reflected waves to be used with tg so gvecedl 
velocity to determine elevations of reflection points. Examples of the results for tests 

with sound shafts and for shafts with known defects of different sizes are presented. _ 

REFERENCE: Hearne, Thomas M., Jr. (Research Asst., Dept. of Civ. Engrg., Univ. of 
Texas, Austin, Tex. 78712), Stokoe, Kenneth H., Il, and Reese, Lyndon C., “Drilled- 

Shaft Integrity by Wave Propagation Method,” Journal of the ~ Geotechnical 
_ Engineering Division, ASCE, Vol. 107, No. GT10, Proc. Paper 16582, October, 1981, 

pp. 1327-1344 


16590 ESTIMA ATION OF SF SPT-N AND DENSITY 

KEY WORDS: testa; Probability theory; Relative 


density; Relative density determination; Sands; Soil density; Soil 7 
_ engineering; Standard penetration tests; Statistical analysis aye 
ABSTRACT: A consistent set of models is developed for inference of the relative © 
density (D;) from the blowcount number of the SPT (N) under both laboratory and 


ss conditions. The model in the laboratory (reference model) is based on a large data 


_ set, which consists of essentially all available laboratory information on joint 
measurements of N, D,; and vertical effective stress. It takes explicitly into account 
errors in the laboratory determination of D; and is of a form which makes its use 
particularly simple. The reference model is not directly applicable in situ. However, the 

_ model can be applied in the field after a simple, but sometimes depth dependent, — 
modification which requires a few measurements of Nand atthesite. 

_ REFERENCE: Fardis, Michael A. (Asst. Prof. of Civ. -Engrg., Massachusetts Inst. o 
Technology, Cambridge, Mass.), and Veneziano, Daniele, “Estimation of SPT-N and 
_ Relative Density,” Journal of the Geotechnical Engineering Division, ASCE, Vol. 107, 

| No. GT10, Proc. Paper 16590, October, 1981, pp. 1345-1359 as ie meer) | 
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16604 STATISTICAL ANALYSIS OF SAND LIQUEFACTION 

= WORDS: Laboratory tests; Liquefaction; Relative ve density; Sands; 
Sandy soils; Shear stress; Shear tests; Statistical analysis > 
ABSTRACT: A consistent set of stochastic models is developed for the liquefaction me 
resistance of a homogeneous mass of sand, such as a laboratory sample or small 
element of in situ material. A probabilistic model for liquefaction in the laboratory for i 
undrained simple shear tests is developed first. This model gives the probability 
distribution of the number of cycles to liquefaction as a function of uniform shear 
stress amplitude, mean initial effective stress, and relative density. Corrections are 
quantified and incorporated into the model. Parameters are estimated from a large 
number of simple shear test results obtained by eight different groups of researchers 
using several sands and testing techniques. The laboratory liquefaction model is Ep 
‘converted into one applicable in the field. Finally, the model is to load 

cycles with nonuniform shear stress amplitude. 

REFERENCE: Fardis, Michael N. (Asst. Prof., Dept. of Civ. Engrg., Mass. Inst. of 
Technology, Cambridge, Mass.), and Veneziano, Daniele, “Statistical Analysis of Sand 

' Liquefaction,” Journal of the Geotechnical Engineering Division, ASCE, Vol. 107, No. 

| GT10, Proc. Paper 16604, October, 1981, pp. 1361-1377 


KEY WORDS: Deformation; Finite elements; 

i Ye 


"ABSTRACT: may element of solution caverns | in rock salt 


_ were performed to study the long-term behavior of such caverns. The effect of cavern 
shape (height-to-diameter ratio) and the internal pressure on the time- -dependent 
behavior of caverns were the main objectives. The pattern of creep-induced stress — 
changes around the caverns were studied, and the main aspects were demonstrated 
through the results of the analyses. The cavern volume closures can be expressed by a 
_ simple equation as a function of time, and the accuracy of this expression was verified 
_ by the results of the analyses. The other aspects of the time-dependent behavior of the 
| caverns studied were: (1) the influence of a drop in internal pressure; (2) the influence 
of creep deformations on the results and interpretation of pressure tests for coy 
: — and (3) the influence of gradual cavity growth due to the leaching process as ~ 
- compared to the analysis of cavity assumed at full size initially Piney 
REFERENCE: Ghaboussi, Jameshid (Assoc. Prof., Dept. of Civ. . Engrg., Univ 
Illinois at Urbana-Champaign, Urbana, Ill.), Hendron, Afred J., Jr., and rr 


Randall E., “Time-Dependant Behavior of Solution Caverns in Salt,” Journal of the 
- Geotechnical Engineering Division, ASCE, Vol. 107, No. GT10, Proc. Paper 16597, 


16588 ROAD ON PEAT: OBSERVATIONS AND DESIGN hoa 
KEY WORDS: Consolidation; Deformation; Models; Peat; — Peat bog; red 
ABSTRACT: In 1953, a — was reconstructed on a uniform peat bog. Certain. lengt s 
were overstressed by an excessive thickness of road material, and in places it was — 
necessary to replace the gravel filling with a light-weight substitute. Road levels and 
settlements have since been recorded over a period of 26 years. Following an 
| investigation in 1963, a simple method of design was formulated which requires 
5 settlement to be separately estimated under three headings: elasiic (immediate), 
BF consolidation, and creep. These are then added together to give total settlements. 
| Satisfactory agreement was found between observed settlements and those values 
‘| estimated by the simple method. This agreement is again confirmed after 16 years. 


_ REFERENCE: Hanrahan, Edward ‘A (Assoc. Prof. of Civ. Engrg., University College, 
Dublin, Ireland), and Rogers, Martin G., “Road on Peat: Observations and Design,” 
Journal of the Geotechnical Engineering ‘Division, ASCE, ie 107, No. GT10, Free. 


; 
H 
= 


hed 


be prepaid. 


U. S. Customary-SI FacTors 


bo compute the SI unit values of measurements. A complete guide to the 


4879 


= (International System) units, the following list contains conversion factors ¢ to enable readers 
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James W. Crosby III," M. ASCE, Byron Konstantinidis,’ A.M. ASCE, 


. = development for nearly half a century. In these applications it has proved to 
be a tool of great value, but one that is excessively expensive by geotechnical — 
engineering standards. During the last two decades, however, the general — . 
_ availability of inexpensive, miniaturized electronic components has placed these 
_ same capabilities at the disposal of the geotechnical c consultant. However, the 
potential for application of borehole geophysics i in . geotechnical engineering is 
_ far from being fully realized. This is primarily because geotechnical engineers, . 
_ in general, are not familiar with the capabilities of geophysical borehole logging — 
techniques. The purpose of this paper is to demonstrate the capabilities and 7 
: anaes’ of geophysical logging as a tool for obtaining both a and 
Geophysical borehole logging has been used extensively by the writers 
‘stratigraphic characterization and correlation. However, quantitative and semi- 
_ quantitative information about physical properties such as in-situ moisture — 
content, density, porosity, and degree of saturation can also be derived from 
geophysical logs. With proper site specific calibrations, data derived from 


geophysical logs can be used with a . limited amount of laboratory test data 
to provide a a detailed picture of subsurface conditions. _ 
_ The methodology described in this paper draws upon a a well-established body 
of theoretical knowledge and a variety of studies conducted by the writers 
over an eight-year period. Most of the field applications have been related 
to siting studies for nuclear power plants, but others include missile site studies, 


a Presented at the October 22-26, 1979, ASCE Annual Convention and Exposition, held 
 ' Prof. and Head of Geological Engrg. Section, Dept. . of Civ. and Environmental Engrg., 
Washington ‘State Univ., Pullman, Washo ond? 
wa Sr. Engr., Fugro, Inc., Long Beach, Calif. faa 
*Chf. Geologist, Medall, Aragon, Worswick & Associates, Inc., Santa Ana, Calif. 
ap _ Note .—Discussion open until March 1, 1982. To extend the closing date one month, z 
_awritten request must be filed with the Manager of Technical and Professional Publications, _ 
_ ASCE. Manuscript was submitted for review for possible publication on February 17, 
r F 1981. This paper is part of the Journal of the Geotechnical Engineering Division, Proceedings — 
of the American Society of Civil Engineers © ASCE, 107, No. GT10, October, 
1981. ISSN 0093-6405 /81 /0010-1255 / $01.00. 
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accomplish various types of downhole measurements. Many are e limited fede 
tion devices, and only a few are routinely used in geotechnical engineering 
studies. Among the most frequently used tools are those which comptins the a , 

a radiation and electric groups (Fig. 1). The radiation sondes can be used in 

: ~ cased or uncased holes which are air or fluid- filled. Electric logging tools are = 
generally run in fluid-filled, uncased holes. 
Radiation Logs. —The radiation logs are of great interest in 1 geotechnical practice, 

both from the of desired information and required flexibility. 


1.—Full Suite Geophysical Logs (Sundesert Boring B-DH- 35) 
photons or neutrons which have the ability to penetrate matter. . The gamma — 
can be generated naturally by -Fadioactive elements in the tested 
formations or by fabricated gamma- emitting sources. The neutrons are generated | 
exclusively by introduced sources. 
Ss In counting radiation pulses, some form | of ratemeter is utilized. Each ratemeter 
i! a characteristic dead time or a time increment between near-coincident _ = 
pulses less than which the meter cannot discriminate between them. The _ 
relationship canbe expressed 


on 


my Dead times for each wn are readily determined in the laboratory as part 
of the instrumental calibeation Procedures (Gray and Preiss, rabies Ratemeters 


= tine time of tl ‘the ‘ratemeter; R count rate; = 


resented 


7 
t pulses are averaged). Pulse averaging is accomplished by an ‘integrator, composed a 
ae _ of a capacitor in series, with a resistance. With such an integrator it is obvious 
ed wd that the longer the time constant, _ the closer will be the registered count rate col 
to the | average count rate. - Similarly, a short time constant will not. give a good : 
sample of the average pulse rate because of the statistical variation of the > 
radiation within the averaging period. If a ratemeter is subjected to an instanta- 
_ neous change from one average count rate to another, the RC circuit will not 
permit it to react immediately. Instead, the change from the first to the second 
rate will occur exponentially with time. In borehole logging the time ‘constant 


associated with a moving sonde has the effect of slurring lithologic a 


_ To compensate for this slurring effect, the digital recording systems utilized — 
in the described studies incorporate a 0.5-sec time constant. This short time 
: constant reduces the effects of contact slurring but introduces statistical errors 
which are most pronounced at low count rates. To compensate for statistical — 
errors, the digital data are subjected to a mathematical filter during computer a 
processing. This procedure reduces the slurring effect, and all indications are 
the integrity of the datais maintained. 
_ The radiation group of iogs consists of natural gamma (Nat- y), gamma- -gamma 
neutron-gamma (n-y), and neutron-epithermal neutron (n-n) logs. Of these, 
the natural gamma log is a passive radiation log, while the others are induced o/ 


‘Therefore, in most environments, the Shows a positive response to clay 


content in the formation. By ‘contrast, the | natural gamma log shows a a low = 


count rate in clean quartz sands. Because of such an intimate relationship with | 
soil and rock chemistry, the log is used extensively in stratigraphic correlations. 
"4 _ The gamma-gamma (y-y) log records the backscattered gamma radiation induced ai 
_ by a cobalt-60 source and detected by a thallium-activated sodium iodide crystal. b 4 
In the process, gamma photons are degraded in proportion to the number of i 
electrons i in the interrogated field. Since electron density is proportional to bulk - . 
~ density when properly calibrated for borehole effects, the log can be used to 
_ derive bulk density. The principle of operation is very similar to that of nuclear | 
_ density devices used for earthwork inspection. Cobalt-60 has a relatively short 
half-life (5.27 yr). Therefore, the ‘decay characteristics must be accounted for 2 
_ The neutron-gamma (n-y) log also records gamma radiation, like the gamma- — 
gamma log, but the source mechanism involved is different. As neutrons emitted ag 
by a Pu-Be or Am-Be source are slowed down to thermal energy levels, they — n 
are captured by the nuclei of atoms in the material being penetrated. At collision, ~ 
gamma photons are released by the capturing atom. These photons have energies _ 
“characteristic of the element. Water molecules are the most effective neutron 


"moderators in a natural system. ‘Thus, gamma photon production becomes a 


function of water content. However, the log i is not always a reliable ramgeredl 


| 
| 
ang The device used for this purpose responds to the potassium 
#0. (K-40) found in most soils and rock, but can also be influenced by small 
ate af _ A thasium Cleve two theese hie 
| 


_ neutron-gamma log is diagnostic of stratigraphy because of chemical anomalies — 
associated with some formations. _ orn & 
_ The neutron-epithermal neutron (n-n) log is the most reliable moisture indicator — 
of all radiation logs. The source is the same as that of the neutron- a 
log, but the detection mechanism is different. Fast neutrons, reduced to epithermal 
energy levels by interaction with materials surrounding the logging device, are 
; detected by a europium-activated lithium iodide crystal, enriched in lithium (3). 
= since hydrogen atoms are the most effective neutron ‘moderators, an 
log is indicative of the water content of the tested materials. In addition, below 
the water table the log becomes an indicator of porosity. This log is not affected 
by formation chemistry and is used most often for determination of moisture 
content and porosity. The sources used in neutron devices have long half-lives, 
thus frequent calibration for source decay is not essential. 2 Si = 
Electric Logs.—In uncased, fluid-filled holes, valuable data can be Lomond 
_ by electric logs which measure spontaneous potential and single point resistance. 
The spontaneous potential log provides a continuous record of natural potentials 
A developed between the downhole sonde and an electrode at the ground surface. Bi 
The single point resistance log records the apparent formation resistance between _ bs 
a close-coupled lead plummet and the shell of the measuring sonde. Both logs — 
are basically indicators of lithology but are also affected by water content. — 
They are best used in engineering studies where water tables are shallow, such © 
Although the electric logging devices can be run in fluid-filled holes above 
a be water table, logs so developed are difficult to > interpret and do not a 


limited Any in geotechnical investigations, and they are not discussed further 
Logging Equipment and Procedures.—The geophysical systems used in the © 
studies described herein were custom-built in university facilities for the express ah 
_ purpose of gathering engineering and geohydrological data. The tools (sondes) 
are less than 7 ft (2.1 m) in length and range from 1.0-2.25 in. (25 mm-S7 © 

- mm) in diameter. The larger diameters ae associated with the most sensitive 
radiation devices. Tools are fabricated from stainless steel and can be used 

_ inany environment where temperatures are below 90° C. Although the geophysical 2 
systems used for the studies described herein were mounted in a two ton truck, | 2 
i: operations have been successfully conducted from the back of a pickup truck — 

or jeep using a portable winching system and generator. 
_ Downhole intelligence from the sondes is transmitted to the electric modules — 


4 generates a a hard copy output for immediate review, while the other — 


| at the ground surface where the signals are bifurcated; one portion of the signal | ; 


downhole data is ‘normally accomplished at 1 ft (0.3 m) intervals, but sampling 
- intervals could be decreased to 0.5 or 0.25 ft (0.15 or 0.075 m) intervals if 
“the additional data warranted the extra processing and storage. In most cases, | 
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se to Stratigraphy. —The primary “geotechnical application 


vo 


| 
: 


of borehole been in in stratigraphic and correlation 


_ parameters, subtle lithologic changes within a stratigraphic section can often _ 


be delineated by geophysical techniques. Such information is invaluable for 
Clear definition of stratigraphic and of a number of 


= geotechnical intelligence represents a response to chemical and physical = 


“faulting i in near-surface strata. Such information can also be used in ‘constructing 

or seepage models of stratified deposits. 

“2 For stratigraphic correlations, the logs provide detailed information ms 
a _top and bottom contacts, thicknesses, and lateral continuity of diagnostic lithologic 
= units. Such diagnostic units must have one or more geophysical responses which 


iv wis 


unique and lend the unit a distinctive “‘signature.’’ Individual stratigraphic 


_ units or sequences which are less diagnostic can often be correlated on the “4 
basis of gross similarity of response, but such determinations are less certain _ 
should be supported by other geologic evidence. 


In establishing stratigraphic correlations, some knowledge | of the depositional os 


environment is essential for proper interpretation of logs. The stratigraphic — 


_ characteristics will be different for stream, lake, or wind deposits. In the writers’ 
good stratigraphic correlations were most often obtained in layered 
fluvial or lacustrine deposits. 
_ An example of excellent stratigraphic cheesieitelin! is s the site of the Palo Verde 
Nuclear Generating Station (PVNGS), located 60 mile = 


sate DESCRIPTION a= 


SuLTY CLaT (CL) piasts city 
hard, ell cemented 


astic 


caliche nodules, lens of CLATEY 


CL) at 47” 


CLAY (CL CH) medium to high 


(he 


sano (sa) very ‘very Genes, wit comenied 


Description and Natural Gamma Log (PVNGS U5- 


| 


west of Phoenix, Arizona. The site is underlain by 300-400 ft (90-120 m) of 
é alluvial and lacustrine sediments, consisting of alternating layers of clays, silts, a 
i and sands. Stratification at the site approximately parallels the existing ground 
that has a a southerly gradient of about one percent. 3 
_ The most diagnostic and continuous soil unit beneath the site is the Palo | 
Verde clay, a deposit of hard plastic-like clay interbedded with minor layers 
ing 
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. FIG. 3.—Stretigraphic Correlations Between Two Borings at PVNGS Using Natural : 


NORTH- SOUTH PROFILE 


ole 


FIG. Across the PYNGS Site Using Natural Gamma 
of silts and sands. Fig. 2 displays the lithologic and natural gamma log of a © 
za continuously sampled boring at the center of Unit 5. The natural gamma 


he Palo Neste clay is traceable a as a continuous deposit not only within 
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thes Fe. 3), but also over the entire sit 
extending several ‘miles (Fig. 4). In the power | block area, area, subunits within and — 


the entire site (Fig. 4). ‘In both profiles, the natural gamma log has been used 
for correlations. Similar results were obtained with other geophysical logs. _ ae 
At this site, ‘the us use e of geophysical | logs was essential in distinguishing 


block a areas, however, where boring spacing was amch closer, adequate strati- 
a graphic correlations could have been es established on } the | basis of of “lithology : and 
Although individual logs have been discussed as valuable instruments for 
roa 4 correlation, it must be stressed | that that the use of a a single log can | be grossly | 


misleading and can an often lead to ) erroneous findings. All logs comprising ‘ the — 


vat SOUL DESCRIPTION, NATURAL GAMMA LOG 
i (SP) with lenses of gravel and q ai 
edt wi bres A od AA. doderod 


ayers and silty sang 


oc 
matrix 


= Angular rock fragments and gravel 
4 al origin in a cemented sand 


CLAY (CH) fine 
d 

a FIG. 5.—Soil Description and Natural Gamma Log of Sundesert Boring B-DH- 35 

guise > from a given boring s should be considered jointly i in the the process of of f composite 
= - mutually supportive of the same samnietien is it likely to be the correct one. | 

a Quite often, however, one log will be more diagnostic than the others and 

will become the basic correlation took, 

- -Anexample of a site where the use of multiple logs was essential for stratigraphic 
correlations is the site of the Proposed Sundesert Nuclear Plant near Blythe, 

California. The subsurface profile consists of approximately 330 ft (100 m) of 
- Colorado River sands, underlain by a fanglomerate (cemented rock fragments — 
of local origin) and the lacustrine Bouse formation (Fig. 5). Of interest at this “ 

I ti - site is the positive response of the natural gamma log to both clays and local 

rock or gravels. Natural gamma activity alone did not, therefore, serve to identify — 
these stratigraphic units. It became necessary at this site to define eochysicel 
units in terms of porosity plus natural gamma activity (Fig. 1), or the vertical — 
distribution of natural gamma activity, or both, within a single unit. Fortunately, 


| 
@ ove the oO Verde clay are also traceable across borings spaced about % 
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the stratigraphic the ‘formation (Fig. 5) displayed a 
cole amount of continuity which facilitated the correlations. e'grhostee — 
Geophysical Response to Physical Properties.—The induced radiation logs can can 3 
also be used in quantitative or semi-quantitative fashion to determine the physical _ 
_ properties of the materials found in the borehole (1). To obtain quantitative ps 
» data from the induced radiation logs, it is imperative to compensate for the ond 
influence of all factors other than the physical property being measured. Such — 
_ 3. Properties of borehole fluid (if any). | tat 
Most of these factors can be filtered out of the record or compensated for 
by careful calibration procedures conducted in pits constructed from materials _ 
A more but potentially le less calibration is on- n-site 


a and the physical properties determined by laboratory tests. The correlations a 
are based on the following two empirical relationships s between recorded count 
rates and physical properties being measured: ieee 
— 


‘in which N. = gamma- gamma count rate; N,, = neutron-neutron count rate; 
Y, = total (bulk) density; WwW, = = water content ‘per total volume (total porosity 
; below the water table); and a, b, c, dys = = constants established by calibration. 
_ The calibration constants can be determined by selecting at least two sets 7 
of data points for each relationship. Other physical properties, including water — 
content (by weight), total or partial porosity, void ratio, dry density, and degree 
_ Of saturation, can be indirectly derived from the above relationships if the specific 
gravity of the earth materials and pore fluid is known. 0 POTOSI 
_ The induced radiation logs are very sensitive to variations in borehole size. . 
_ Therefore, it is imperative that a caliper log be obtained to determine the variation — 
of borehole diameter. Zones with significant variations in borehole diameter 
should be excluded from quantitative analyses. evade cA 
Induced radiation logs have been used by the writers in several siting studies 
to supplement the stratigraphic correlations generated by the natural gamma 
log. No express attempt was made during these studies to obtain geophysical © 
logs for quantitative assessment of physical properties. However, review of 
_ the data obtained indicates that the logs contained information that could have 
_ been used for quantitative or semi-quantitative evaluations of physical properties. 


At the site of the Palo Verde Nuclear Generating Station, Subtle differences 
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im 
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the neutron-epithermal neutron (n-n) and neutron-gamma (n-y) logs together 


pes “BOREHOLE GEOPHYSICS = 
with | porosity and moisture. content. The response of the Palo’ Verde 
clay is clearly recognizable on both the n-n and the n-y logs but is more re 
on the n-n log. The difference in response on the n-y log is believed to | be 
due to the influence of chloride ions in the soil. 
a Attention | is also directed to the neutron log baseline shift associated a 
the Palo Verde clay section (depth range 200-280 ft). This shift is not as 0s men 
in laboratory porosity and moisture data and is illustrative of an important 
characteristic of the neutron response. The neutron tool responds to both 
ree- -water molecules and chemically- combined water molecules. In keeping 


will display ‘neutron moisture c or porosity responses. Therefore, two neutron | . 
calibrations would be required in boring U5-B43, one for the Palo Verde clay 
and one forthe remainder of the section 
Fig. 7 displays quantitative correlations between geophysical logs obtained a 
in Boring US- B43 and physical Properties determined ad acaba tests on 


WATER CONTENT 


»pither- 

a neutron (n-n) log with an overlay of porosity data, exponentiated in accordance > 

with Eq. 3. The solid curve shows the log calibrated for the portion of the 

= above the Palo Verde clay. The dashed curve represents a separate 
calibration | for the Palo Verde clay that removes the baseline shift associated as 


with the variation in hydrous clay y minerals. 
Plot (6) presents bulk density data from tests plotted over the 
gamma-gamma (-y-y) log, exponentiated in accordance with Eq. 2. The log still 
contains the effect of all borehole irregularities. Plot (c) presents a plot of 

a dry density data from laboratory tests over a hybrid log derived from the ff 

combination of the n-n log and the y-y log. The effect of moisture content, 


represented by by the log, has from the bulk density, represented 
by the y-y log, to produce the dry density log. 


= The three overlays show excellent correspondence between | physical parameters 
determined by tests: and top scaled to ‘Span the range 


¥ 
| 
6.—Comparison of Water Content and 
Porosity Data (PVNGS Boring U5-B-43) q 
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of Such correspondence i is remarkable consider 


ing the following: 


1. The geophysical logs represent 8 a continuous record of responses while a. 
“the data are determined on of finite ver (typically 
The effects of borehole size. variations (agoeity) reflected. on the caliper 


log (Fig. 7) have not been removed from the record. oe ca, as 


Extreme variations in physical parameters were faithfully detected by geo- 
physical logs in a cored boring at Cheney, Washington. The boring penetrates 
a sequence of basalts and | interbedded sediments of the Columbia River Basalt 
_ Group. As might be expected, the basalts exhibited much higher density and — 
much lower porosity than the sediments. At the time of the logging, the boring — 
diameter was 3 7/8 in. (NC size) to about 600 ft (183 m) and 2 63/64 
in. (NX size) to 1,174 ft (358 m). The boring was relatively true to gauge and 
breakouts were minor, so the only borehole effect was the change in diameter 
Correlations between the gamma- a-gamma log and bulk density, and between _ 
_ the neutron-epithermal neutron log and porosity are presented in Fig. 8. The § 
_ close correspondence between laboratory-derived information and the two 
geophysical logs is evident despite the extreme range in physical parameters — 
. displayed. The effect of the change in borehole diameter can be seen in both — 
logs, but it is much more pronounced in the gamma-gamma log. This significant _ 
effect would necessitate a separate calibration below 600 ft (183 m) for the ~ 
a S. while a single calibration appears to to be adequate 1 for the porosity 
The preceding examples demonstrate that, with careful calibration, — 
determinations of porosity, density, and related parameters can be made with 
- reasonable accuracy using geophysical logs. Once the logs are calibrated for af 
typical boring at the site, correlations can be extrapolated 
7 as long as physical and lithologic conditions remain vasa 


AND 


tions has been demonstrated by several examples from practice. Based on the 
studies described herein, the following can be concluded. 
__ 1. Geophysical borehole logging is a useful tool for stratigraphic | characteriza- 
tion and correlation, particularly over large distances. pa a) "ee 
2. Quantitative determinations of physical parameters sucha as moisture content, 
density, porosity, etc. can also be made using geophysical logs. At present, 7 


such determinations require careful on-site calibration to compensate for borehole 
effects, decay characteristics of the radioactive Sources used in the tools, and ~ 


| 
| 

i 

The methodology described herein has been us ry for 

_& over half a century. Tools used for such applications are more sophisticated | 

and much costlier to operate than those used by the writers. » . ee 

= 


At present, geophysical logs “of limited value for small-scale 
geotechnical investigations, consisting of a few, shallow borings. However, as’ 

. « number and depth of borings is increased, savings can be realized by ee 


program and substituting geophysical exploration. 


The writers would like to thank Fugro, Inc., for its support in the iia 
of this paper and Eileen Poeter for providing the data from the Cheney, Washington 
a The rest of the data reported in this paper were generated during siting 

investigations for two nuclear power plants (Sundesert and Palo Verde) and F 
ss part of the public record. The writers wish to thank San Diego Gas and © 
Electric Company and Arizona Public Company for permitting the the use 


these data for this 
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and Joseph R. Curro, 


“The of the dynamic response of a soil or structure 
on soil, requires” knowledge of the stress-strain properties of the foundation 
material. When using current computer analysis techniques, the initial elastic 
moduli low strains, generally cm/cm (or less) are 
the variation in moduli asa function of strain. | 
_ These initial moduli (shear and Young’ s) can be determined in situ using — 
bs various geophysical techniques. Or ‘resonant column tests can be conducted 
in the laboratory on undisturbed samples to obtain moduli as functions of strain 
i. Since the moduli previously mentioned would be needed for dynamic progam ; 


y design of four proposed buildings and a test loop structure (13), a test program 
was formulated that consisted of field and laboratory investigations to provide 
required moduli. Ih addition, -Poisson’s ratios were determined from the 
_ investigations to aid in the dynamic. design. 
_ Therefore, the primary purpose of this project was to determine values of — 
shear and Young’s moduli for the foundation materials underlying the five _ 
— structures (designated as X-340-1 through X-340-4 and X-774), in Fig. r 
. Additionally, building X-340-4 was arbitrarily selected as an example case a 
« illustrate the procedure used for determining dynamic. design values. The 
7 ‘purpose of this paper is to compare dynamic soil properties determined i in situ a 
and in the laboratory and to examine the rationale used in obtaining design 
‘ Site Description.—The construction site is located in the vicinity of Portsmouth, a _ 
"Research Civ. Engr., Geotechnical Laboratory, U.S. Army Engineer Waterways 


Note.—Discussion open until March 1, 1982. Separate discussions should be submitted 
for the individual papers in this symposium. To extend the closing date one month, 7 
4 written request must be filed with the Manager of Technical and Professional Publications, F 
ASCE. Manuscript was submitted for review for possible publication on February 17, 
ont This paper is part of the Journal of the Geotechnical Engineering Division, Proceedings Ste 
of the American Society of Civil Engineers, © ASCE, Vol. 107, No. GT10, October, 
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at the site were primarily Pleistocene river-bed and lake- bed S decetin. ‘Fig. 
isa geologic map of the site (13). The map outlines the areal distribution of . 
the glacial lake-bed deposits and early alluvial deposits. The lake-bed deposits, . 


- consisting of fat, thinly laminated clays, are largely confined to the old Scioto 


- The area outlined as lake d deposits in Fig. 1 shows the occurrence of a 
la 


laminated clay either at the surface or at some depth beneath the surface. 
Lean clays, silts (alluvium), and a fill material, which was placed during the " 
period 1953-1955, lie above the fat clay in various thicknesses over a large a 
portion of the outlined area. Early river alluvial silty th and oh occupy sla 


portion of the valley fill outside the fat clay boundaries. i 


th 
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+ ELEVATIONS ON BASE OF FAT CLAY, MSL 

SHALE, MISSISSIPPIAN (CUYAHOGA) 
ALLUVIAL DEPOSITS OF LEAN 

AND SILTS (PLEISTOCENE) 

LAKE BED DEPOSITS (FAT 
CLAY) AT OR BENEATH THE SURFACE 


wie FIG. 1 .—Proposed Building Locations and Geologic Map of the Site one é 

Up to 60 ft (18 m) leveraging 35-40 ft (10. 5-12.0 m)] of early river-bed z= 
and lake-bed deposits in the old valley still remain at the site. These deposits — 
_ are mostly fat, laminated clays, lean clays, and fine silts, with a2-to 10-ft 

(0.6- to 3-m) thick layer of weathered clayey gravel at the base. l-to 10-ft 
60. S to 3-m) thick layer of residual soil and colluvium mantles a large part 


7 of the valley. Rock, consisting, of shale and sandstone, underlies the soil (13). 


‘The field investigation w was s comprised of two o phases: (1) Drilling : and ‘sampling; 

(2) geophysical investigation, 
ie Drilling and Sampling.—During this phase of the field investigation, 27 holes . 
were drilled at the site for ‘in situ seismic test use and to obtain undisturbed 
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cag 


>. 
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one set near each end of the four sapatad eciaiins locations and one se 
at the test loop location (X-774). The three holes were oriented to form a 
L-type configuration with the distance from the boring at the vertex of the 
& Lto the boring at the end of each leg being approximately 20 ft (6 m). Undisturbed _ 
samples were obtained from the boring at the vertex of the L for each set 
with a 5-in. (12.7-cm) Hvorslev fixed piston sampler. The remaining holes were 
drilled with a Waterways Experiment Station (WES) modified fishtail bit with 
upward baffles. After the holes had been drilled to ) the desired depth 
ft (0.6-2 m) into rock], 3-in. 


coment Test 


NOTE THE ELEVATION (TOP OF BORING 821 
WAS 67) FT. wiTh THIS 


THE READER 
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«FIG. 2.—In Situ Seismic Test Layout; ; Proposed Building X-340-4 Location ay 


FIG. 3.—P- and S-Wave Velocities Determined from Crosshole Tests; ; Proposed ~ 


4 Building X-340-4 Location mow diged amit levind—2 

_ installed in the boreholes, and the annular space between the casing and walls " 

% of the borings was grouted with a mixture of cement, bentonite, and water 
that, after setting up, had the approximate consistency (modulus) of soil, = = 


_ The geophysical investigation conducted at the site consisted of nine series 


of in situ seismic tests. Each series was comprised of crosshole, ‘downhole, 


| 
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Crosshole Test Procedure. _—The crosshole tests were using a portable 
- 24-channel refraction n seismograph. A recorder displayed the data while operating 
at a paper speed of approx 35 in./sec (89 cm/s). Timing lines were od 
on the oscillogram at 10-msec intervals. Resolution time with this equipment 
and the above oscillograph speed was 0.5 msec. The test procedure to obtain 
- shear wave (S-wave) data was a modified technique that used the same vibratory 
_ source and control package as that described by Ballard (1). The primary difference * 
is the two procedures was concerned with the data acquisition package, i.e., 
~ no signal enhancement was used. Rather, the vibratory source was swept through 


Incremental velocities 


FIG. 4 4.—Arrival Time Versus Depth from Downhole 


msec 


5.—Arrival Time Versus Depth from Downhole Tost; | 824; Proposed 
a frequency range while monitoring the output of the two geophones placed — 
at the same elevation in the receiver boreholes. (The vibrator was producing — 
vertical oscillations, thereby transmitting a vertically polarized shear wave.) 
When an acceptable response was received at a specific frequency, this frequency — 
was then interrupted by a tone-burst generator to send a specific number of 
cycles of that frequency to the receiver units. The source geophone was displayed — 7 
along with both receiver r geophones on In 


Proposed 
! 
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“time difference ‘that and the signal at each of the two 


The tests to obtain S-wave data were conducted by first placing the source | 


and receivers at the same elevation near the top of the boreholes, then pulsing 
a the source unit several times, recording both transmitted and received signals. 
_ After a satisfactory record had been obtained, the units were then repositioned 


9 5 ft (1.5 m) deeper. The procedure above was repeated at this and each succeeding» 


= 1545 


FIG. 6.—P-Wave Arrival Time Versus Diste Distance; Prope X-340- 4 ‘Location; 


an 7.—P- Wave Arrival ‘Time Versus Distance; Proposed ris X-340-4 Location; 


“within +3 in. (47. 6 cm), thereby, eliminating the need for a formal ‘svaton 
The compression wave (P-wave) data were re obtained in a similar manner, 
except that the vibrator and its associated instrumentation were not used. Instead, 
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relatively shallow, their verticality could be checked by visual observation 


bridgewire ¢ detonators (EBW's)w were ere used at 1.5- -m) 


crosshole tests were the times required for P- and — 
_ S-waves to propagate from the source to a point of detection. These times 


— 
Borings 82k, 825,826 rs Dorings 821, 822, 823 


All velocities aps asterisk 


denotes velocities from crosshole tests. 
One asterisk denotes velocities from 
a downhole tests. Tvo asterisks denote 
velocities from refraction seismic tests — 


FIG. 8. —s- Wave Velocity Profile from in Situ er ae : Proposed Building X- 340- 4 


Borings 621, 622, 823 


q FIG. 9.—-S-Wave Velocity Profile from In Situ Seismic Tests; ; Proposed | Building X-340-4 
t 

é interpretation, based on Snell’s law of refraction, is used to determine true a 

velocities by accounting for zones of high velocity contrast. 
“ Downhole Test Procedure.— This test, designed to to 0 provide data for determination 7 


seismograph and ‘recorder as used for the crosshole tests. 

‘The receiver geophone was placed at a depth of 5 ft (1.5 m) in a borehole, 
and a vertical hammer blow to a steel plate positioned approx | ft (0.3 m) £ 
from the mouth of the boring on the ground surface was used as the P-wave — 
‘source. A geophone adjacent to the hammer impact point provides the time 


 &§ 
_were then divided into the distance between source and receiver geophones 
: 7 to provide apparent velocities. If a nearby higher velocity layer exists, the Jf 
: ] wave may refract and travel along that layer, th " 
— path than the direct distance path. A computer prog} 
= 
= 
3 
© Seismic source and geophone denotes velocities 
\ | 
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SOIL MODULI DETERMINATION 
Of seismic input. S-wave tests were then conducted by placing a large wiles > 


plank on the ground surface and striking the plank at each end, thereby, reversing 
_ polarity of the horizontally polarized S-wave, which facilitates identification. 
The procedures for obtaining P- and S-wave data were then repeated with the 
receiver geophone at 5-ft (1.5-m) deeper increments until the bottom of the 


Surface Refraction ‘ Seismic Test Procedure. —Conventional surface refraction 
seismic tests were conducted at the site using the 24-channel unit. Refraction — 

- lines run at the site were 240 ft (73 m) in length. Geophones were idaaet 
at 108 (3-m) intervals with the first geophone positioned 5 ft (1.5 m) from © 
an explosive source. Both forward and reverse traverses were made at the _ 
‘- site so that true velocities, in addition to apparent velocities, could be Meplne glk 

depths to refracting interfaces computed. ps sie ni 
Surface Vibratory Test Procedure.—In order to provide a necessary input 
parameter (near-surface S-wave velocity) to the crosshole computer program, 
‘ surface vibratory tests (5) were conducted near each crosshole set. Rayleigh 

(R) waves may be generated using a controlled energy source such as a vibrator. 

_ In this case, the electromagnetic vibrator used to conduct the cross-hole tests 
was also" as the R-wave so The vibrator at 


FIG. 10.—Shear and Young’ s Moduli Profile oon In ‘al Seismic Tests; Lal 


f frequencies, and the R-wave was monitored by the geophones used for the 


_ surface refraction seismic tests, , which were placed in a line at intervals o q 
1 ft (0.3 m). In this manner, the velocity of the R-wave can be determined 
from: a plot of R-wave phase arrival time versus distance for each "specific 

¥ frequency. Then, a corresponding wavelength can be calculated by dividing ay, 
- the frequency into the velocity. Wave velocities thus derived are considered 

, to be average values for an effective depth of one half the wavelength (10). 


 S-waves, for practical purposes, can be considered to have the same velocity © 
as 


va Since the vibrator | used at this site was of | limited force output 50 Ib (222, 
_N), the depth of investigation was limited to a approx 5 ft (1.5 m). These data 


could then be compared with ‘those obtained for the shallow crosshole- test 


/ 
q 
| 
q 
| 
procedures should be converted to elastic parameters, such as shear modulus, 


G, Young’s E, and Poisson’ s ratio, use in design. 


; This conversion can n be accomplished using theory of wave propagation in eetie 


media (8) and the total mass density, p, of the medium. Yo 
ew sted! stew ave @ 


_ The results determined from the in-situ seismic tests conducted at the proposed _ 
location of building X-340-4 are presented and examined below. It appeared 7 
_ feasible to the authors to select one location as “‘typical”’ and to to examine the i. 


test results obtained ithere, cout ‘ody f not 


Proroseo Buioinc X-340-4 
Two series of in situ seismic tests were conducted at this location, « one near 


end of the proposed building location (Fig. 2), 
TABLE 1.—Samples Obtained in Boring 821 UD from Which Resonant Column +; - 


Depth of layer PL. consolidation 
Deptt the | stress, in pounds 
Sample | in feet 


A 
represents, in me per square inch 


number | (meters) 


(27-42 
1.3) 


| 14-22 


(61-66) 43-6. (41.4) 

24.5-26.6 9.5 

9.3) 


, and 823 couprieed one set near south end of the building, by 
_ whereas the second set, near the north end of the building, was made up of 
borings 824, 825, and 826 (Fig.2). 
Borings 821 and 824 were used as the seismic source boreholes for yr their 
- respective sets, and borings 822, 823, 825, and 826 were employed as receiver — 
holes for their respective sets. The seismic source and receiver (geophone) — 
placements used in the performance of the crosshole tests are shown in Fig. 
_ 3. True P- and S-wave velocities determined from the crosshole data are presented 
a alongside the | receiver and source locations. Note that only one oe and S- wave 


were determined from the seismic source borehole to each of the receiver 
boreholes for each crosshole set. In comparing the velocities from the source 
to each receiver, the spread was so narrow, less than 10%, that it was thought 
a average the tw two | P-wave velocities am and, ,the two “wav ye 
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downhole tests located i in Fig. 2¥ were conducted, one 
using boring 821 and the other using boring 824. The results of the downhole 
A tests conducted in borings 821 and 824 are presented in Figs. 4 and 5, respectively. = j 
As shown, average and P- and S- -wave were determined 
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from the downhole em. 
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FIG. 12.—Damping Ratio Versus Dynamic Resonant Column 


Tests; 821 UD No. 2-1; o,., = 7 kPa; B = 0.18; Nonsaturated 


_ Surface Refraction Seismic Tests.— te tests consisted of four traverses (two: a 
‘Fig. 2. The time versus distance plots for these traverses, designated as S-15- S- ie 
3 shown in Fig. 6 and 7. These figures also present naa a true 


| 
11.—Dyna 
Sample 


Surface Vibratory Tests. —One surface vibratory test was conducted at 
- 821 and another at boring 824. The results of these tests yielded average R-wave 


velocities of 375 and 235 _ (114 and 72 m/s) for the near- r-surface ‘material, 


{ FIG. 13 Versus Dynamic Strain from Resonant Column 7 
Tests; Sample: 821 UD No. 9-1; 41. 4 kPa; B 0.68 ‘ bolt 
ate 


iL 
a 


STRAIN «,, RAD/RAD, 


were used as true S-wave vel: the 

and S-Wave Velocity Zoning.—The P- and S-wave velocity results determined 
from the crosshole, downhole, surface refraction seismic, and surface vibratory 


14.—Damping Ratio Versus Dyharm Resonant Column 


Sol L ‘MODULI DETERMINA 


oy + 
at ‘this. building location were analyzed and to 
P- and S-wave velocity profiles shown in Figs. 8 and The 


due to’ a large increase in -wave velocity as percent saturation of the soll 
reaches 99-100%, whereas S-wave velocity isnot affected. = = 
4 The Pweve is very percent 


reason, the 
— 


FIG. 15.—Dynamic Modulus Versus Dynamic Strain Amplitude from Resonant Column 7 : 


Yom; Sample: UD No. 11-1; 0,,,=65.5kPa;B=083 


r 


DAMPING RATIO, PERCENT 


payee 
esis 


Tests; Sample: 821 UD No. 11-1; o,., = 65.5 kPa; B = 0. oo 120 lalgine? jeans 
“Referring to Fig. 8, five velocity z 

data. The velocities presented for the z zones were e usually weighted te onl 

based on data quality. A 1,400-fps (425-m/s) velocity was established for the 
near- -surface zone, which sen eee about 7 ft (2.1 m) in thickness. The et 


4 | 
.. 
q 


‘ ranged from 5-13 ft (1.5-4 m) in thickness and is best delineated by the ——— ; 


seismic data. The next two velocity zones [7,100 and 9,400 fps (2,165 and 
2,865 m/s)] could be collapsed into one. However, since the top of the 9,400-fps i 
Qs 865- —-- zone coincides with the top of shale, it was believed angie 


= LEGEND 


Lone: 


FIG. 1 —Dynemioh Modulus Versus Strain Amplitude from Resonant Column 
Tests; Sample: 821 UD No. 13-1;0,,=74kPa;B=091 


rer 


FIG. 18. —Damping | Ratio Versus Dynamic Strain Amplitude from ee Column 

‘Tests; Sample: 821 UD No. 13-1; o,., = 74 kPa; B = 0. 

for analysis purposes to define top of rock and to keep the zones separate. 

The 7 ,100-fps (2,165-m/s) zone interpreted from the crosshole and downhole | 

tests in borings 824-826 is probably indicative of a silt, sand, and gravel mixture 
cemented to form the conglomerate noted in the laboratory classifications of 4 


| 
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(Of particular interest is the data ‘obtained for the third rom 
‘surface (Fig. 8). This layer is below the groundwater table, yet it has a P-wave 7 

Velocity of 2,800 fps (850 m/s). One would expect a minimum P-wave velocity 

g approx 5,000 fps (1,530 m/s) which is the P-wave velocity of water. Since | q 
such a velocity was not obtained, _ the layer m must be less than fully saturated. hs 


saturation is in 1 the range of 98- -100%. It is quite possible, « even likely, that “ 
the material in situ contains entrapped air or gas. For whatever reason, the 


indicates that this zone is not fully saturated. 
aie to Fig. 9, five S-wave velocity zones were interpreted from he 


4 a 500-fps (150-m / s) average a thickness of about , 15 ft (4. 6 m). The third 
- zone averaged 700 fps (215 m/s) and was approximately 12 ft (2.7 m) thick. 


"The two deeper zones [1,750 and 2,000 fps (535 and 610 m/s)] could probably 
samples from various deaths borings UD, 
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if 4 


FIG. 19.—Dynamic Shear Modulus Versus Dynamic Shear Strain Amplitude wom 

have been combined; however, for the same reason given for the P-wave zoning, 

they were not. It will be noted that the value of wet unit weight applicable 

for each zone is also presented in Fig. 9. These values were obtained from 

the undisturbed samples used in the resonant column tests and from laboratory 

tests conducted on samples from borings near the in situ seismic test locations 


_ Shear and Young’s Moduli and Poisson’s Ratio (G, £,v).—The interpretation 7 


of the data was continued by establishing zones from the velocity profiles and - 
determining values of G, E, and v for these zones. This was accomplished . 
by overlaying the P- and S-wave velocity profiles (Figs. 8 and 9). This overlay 
produced seven zones. Then, values of G, E, and v were computed for each 
zone using the P- and S-wave » velocities and wet unit weights that fell into q 
each zone. After analyzing G, E, and v values for the various zones, it was 2 : 
evident that the controlling factor was the S- “wave because the 


— 
4 
= 


P-wave velocity variations hed ‘little “effect on v course, 
‘ no effect on G values. Therefore, the seven newly established zones could ee 
be reduced to the five original zones shown on the S-wave profile (Fig. 9). *, 
The final field interpretation for this building location showing the five zones ‘4 
with their associated G, E, and v values is presented in Fig. 10. 
Undisturbed samples taken in the field were bela by truck to the WES ‘et 
pen in the WES laboratory, the samples were opened and classified. The 
ale _ samples were evaluated to see which were representative of field conditions ee 
Bema: te the boring logs a and general knowledge of the site (13) and which 3 


‘Resodent column tests were performed on the representative samples found 
be suitabie for testing, which did particles of gravel 


The equipment used in this series of resonant column tests is ne as the j 
a Long-Tor resonant column apparatus. Resonant column apparatus >a 
normally excite the specimen only torsionally; however, this equipment differs . 
in that it also includes an electromagnetic oscillator for longitudinal vibration — We 
: Ec of the specimen. A more detailed description of this ee is given elsewhere 
» 


‘Test Proceounes 
all tests were performed on 2. sin in. diam x x 1- -in. long (7.1-cm x 17. 
“ undisturbed specimens. These specimens were ‘isotropically consolidated | so as” i 
Pa to correspond to the in situ mean principal effective stress or octahedral stress, __ 
— &.... The in situ mean principal effective stress was computed assuming normally — 
"consolidated conditions (13) and a coefficient of earth pressure at rest of 0.5. ee 
hee primary consolidation was complete, the drainage valve was closed 


to sinusoidal excitation in both the longitudinal and torsional modes. Longitudinal 

and torsional vibrations were not conducted simultaneously. The frequency of 
* vibration was varied until resonance was obtained and data were recorded such a i 
that the modulus and strain amplitude values could be computed. In order to 
obtain the Variation of modulus as a function of strain amplitude, the 


and the B valve was checked and recorded (3). The specimen was subjected r. 


Consolidation. of all 1 resonant column tests was started at the end of the 
day and specimens were to general, primary 


to sinusoidal excitation the morning. testing required approx 
_ 2h. Thus, the data are plotted for a time of approx 12h after primary consolidation asa J 


was complete. Evaluation of the effects of time on the modulus and damping 


| 
: 
ne SOU Was beyond scope 0 nis study. Is recognize a ime errects 
es may be significant for the clay CH materials, = 


SOIL MODULI DETERMINATION 
Lin 
A At the initiation of ‘this investigation, it was decided that the date obtained 


during torsional excitation were of primary importance. These resonant column 
tests were not | conducted on back ‘Pressure satu saturated specimens 


saturation; and @ the equipment does nc not lend itself to easy back pressure . 
j saturation. The top cap has no provision for vacuum or drainage which is needed ; 
for seepage saturation and an air fluid interface is required in the chamber; oa 
consequently, lower chamber pressures are desirable. 
_ The procedure above is a brief description of that used by the WES. A a 


_more detailed testing procedure, along with a complete examination of the analysis 
and calculations required to obtain shear and Young’: s moduli from the raw 
data, is given elsewhere (4). ‘ 


Resonant Co.umn Test Resuits 


_ Undisturbed samples from \ various depths | in if Sed 800 UD, 803 UD, 806 
UD, 809 UD, 812 UD, 815 UD, 818 UD, and 821 UD were tested in the resonant — 
column device (3). A total of 37 specimens were tested. tet teed 
_ Ref. 3 contains plots of modulus and damping versus dynamic strain amplitudes 
for all the resonant column tests. In order to familiarize the reader with the 
format of the resonant column data and its interpretation, the t test results obtained 


_ The samples taken from boring 821 UD, which were tested in the laboratory, 7 
are listed in Table 1. Also mepeeene 0 are the depths from which the = | 


‘The results of these tests are presented in Figs. 11- 18. As expected, the — 


- modulus decreases and the damping increases with increasing strain (6,7,9,12). 
4 The data in Figs. 12, 14, 16, and 18 show a trend indicating that the damping 
obtained’ in the torsional vneutens mode is on the order of of one- half that determined 
i, Below a depth of 31 ft (9.4 m), the soil is extremely gravelly. Samples obtained 
were disturbed to such an extent as to be considered for 
Rock was found at a 1 depth of about 36 ft (il ‘m). tlt 


resonant data obtained from specimens from samples 
taken in boring 821 UD have been presented. A soil profile will be developed 
for this boring i in terms of shear modulus, Young’s modulus, and internal damping 
Shear Modulus. —Fig. 19 summarizes the shear modulus data obtained from 
testing four specimens taken from samples of boring 821 UD. Also shown in 
® Fig. 19 is a soil profile depicting the four layers that were observed in the 
a boring. Note that the near-surface clay layer is 14 ft (4. 3 m) thick and has 
a shear modulus at low strain levels of about 3.5 x 10° psi (24 x 10° kPa) 
10° kPa = 1 MPa). As the shear Strain is increased to a dynamic strain level — 


- 
7 | 
| 
q 
EE 


10° ‘psi 6. 9 MPa) 
as the ‘strain level j increases to 10~° rad /rad. Below a depth of 22 ft (6.7 m) 
there exists a 4-ft (1.2-m) thick silt layer. This material exhibits a shear modulus 
of 17 x 10° psi (117 MPa) at low strain levels. As the strain increases to . 


OLSCRIPTION 


Column Tests; Boring 821 

a strain of rad /rad, dynamic dene 
decreases toa value of 5.6 x 10° " psi (38. 6 Mme. Below the 4-ft (1.2-m) thick 


and 31 ft (7.9 and 9.4 ~~ This depen exhibits a shear modulus of about 5. 6 
x 10° psi (38.6 MPa) at low strain levels. This modulus" decreases to about a 
2.8 x 10° psi (19.3 MPa) as the strain level becomes greater than 10“  rad/rad. 
In Fig. 19 several specific items should be noted. These are 


_ about 107? rad/rad, the shear modulus decreases to about 1.7 x 10° pi_—iég 
4 
= 
| 
— 
a 
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ie rom 10~°-10~* rad/rad, the shear modulus 
a is generally reduced by about 50%, oa in the silt layer at a depth of 22-26 4 
_ ft (6.7-7.9 m). In this layer, the shear modulus is reduced to about one third © 
---2. One would expect that the modulus would be a function of the mean 
F principal effective stress. Thus, the modulus should ‘increase as a function of 
; depth, other things being equal. This is not the case in boring 821 UD. The 
clay layer at a depth of 14-22 ft (4.3-6.7 m) is not stiffer than the surficial 
a 3. The 5-ft a. 5-m) thick clay layer immediately below the silt has a modulus eg 
that is significantly lower than the silt layer at a depth of 22-26 ft (6.7-7. £ 
 m). This is not unusual, but should be noted because it may be important _ ‘ 


in any future dynamic analyses or calculations. ime old 


burebiencs of b bus papal 


wol 


soupi 
‘Tests; Boring 821UD mov: bas ble no 
 Young’s Modulus.— —Fig. 20 shows Young’s modulus (£) versus dynamic strain — 
amplitude as a function of depth in boring 821 UD. The 14-ft (4.3-m) thick 
surficial layer has a Young’s modulus at low strain levels of about 26 x 10 , 
psi (179 MPa). This value decreases to approx 15 x 10° psi (103 MPa) when — 
_ the strain level increases to 5 x 10 ~* in. /in. (cm/cm). The clay layer between 
a depth of 14 and 22 ft (4.3 and 6.7 m) is similar to the surficial layer. This 
8-ft (2.4-m) thick clay layer has a Young's modulus of about 26 x 10° psi a 
(197 MPa) at a strain level of 4 x 10~’ in. /in. (em/cm). As the strain increases _ 
to a level of approx 10~* in. /in. (cm/cm), the Young’s modulus decreases 


7 | 


_ e silt layer that lies 
@a pe of 22- “26 ft 6. 7- 1, 9 m):i is approx 56 x10? psi (386 MPa) at strain 
levels below 10~° in./in. (cm/cm). As the strain level increases above 10~° 3 
in./in. (cm/cm), this Young’s modulus decreases to approx 20 x 10° psi (138 2% 
- MPa) when the strain level is about 4 x 107° in. jog remnant The 5 be (1.5-m) 


. layers. The two clay layers” above 22 ft (6.7 m) have similar characteristics, 
" with the layer located between 14 and 22 ft (4.3 and 6.7 m) being more sensitive | 
strain than the surficial layer, 

_ Damping Ratio.—Fig. 21 presents the damping ratio versus dynamic strain | 
amplitude as a function of depth for boring 821 ‘UD. In dynamic analyses and 
_ in the calculation of dynamic response of soils, ‘the predominant energy input 
comes from shear waves. Consequently, shear or torsional damping is the most 
important consideration, and only shear damping will be examined here. The od 
; Se damping data were sometimes atypical (3) and should be considered = 
suspect (see Fig. 14). Close study of Fig. 21 reveals that damping is relatively a 
_ constant as a function of depth for boring 821. For strains less than 10~* rad /rad_ 
the damping ratio is on the order of 3%. As the strains increase above 1o~*, 
then the damping ratio approaches 12%. Based on WES 
be reasonable and consistent. 


a depth of 22-26 ft (6.7-7.9 m) is stiff compared with the other 


~ column tests. In general, shear and compression moduli decrease as a function 
of i increasing strain (4,6,7) while damping increases. 
_ Fig. 22 shows a boring log for boring 821 UD along with iiiaiiniass ee | 
_ s of shear modulus, Young’ s modulus, and damping as a function of depth a 
at low (~10~° strain levels. The subdivision of the ‘soil was based 
on soil type and laboratory visual classification. 
For comparison purposes, the values that were determined in ‘the field by 
_ geophysical techniques are shown immediately to the right of the laboratory 
values. The reader is reminded that the subdivision of the soil profile is based 
on field techniques and results from changes in P- and S- -wave velocities and 
. independent of material types or visual classification. Consequently, it is” 


not surprising that the soil profile is not subdivided identically using data obtained ~ 
_ Generally, moduli values obtained in the laboratory are less than in situ 

determined values. In fact, the values obtained for laboratory tests may be 4 
as slow as one half the field values (1 1). The: fact that the ; laboratory data obtained g 


trend can un be ‘obtained from a “comparison of the data presented in | Fig. 2. < 
For some zones (soil layers) the laboratory values are approximately equal to | 
the field values. For other zones the laboratory values are greater than the 


S:—C 
x has a Young’s modulus at low strain levels of about 43 x 10° psi (296 MPa). hi ; 
Again, as the strain increases from a level of around in. /in. (em/cm), 
the Young’s moduli decrease to approx 16 x 10° psi(IIOMPa). 
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fil eld values and vice versa. Because of the scatter in n the data, data comparison _ 


and the selection of design values are presented below. 


a -Near- Surface “Material. —In comparing the laboratory and field values, the 
thickness of the near-surface layer was determined to be 14 ft (4.3 m) in the 
laboratory. This was based awed on the boring logs; however, in the field, 


m) ‘thick. The modulus of 3.5 x 10° p Psi (24 MPa) 
compares very favorably with the field-determined value of 3.3 x 10° psi (22.7 


computed from the field date was determined to be 9.7 7 x 10° psi (67 MPa). 
This value seems reasonable; however, it is about one third of the laboratory-de- 
termined value, which is 26 x 10° psi (180 MPa). This value is judged to be 
too high. For the top layer, v values of shear and Young’s moduli selected = 
design were 3.4 x 10° psi and 10 x 10° psi (23.4 and 69 MPa), respectively. aa 
Second Layer.—Immediately below the near-surface layer is a second clay 
layer that ends at 20-22 ft (6-6.7 m) depending on whether the field ae 2 
_ data or the laboratory classifications of samples obtained from the boring are 
_ accepted. The shear modulus as determined in the laboratory for this layer ” 
is 2.8 x 10° psi (19.3 MPa). This value is slightly less than half of the shear 
modulus 6.6 x 10° psi (45 MPa) that wes determined in the field. Because 
the soil is predominantly clay to a depth of 20 ft (6 m), and because the shear 
o- modulus should increase with increasing Gepth, a design shear modulus of about — 
oe 10° psi (34 MPa) was selected. The values of Young’s modulus range > 
from 19.7-24.4 x 10° psi (136-168 MPa) depending on whether one chooses 
"the field or laboratory values. These values are relatively close (within about 
20% of each other), and this is considered good agreement. A value of 20, 
x 10° psi (138 MPa) has been s selected for design purposes. 
_ Third Layer.—The third layer selected for design purposes combines the tw 
i layers within the depth interval of 22-31 ft (6.7-9.4 m) in the laboratory- aes 
- values and is assumed to be at a depth of 20-30 ft (6-9.1 m). If one averages” 
> the shear modulus values of 16.7 x 10° and 5.56 x 10° psi (115 and 38 MPa), 5 


one obtains 11. 1 x 10° psi (76.5 MPa). This value is in close agreement with 
12.6 x 10° psi AL MPa) that was obtained in the field, and consequently a 
value of 12 x 10° psi (83 MPa) was assigned for design purposes. In the laboratory, 
averaging the Young’s moduli for the two layers between 22 and 31 ft (6. 7 
and 9.4 m) yields a value of about 50 x 10° psi (344 MPa). This is considerably 
higher than Young’s modulus of of 37.1 x 10° psi (256 MPa) as determined by = 
geophysical methods. Based on engineering judgment, a value of 40 x 10° 
psi (276 MPa) has been selected for design purposes. 
j Fourth Layer.—From a depth of 30-36 ft (9.1-11 m), the soil consisted of 


sandy silts with gravels and consequently was not suitable for laboratory testing. 

Shear and Young’s moduli of 86 x 10° and 225 x 10° psi (593 and 1,550 MPa), 7 au 

respectively, were assumed for design purposes. These values are based on 

* _ Foundation Rock.— —Finally, the shear and Young’ s moduli of the foundation 


- rock were determined to be 125 x 10° and 370 x 10° psi (860 and 2,550 MPa), ; 


4 
= 
4 
| 
| 
Tespectively. ese values are Dased solely On the field geophysical data. 


The data presented in Fig. 22 show that the shear modulus Aseiecitaall in 
the laboratory is approximately equal to or less than the shear modulus obtained : 
_ by in situ geophysical techniques. This trend has been reported by others (10) ol 
: can be attributed to the following: (1) The data obtained in the laboratory | 
are influenced by such things as sample disturbance, equipment limitations, 
and boundary conditions; (2) the geophysical techniques produce data at a much 
smaller strain than the laboratory tests; and (3) time effects that tend to increase _ 
the laboratory Gutormined shear have not been included. 


= soil anisotropy | is “responsible for ‘the: laboratory values being greater il 
~ field values. In the resonant column technique, Young’s modulus is measured 
_ perpendicular to the bedding planes. The geophysical tests measure Young’ s 
modulus both parallel and perpendicular to the sample bedding planes. The 
crosshole tests and the surface refraction seismic tests provide P-wave velocities _ 
_ for travel paths parallel to the bedding planes. The downhole tests provide 
“= P-wave velocities perpendicular to the bedding planes. Comparison of the P-wave — 
velocity data obtained by the three geophysical techniques does not substantiate 
- the hypothesis of the soil being anisotropic (Fig. 8). On the other hand, these 
= materials are glacial lake-bed deposits, and the geology and mode of deposition : 
support the hypothesis that the soil is anisotropic. This cross-anisotropy is believed 
to be responsible for the data obtained in the laboratory being higher -_ 
- An additional explanation for the laboratory determined Young’ s moduli being 
higher ow the field values is the possibility that the laboratory specimens 
did not adequately represent the entire in situ soil layer. The qroghyuicel 7 
techniques evaluate several cubic meters of soil while the laboratory technique 
is limited to several cubic centimeters of soil. As a result, the specimens tested : 
__ in the laboratory may or may not adequately represent in situ conditions. = — 
A third possible explanation for the laboratory determined Young’s — 
being greater than the field values could be the difference in percent saturation. 
previously the third layer from the ground surface which has 
s a P-wave velocity of 2,800 fps (850 m/s) is not believed to be | > fully saturated, 
; and the resonant column test specimens were not back- -pressure saturated. 
_ Consequently, it is difficult to say whether the laboratory specimens were more 
or less saturated than the in situ soil. Nevertheless as the percent saturation 
Se 100%, a change in percent saturation of 0.1 or 0.2% can cause : 
a a change in -% -wave > velocity of several hundred feet per second. One cannot 
a ~ isolate the determined cause of the high laboratory Young’s moduli, but it could — 
be one of the hypotheses mentioned above or a combination of these ypamee. > 
be used for low strain calculations. This value is based on resonant column — 
data and is valid only to a depth of 31 ft (9.4 m) because the soil below that 


depth was not tested. However, because the damping appears to be uniform 


| 
| | 
i 
a 


with to the bottom layer [30-36 ft (9.1- ‘m) in 


Lait 
Based on the data the at 821 UD 
can be subdivided and zoned for the purpose of dynamic analysis of building _ 
or machine foundations Fig. 22. This figure shows four soil layers above rock. : : 
The near-surface layer is 8 ft (2.4 m) thick, the second layer is 12 ft (3.6 _ 
2) thick, the third layer is 10 ft (3 m) thick, and there is a 6-ft (1.8-m) thick — 
layer immediately overlying the rock. This is an idealized soil profile that would | 
suitable for design calculations. The values of shear modulus, Young’s 
and damping ratio for each of these layers to be used in design are listed — 
in the right hand three columns of Fig. 22,00 
_ When there exists considerable scatter in the laboratory and field data obtained 
: ‘for a specific soil layer, WES usually places more emphasis on the field data. 


en The field investigation usually involves the use of several procedures (such 
as the crosshole, downhole, surface refraction, and vibratory techniques), which — 
7 2. The laboratory investigation is usually limited to resonant column testing. 5 
3. The data obtained in the laboratory are influenced by such things as s sample ‘ 
disturbance, equipment limitations, and boundary c conditions. 
a Furthermore, the laboratory data are only as valid as the specimen tested 


by a volume of soil ‘mensesed in terms of cubic meters not cubic centimeters | 


3 representative of the soil layer. In this regard the field investigation is iefinenced— 
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wrieg By Arshud Mahmood,’ Clarence J. Ehlers,’ and Blase A. Cilweck, 


are used in in on land to 
‘supplement borehole data and other site investigation information. In studies 

_of soil conditions offshore, high-resolution geophysical tools and methods are 
“used on a nearly routine basis in conjunction with bottom sampling and soil ; 
_ borings. Although the water at an offshore site inhibits the familiar and reassuring — . 
experience of a visual reconnaissance. Since | water isa medium i in which 


a a rapid and comprehensive geophysical investigation. 
_ In 1973, acoustic (geophysical) methods were used to investigate » the e seafloor _ 4 
geological conditions, particularly sand waves, in an area of potential offshore _ 

development, lower Cook Inlet, Alaska. The bedform aspects of the survey 


are reviewed here. The term “sand waves” is generally used for 
_ wave- shaped bedforms composed of granular m material and formed at the boundary > 


between flowing water and a mobile bed of unconsolidated sediment. Sand — 
waves are a potential foundation hazard, particularly to submarine sein be 
due to the possibility of their | 
Regional Setting. .—Cook Inlet is a 190 mile (300 km) long basin, a Sneonatint j 
expression of a structural trough bounded d by mountain ranges. The portion | 
of Cook Inlet south of Kalgin Inland, generally known as lower Cook Inlet, 
is underlain by a sequence of thick sedimentary rocks that extend across much 


Of southern Alaska. Several quaternary glacial episodes prior to about 15,000 
: years ago have affected lower Cook Inlet and resulted in the deposition of | 
5 considerable thicknesses of granular sediments, which ae nie ef have been 

Geophysical Survey.—The acoustic data were in 1973 consist 
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"of 3,400 line-miles (5,440 line-kilometers) of multisystem, high-resolution, marine P 
7 j acoustic survey data collected on a variable grid shown in Fig. 1. The survey __ 
a _ was performed using a 4 kilojoule sparker, as acoustipulse, a tuned transducer, Te 
fathometer and dual side-scan sonar. In each system, a sound source in water — 
produces sound energy, and this energy is transmitted down to the — ; 


| When this sound energy strikes the seafloor (or any other object having acoustic 


_ properties different from water) a part of it is reflected back to a sound receiver 
7 q an echo. This sound energy then is reconverted to electric energy, amplified, __ 
- and recorded on a graphic recorder. The amount of time that lapses between ¢ 


“ 
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FIG. 1.—Multisensor Acoustic Survey (1973) Lines in Lower Cook Inlet, Alaska 


_ Generally, several acoustic devices can be deployed simultaneously from a — 
sufficiently large boat. A typical equipment layout is shown in Fig. 2. The i 
_ three systems most useful in interpreting sand waves were acoustipulse, fathome- 
ter and side-scan sonar. Specific features of these three systems are summarized 
in the section entitled a pce Systems. Characteristics of all five Systems 
a All of the acoustic systems used in this study operate on the principle whereby 
transmitted sound energy, incident upon an acoustic interface, is reflected from — 


q _ the pulse transmission and the echo reception is a measure of the distance 
te ssion eceptic as 

t 

A 


| / § : 
| 


the: interface An acoustic interface can be defined as any interface across 
which there is a contrast in acoustic properties. Such contrast is dependent 
on the acoustic impedance of the materials on each side of the interface. The © 
characteristic acoustic impedence of a soil layer is ] = p,V, in which p, = 


the soil bulk density and V = the velocity of a compressional wave, such 


as a sound wave traveling in soil layer. adi Ww 


Navigation Equipment 
and Antenna 


Side Winch 


‘a 
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The amount of energy that is transmitted to the next interface depends on | 
the contrast in acoustic impedance between the two ‘materials. ~The contrast 
is symbolized by the Rayleigh reflection coefficient, R, which at normal comet id 
is expressed as R = J, — I,/1, + 1,, in which J, and J, = impedance values 
for two layers on either side of an interface. Thus, boundary between 
very dense lange acoustic and an am overiying soft clay 


| | 
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‘with a small impedance will act as a strong reflector. ‘The energy 
within each layer is frequency-dependent, 
_ expressed as a = kf” in which a 
f= frequency in kHz and n 


"increases. For this r reason, “the depth scale on the profiles « cannot be ‘considered 


attenuation in db/m, k = a constant, 
= the exponent of frequency 


a 


and the relationship is “generally 
4 


truly linear, and hence additional information is required for accurate calibration. — 


TABLE 1 —Comparative Description of Marine High- Resolution Geophysical Systems =a 


pon. 
System 
— 
provide map view 
of seafloor to 
(a) locate 
wrecks, pipelin- 


es/debris, (b) 
determine sea- 


Frequency 


of output 
pulse 


Operational 


resolution 


| 


variable, typically 
lto several 
feet (O3toa 


<J 


105 + 10 kHz 


‘Type of 
acoustic source 


Power output 


128 db (Peak) 
Ref. | microbar 


piezoelectric 
transducer 


few meters) 


floor topogra- 


phy. 
detect gas bub- a 

bles in water. 
determine water 


depth 
detect gas bub- 


bles in water 
column 

provide profile of 

Tuned Trans- | provide profiles 

| of very shallow 
subseafloor to 
about 100 ft (30 

determine water 
depth 

detect gas bub- 
bles in water 

* | provide profiles 

of shallow sub- 

seafloor to 500 


Fathometer 


determine water 
depth 


bles in water. 
provide subsea- 
floor profiles to 
intermediate 
depths about 
3,000 ft (900 m) 


=| 


detect gas bub- 


2 ft (0.6 m) 


variable, typically 
(0.6to 


iss 
variable; typically 
2h 
depending on 
application 


about 25 ft (7.6 


& to 4 kHz 


Alternatively, , a close epgvidianation on of shallow _— depths can be made 


piezoelectric 
‘transducer 


up to 2 joules 


provides excellent 
resolution of 
shallow faults, 7 
bedding. gas 
zones, etc. 


200-1 ,500 


electromechanical 
twaneducer 


most 


penetration but 
least resolution 
of the subsea- 
floor profiling 
systems used 
for site sur- 
veys. 


by assuming a velocity previously determined for other similar materials. ac 
All of the acoustic systems except the side-scan sonar system are designed “S 


to provide data essentially from beneath each sound source; no data is received 
from either side of the line of profile, except for occasional spurious 


‘ide 


pemarenere In contrast, the side-scan sonar system provides data from several 


3 
§ = 
alsoreferredtoas 
40-200 kHz about 0.5 joules to 
ric by also referred to as 
broadband 2¢ 
¢ broadband 10 electric Spark | 3.3, 4, 0.0 oF 
| 
ch *Tradename. 
4 
_ geophysical systems from which data are used in studying sandwaves, Jf 
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“AG. 3.—Acoustipulse Data Showing Sand Waves (A and B Refer to Corresponding x 
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FIG. §.—Acoustipulse Data Data honing, Sand Waves and Old Ere Erosion on Surface. (Dashed 
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acoustipulse, "side-scan sonar, fathometer, are described briefly in a later 
Interpretation of Geophysical | Data. —Interpretation « of f geophysical data from 
> lower Cook Inlet survey showed that sand waves were present in the northern — 


portion of the study area. The height of sand waves generally observed ranged © 
from less than 2 ft (0.6 m) to over 50 ft (15 m), and the wave length ranged 
= less than 20 ft (6 m) to over 1,000 ft (305 m). The height to wave-length 
ratio was generally between 1:15 and 1:50. Several shapes of ripples and dunes — 
were observed. Examples of acoustic data showing sand v waves are given | in 


Es 


AG. 6. —Sand Wave in Lower Cook inlet, Alaska 
oo smallest detectable sand wave was about 2 ft (0.6 m). Sand waves of 
5 ft (1.5 m) in height or greater were considered to be significant, since local _ 
seafloor elevation changes exceeding 5 ft (1.5 m) over a short distance could _ 
affect the installation and performance of structures. Sand 
with a height of 5 ft (1.5 m) or greater were ‘mapped, , and their distribution 
is shown in Fig. 6. The area of most numerous sand waves contained waves - 
in over 80% of the geophysical profile lines, and is about 1,200 sq mile (3, 100 
s Sq km) in size. The area of moderate numbers of sand waves had about 50% © 


 Seafloor materials were sampled using a grab 


= were with significant | sand waves. 


— | 
| | 
| 
igher than m). In = 
n 20% of the survey line lengths _ 


weighing 2 about 100 Ib (45 that achieved a of about 6-8 in. 


ae in five sampling areas selected on the beni of features seen on the geophysical _ : 
data. Some of the 65 sample sites of area No. 3 were in | the area of sand 


& The survey vessel was located over each sample site by means of an electronic 
positioning system. The grab sampler was lowered to the seafloor and ayrrner-l 
_ by means of a winch onboard the survey vessel. Each sample was el 
observed, logged, and then a representative portion placed in a water- a a. 


“were mixed with an equal ‘amount of water, stirred thoroughly, and then 
centrifuged. The water was decanted and, using a 10 power binocular microscope, — 
the apparent clay/silt and silt/sand interfaces were visually determined. The 

clay, silt, and sand volumes were expressed as percentages of the total sample 

a All of the 64 samples analyzed from sampling area No. 3 were predominantly 
r ‘sand. About 50 of the 64 samples contained more than 90% sand by volume. _ 
Visual examination placed the sand mostly in the fine to medium grain size 
range with a 50 percentile size, D,,, of about 0.04 in. (1.0 Sos 

- Sand Waves—Terminology.—A variety of terms and definitions have been : 
used in the literature to describe the various shapes and origins of sand waves. a 
_ The two disciplines that study sand waves are sedimentology and hydraulics. 

_ There is no general agreement on terms and definitions between the two areas _ 
of study, nor even within each discipline. Many comprehensive publications — a 
present an ad hoc set of definitions to establish a basis for examination. One - 
€ the main causes of confusion is indiscriminate mixing of terms based on 
‘origin, shape, or size. No attempt is made here to resolve all the differences, 
as some of them result from geologic settings or other specific ‘study 


‘* In this paper, the term “sand waves” is applied in a generic sense to all 
wave-shaped bedforms composed of granular materials. The description of I 
seafloor features in lower Cook Inlet is mainly confined to their size. A variety 


shapes and sizes were observed, some of which fit the Shapes and sizes 


implied by ripples, dunes and ridges, without implying any y specific origin, 
_ In the theoretical analysis on sediment transport flow regime, a clear classifica- 
tion can be made on the basis of origin. Dunes and ripples are asymmetric — a 
sand waves that are formed under tranquil flow conditions (Froude number, 
F, is less than 1). The symmetrical wave forms termed ‘‘antidunes’’ are formed 
under shooting flow conditions (F greater than 1). tb wo q 
= Influence of Sand Wave Movement on Offshore Facilities—Because of the 
- possibility of their movement, submarine sand waves constitute a potentially 
hazardous foundation condition and place constraints on offshore exploratory 
drilling and construction. The movement of bedforms is of primary concern 
. the construction and performance of offshore pipelines. As extreme sand — 
waves with side slopes representing over 50 ft (15 m) height change in a shorizontal 
_ distance of 200-300 ft (60-90 m) were observed in the study area, burial of 
Seen below the potential depth of scour would be a teem Movement 


| | 
| 
q 
an 
| 
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areas of numerous waves, of pipelines, use 
q heavier pipe and periodic monitoring of the seafloor may t be required. Movement _ 
«Of sand waves in the upper Cook Inlet has caused several pipeline failures, 
_ and vigilance is maintained by means of yearly side-scan sonar surveys to identify 


areas of bedding support loss(16), 


Fer pile-supported structures, partial loss of lateral support would result due idl 


to sand movement away from the structure, analogous to deep scour. Such — 


: loss of support would have to be anticipated in lateral load design Of owl 


structure piles. For gravity structures, local relief associated with sand waves | 
would influence site selection. Such relief would also influence the siting of | 
mat-supported mobile jack-up rigs. For leg-supported mobile rigs used in oil 
exploration, leg penetrations are expected to be shallow if the seafloor soils 
consist of sands. Thus, scour due to possible sand wave movement could 
undermine the leg footings and present a serious threat to rig stability. __ ont? 
‘The presence of sand waves in an area generally implies that some sediment _ 
" movement may be occurring. However, the rate and frequency of movement — 
on be inferred by observing wave forms. All sand waves are in a dynamic 
condition during formation, and the movement continues as long as similar 
——- prevail. The absence or presence of sand wave movement under — 
_ the prevailing conditions in lower Cook Inlet is difficult to assess. An indication _ 
of movement has been sought through an approach using the following steps: Tv 


1. Examination of the two phase (water-sand) flow regime from a knowledge _ 
of grain size and current data. 


2. Direct observations (diver or TV) to ascertain bed movement. __ ae 

geophysical to examine the bedforms for changes. 
Flow Regime. —The 64 bottom samples described earlier showed that the shallow _ 
: seafloor materials consist predominantly of fine to medium sand, with an = 
_ approximate mean grain size of 0.04 in. (1 mm). Data on bottom ‘corent 

_ observations is quite limited, and shows a mean velocity at maximum flood 
tide of 1.0 knot (0.5 m/sec) at one location, although some measurements bie. 
exceeding 1.4 knots (0. 7 m/sec) were also recorded (16). These velocity 


q 


4. Based on relationships between flow velocity, bed shear, and grain n size, a 7 
l 


0-knot (0.5 m/sec) current imparts a bottom shear stress of 0.25 lb/ft? (dl. 2 


N/m? or 12 dynes/cm’) as computed from Hjulstrom’s and Shield’s diagrams in 


(1). This shear stress value can be used to assess whether the flow regime — 
: is conducive to the formation of sand waves, which would imply movement. 


- By computing grain size, Reynolds number, and the mobility number for 0.04 a 


in. (1.0 mm) size sand sampled in lower Cook Inlet, and assuming a 1.0 knot 


(0.5 m/sec) current, the lower Cook Inlet flow regime can be compared with : 


data obtained under controlled conditions. A plot of Reynolds number and mobility 
in Fig. 7 shows that the Cook Inlet data ‘fall in a zone 


] 
I 
The sand wave formation is also a fi evelopme . 


duration can be as long as several days at small flow velocities (17). Flow 

a velocities of 1.0 knot (0.5 m/sec) acting for short periods and equally in both 
F| directions of tidal flow, are probably insufficient to form significant sand waves. 

Under these conditions, bedform development might take place only during 

severe storms when current velocities are high. The complex seafloor and shoreline — 

topography could also affect current velocities during storms and periods - 
_ wide tidal ranges. The applicability of the laboratory- -based relationships —— 
a in Fig. 7 to field conditions such as those found in lower Cook Inlet is uncertain. > 
| However, sediment transport concepts have been applied to the examination 
4 of phenomena ranging from desert dunes to sand waves (17). Detailed computa- 
tions of flow regime appears in a later 


— 


” Direct Observations.—TV camera observations of the seafloor at the end of 
a flood tide by a U.S. Geological Survey (USGS) study team (2,3), showed 
_ bed load transport occurring in the form of rolling and saltating grains. This 
*s to the conclusion that active transport takes place in lower Cook Inlet. | 
As only a limited area could be observed for a brief time period, movement _ 
of bedforms s larger than ti the smallest — could not be ascertained by direct 
e Repeat Geophysical Surveys. —During +s 1973 survey reported in this paper, 
duplicate fathometer coverage of portions of the survey grid was obtained while 
performing separate sparker and acoustipulse a fathometer records 


= sedimen ovemen 
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obtained two ture surveys 3-52 52 days apart. A compatison 
of fathometer profiles provides information concerning the short-term movement — 
Comparison of these ‘repeat. fathometer records ‘shows that ‘the ‘seafloor 
configuration is generally stable, but locally there are differences in the seafloor | va 
profile. Each limited area of possible movement is bounded by other areas | 
that did | not appear to move yee the survey span. An example of poets 


wave movement were re found. in "which postplots indicated that the later and 


earlier survey lines coincided. Thus, navigation errors can be ruled out. io 


USGS teams who resurveyed six lines from the 1973 survey reported here a 
| (2,3,4,16). They concluded that the sand-wave fields in lower Cook Inlet appear 


0 be stable with very limited changes in bed-forms over five years. aoe 


= 


a = _ An assessment of long-term movement was attempted in 1977 and 1978 by 1 


FIG. 8.—Examples of Sand Wave Movement and 
due to navigational problems, data 


This data from a Suite of marine high- resolution 
refers 
g 


= define the water depth and provide a cross-sectional display of the sediment 
_and rock layers below the seafloor. A more appropriate term is high- — 
acoustics or continuous acoustic profiling. The three devices used in the — 
wave study reported here are described in this section. The reader is referred — 


to other publications for a more complete description of the systems (12,13), 2 
for details of data processing techniques (14), for application of methods © 
@, 10,11,12, 15), and for marine socunent properties that affect the Propagation 


| 
— 
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and attenuation of sound energy (5,7). The following denctigtions are adapted 


Fathometer or Depth- -Sounder.— The depth-sounder system consists of a power 
- ‘supply, a transducer (a transceiver that alternately transmits and receives sound), _ 
_ and a graphic recorder. The depth-sounder transducer usually is installed in a 
the hull of the vessel, midship, and 6-10 ft (2-3 } m) below the water surface. — 
_ ‘The transducer converts electric energy to sound energy, and the sound energy 


is transmitted downward to the seafloor. When this energy strikes the seafloor, . 


« any other object having acoustic properties different from water, a part 

of it is reflected back to the transducer (transceiver) as an echo. . This sound 

= then is reconverted to electric energy, then 1 recorded on a graphic r recorder. ; 
Assuming that the speed of sound in water is nearly constant, the amount am 
3 time which lapses between the pulse transmission and the echo reception — 


from an earlier paper that describes the geophysical systems used in this — F 
(l 


isameasure of the distance traveled. 


Because of the high frequency of the fathometer, virtually no subbottom — 
penetration is obtained. This system is is suitable for operating in water r depths — 


_ The accuracy of a depth-sounder system is largely dependent on two variables: 
(1) The velocity of sound in water, which can vary with changes in temperature 
and salinity; and (2) the precision of the electronic timing circuits used. Both 
these variables necessitate calibration the system for accurate 


Side-Scan Sonar.—Side- -scan sonar systems provide graphic records that show _ 
a two-dimensional (map) view of seafloor topography. These records are closely 


: _ analogous to a low-oblique aerial photo. Features shown on the side-scan sonar _ 


commonly appear similar to their natural perspective. Sr) 


_ A typical side- “scan system consists of a towfish containing two arrays” of 


- transducers and a shipboard graphic recorder. The transducers emit 105-kHz _ 
_ pulses as fan-shaped beams on each side of the towfish. The beams are 

perpendicular to the direction of vessel travel and are broad enough in the 
vertical plan to extend from beneath the towfish to the full system range of — 

: 1,500 ft (450 m) to each side of the towfish. The beam width is only 1. 2° = 

} in plan view. Each pulse lasts 0.1 msec and the repetition rate varies | depending . 

4 Signals reflected from the bottom and from objects on the bottom are e displayed _ 
on a continuous graphic record produced on a two-channel, wet paper recorder. 

Rs The intensity and distribution of reflections received depend on the composition 

and surface texture of the reflecting object, its size, and its orientation with 
respect to the transducers in the towfish. Due to the beam shape and the : short 

; length of the transmitted pulse, the side-scan sonar can resolve small objects 
on the seafloor and details of topographic irregularities. = 
_ Acoustipulse.—The acoustipulse system is an electromechanical acoustic pro-— 
filing system consisting of power supplies, trigger banks, shipboard electronics — 
- package, filter graphic recorder, and a towed transducer array with an ac- 
companying towed hydrophone. Acoustipulse is a single channel, high resolution, z 


S a sil 


marine reflection seismic system developed to provide both high resolution and 


4 up to 300 ft (90 m) in by and soft-rock strata. 
= 


| 
| 
— 
f 
a 


= 
aa seafloor penetration is aon 250-300 ft (70-90 m). Penetration of as Se 
‘a much as 500 ft (150 m) can be obtained under the most ideal acoustic conditions * 
‘(where the material being profiled is soft). “4 
3 The acoustipulse waveform is produced by discharging stored electrical energy 
zz in a potted wire coil. Eddy currents set up in the coil repel an adjacent — 
aluminum plate. A partial vacuum created by the flexed plate returns the plate 
q its original position. The return motion is highly dampened | by a rubber Os 
diaphragm, thereby minimizing source reverberations. The source is driven by 
a capacitive-discharge power supply capable of producing up to 1,000 mend 
_ Normally, a source consists of three electromechanical transducers, all of a 
which are fired simultaneously every 0.5 sec. Energy output of each transducer 
is variable, all of which are mounted on a 8 ft (2.5 m) long catamaran and = 
- towed at or near the water surface 100 ft (30 m) behind the ‘survey vessel. 
Reflected signals are received by a 10 element linear hydrophone array — 
at or near the water surface and abeam of the source transducers. . Signals 
outside the desired frequency range are removed by frequency filtering. —_ 
_ data are displayed graphically on a 0.25 sec record. Vertical exaggeration of s 
the subbottom profiles is typically 10-20 times. Signals also can be recorded _ 
on an analog tape recorder for laboratory playback. edt 
Experiments indicate that a mobile granular sediment bed acquires sand wave > 
_ shapes for certain intervals of the dimensionless variables of two-phase motion 
i (17). If the set of variables expressed as Froude number recuces to a value 
less than unity, then the flow is considered tranquil, and sand waves formed _ 
_ take the asymmetrical shapes of ripples and dunes. These may be stationary 
or move along the current path. If the Froude number is greater than a 
_ the flow is shooting, and symmetrical sand waves termed antidunes are formed. 
_ These may be stationary or move upstream against the the current path. The Froude. 


in which u = average velocity of Me flow; gus = acceleration due to gravity; he 


_ In the lower Cook Inlet study area, sand waves are found in water — 
ranging from 100 ft (30 m) to 426 ft (130 m), and flow velocities of one knot 

— (about 0.5 m/sec) have been measured. For these conditions, u = 1. 0 knot ; 
(0.5 m/sec), h = 100 ft G0 m), and g = 32 ft/sec? (9. 81 m/sec” ). F, = 

0.5/V(9.81 x 30) = 0.03, which is |. Therefore, if the full water 

; depth is involved in tidal flow, then the flow will be tranquil, possibly conducive 

_ to the formation of asymmetrical ripples and dunes, but not the symmetrical : 

_ antidunes. As some symmetrical forms are present in the study area, either 

_ the local flow velocities are over 20 knots (10m/sec), which would be improbable, 
or the flow depth involved is that of a shallow current about 6.5-10 ft (2-3 


‘m) deep, and flow velocities exceed 5 knots (2.6 m/sec), probably during severe 


assess whether the present flow regime is of moving, or forming, 


 S& bedforms observed in acoustic profiling data, the field measurement data 
can be plotted on the X-Y plot in Fig. 7 (17), in which X= grain size aware 
~ number and Y = mobility number, as defined below: Grain size Reynolds number, 
ae = vD/v in which D = ‘om sine;,s v = kinematic viscosity; v= shear velocity 5 


=v (t,/p). In which +, = shear sti stress. between the flow and bed 


surface’ = density of fluid = a in which y = specific unit weight of 


Mobility number Y= 


which 7, = specific. unit weight of grains, and v and D are as 


(quartz) = | 165 (2.6 64 g/cm’) 4 


Current t velocity = 1.0 knot (0.: 5 m/sec which, and 
Shield’s s diagrams (1), imparts a bed shear s stress 1 = 0. 025 lb/ft? (1.2 N/m* 
12 dynes/cm’). Therefore, shear velocity v (*,/p) = VO. 025/2) = 


0. cm/sec). water at 39. 2° 4 which is close to 


——=216.... 
=2x——x33x 107 = 0.044 


: a _ Considering Fig. 7 with the above values of X and Y, the conditions are iaaibble 


“for simultaneous formation of ripples and Gunes. 

are an sor part “of offs offshore site 
Selected high-resolution geophysical systems can be used to study the seafloor — 

and shallow subseafloor zones to identify the bedforms. A portion of the seafloor | : 
in lower Cook Inlet is covered with waveshaped bedforms composed of — 
‘een ranging in height from less than 2 ft (0.6 m) to over 50 ft (15 m), a 
with | wave lengths from less than 20 ft (6 ‘m) to over 1,000 ft (305 ‘m), an 
with height to wave- length ratios generally between 1:15 ise 1:50. 
_ Wave-shaped bedforms are generally associated with sediment movement in 
‘strong current areas. Since sand wave movement could put constraints on offshore 

construction, a multiplicity of techniques’ was applied to assess the movement | 

i of these bedforms in lower | Cook Inlet. An examination of the flow regime — 

d shows that the prevalent currents in lower Cook Inlet are capable of ‘initiating nit 


some movement, but apparently the tidal a are too short to form the 7 
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storms when current velocities are high. Reported direct TV observations support 
= the view that there is some facalinn® sediment movement. Repeat geophysical 


"areas with some movement in 1 localized areas. Further | investigations, s, including 4 
~ post- -storm inspection of actual construction sites, is required for more definite . 
In siting offshore facilities in an area of sand waves, a map similar to Fig. 
“s presenting sand wave distribution on will be very helpful in the preliminary selection — ¥ 
of the most suitable or feasible site(s) < or route(s) for which | the incidence of ‘a 
sand waves is lowest. Final selection and evaluation of a structure site or pipeline 
route will - Tequire consideration of the potential for sand wave movement and _ 
the depth of scour associated with the movement. The geophysical surveys 
~ conducted in lower Cook Inlet have demonstrated that high-resolution 1 geophysical — 
can be used to to identify bedforms, determine the ‘size, shape, 


by ; former BBN- Geomarine Services Co., now a part of McClelland Engineers. 
An ennames cleat permitted publication of Fig. 6, which is based on exclusive | 
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of each for laboretcry use. 
_ The following symbols are used in this paper: 


grain diameter; 
= Froude number = w u/V (gh); 
acceleration due to gravity; 


characteristic acoustic impedence; of well 
a constant relating acoustic attenuation to frequency; : 

sand of the sand 
exponent of acoustic 


i: 
= V(t./e); 


wave 


Shields function. The value = Y.., signifies sediment motion; 


specific unit weight of fluid; inte the vibrator 
specific unit weight of grains; ibe cond, at 


density of fluid = 
To = shear stress interacting between flow and bed surface. 


q 
| 
ware h — 
" 
ig 
We 
L = grain size Keynoids number = v v; 


' 


Mom 


ate 


| | 
BENG: 
es | 
—- 

| | 


od} 


SEISMIC TECHNIQUES IN THE LABORATORY“ 


Several seismic exploration techniques are used routinely in 
engineering to characterize earth materials at field sites. Some of these seismic _ 
_ techniques can be used to advantage at a laboratory scale. While the principles 
_ remain the same, some of the details of performing scaled down tests are different. _ 
- Four applications of seismic techniques in the laboratory will be described — 
an emphasis on differences between the normal field application of each techniqu 


and the adeptation of each for laboratory use. to tet 
< A facility was built by Woods (13) to perform steady-state elastic wave screening 
_ studies. This facility, in Fig Wi was about 7 m in diameter and 2 m deep and 
was s constructed with an absorbing perimeter of wet sawdust tor minimize reflection. 
of wave energy. The central testing zone consisted of well graded sand (common | 
mortar sand) at a void ratio of about 0.5. The moisture content of the sand = 
was maintained at between 2.5% and 6% to provide apparent cohesion and 
_ to permit excavation of open trenches and holes for wave screening studies. 
— Itw was necessary to know the shear wave v velocity y profile in the half-space 
facility to analyze adequately dynamic model footing and screening studies. 
Woods (13) used the steady-state Rayleigh wave method described in Richart, 
Hall, and Woods (8) to determine shear wave velocity as a function of of depth. 


.. solid circles in Fig. 2 represent the shear wave velocity profile. ete: we 
_ The electromagnetic vibrator, shown i in Fig. 1, was used to generate Rayleigh 


waves, and the distance to in- phase points was determined by comparing, on a 
a dual beam oscilloscope, the signal from a transducer built into the vibrator xa 
with the signal from a geophone moved along the sand surface at increasing ~ 


2 “Presented at the October 22-26, 1979, ASCE Annual Convention and Exposition, held , 
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- ASCE. Manuscript was submitted for review for possible publication on February 17, 7 
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radial distances from the exciter. Varying the frequency of the exciter provided » 
energy with different wavelengths which as described 
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FIG. 1.—Model Half-Space Facility (13) 
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2.—Shear Wave Velocity Profile in Half-Space 0) 


Ren major difference between these s' steady-state tests and those normally 
_ performed in the field (2,3,5) was the size of the exciter and geophone. In 

both instances smalier units were used in the half-space facility. The exciter 


was an M-B Electronics, Model C C-31 with a maximum force seme of 112 


J 
[AI 
is 
300 
q 


SEISMIC TECHNIQUES 


N while field units at least 1.4 kN capacity. 

A an Electro-Tech, 8.5-Hz unit which was chosen because its base size, = 
diameter, was small compared to the wavelength of the shortest Rayleigh wave A 

z about 225 mm. In addition, because the wavelengths were short, more care | 


was used in matching phase and measuring distances. 


_ Stokoe (10) also performed small scale cross-hole seismic tests to determine 


the shear wave velocity profile in the same half-space The cross- -hole 


& Results of Stokoe’ s ; tests are also plotted « on 1 Fig. 2. The cross- “hole velocities 
I are consistent with the earlier steady-state velocities and with empirical equations 
os proposed by Hardin and Richart (6) at depths greater than about 0.3 m. — 


AP 


"total: unit Ln gl times depth; and K, was assumed to be 1/2. 

, the shear “wave v velocities by cross-hole were al) 
higher than expected. This 1 was ‘nelle due to waves fotiawing a curved vo 
= deeper, more confined soil as described by Haupt(7), 

_ The major differences between usual field cross-hole tests and ieee performed — 
- the half-space facility were the spacing between bore holes and the amount _ 
x of energy (size of hammer used) introduced into the soil. ‘Spacing shorter than 

usual of boreholes necessitated the use of faster sweep rates on the oscilloscope 
to determine travel time and careful measurement of the distance between source ~ 
_ and receiver bore holes. The center-to-center distance between bore holes was _ 
_ checked with plumb lines at each depth at which data was obtained. Plumb 
+ lines were possible because the boreholes were only a maximum of 1.22 1 m 
deep and they could be kept quite straight for that depth. 
This half-space facility having been ‘‘calibrated”’ by the described 
_ tests was used later for initial = of the in-hole Po for field work described — 
“by Auld (1). ~ 
Move. Facury 
Woods, Barnett, and Sagesser_ (14) used a 
~ containing moist, crushed quartz, fine sand at a void ratio of about 0.7 in 
_ model tests to study the effectiveness of cylindrical holes as barriers to screen 
surface waves. A moisture content in the sand from about 2.4%-6.3% ‘provided = 
adequate apparent cohesion to permit the excavation of open barriers. 
Circular footings 30 mm in diameter were » placed on the surface of the sand 
or embedded in the sand and excited in steady-state vertical vibration by an 
- electromagnetic exciter (Fig. 3). Stopped motion interferograms of caveling 
: Rayleigh waves were obtained by a stroboscopic, holographic interferometry _ 
_ technique (Fig. 4). In these interferograms the dark and bright fringes can be — 
- interpreted as contours of equal vertical displacement. Closely spaced contours © 


represent steep inclines, while spread contours are crests or troughs of the 


| 
2, was computed from 
| 


To determine the dynamic shear modulus of the sand 
were made for various frequencies of excitation and the wavelength of the — 
Rayleigh wave was scaled directly from the interferogram. From the measured 7 
wavelengths and known excitation frequency, the Rayleigh wave velocity pwr ; 
wavelengths ranged 


| 4 be computed. Using frequencies from 575 Hz-1700 Hz, 
from 75 mm-165 mm. Using the half-wavelength effective depth as described 


im Richart, et al., (8), the Rayleigh wave velocity versus depth was plotted 

(Fig. 5). Note that the depth sampled in these tests was very shallow, less 
than 85 mm, but, more importantly, the velocity een rather — 


al 
aM 


bad 


than the as be expected. A explanation for this behavior 

, a based on the relationships developed between water content and apparent — 
= in the fine sand. Table | shows the results of direct shear tests on 

this sand at different water contents. The sand in the box was placed at a 


wate content of about 5%, and the surface was covered with a plastic sheet 
retain the moisture. However, the plastic had to be somoves pe some — 
the 


—_— 
3.—Holography Model Facility 
water content decreased near the surface, the apparent cohesion if 


therefore, near the surface apparent co cohesion became the dominant ie 
providing confinement. ‘Using the cohesion intercept as effective confining - 


pressure to hes shear wave velocity the Hardin ‘equation, ped 
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wave velocities in the right column of Table 1 were obtained. At a water content 

¢ near 4% the empirical shear wave velocity was close to the measured Rayleigh © 
- wave velocity and a water content of 4% near the surface was very probable. _ 
Because of the shallow depth of the box, the small size of the footing, and 


— 
| 
| 
| 


_ the properties of the soil, it was not pnts to produce interferograms at 


lower frequencies which would have sampled deeper into the sand. It is expected, 
_ however, that the wave velocity would have started to increase with depth 
again when confinement from pressure became large compared with 


_ The major difference in . these nial. state tests and those in the field = 
that no transducers were needed to determine in-phase points for wavelength 
_ determination. Wavelengths were scaled directly from contours on the interfero- é 


and could be determined accurately. transducers, except 


TABLE —Direct Shear T Tests on Crushed Querte Send 


intercept, 


Water content, in pounds per 


asa percentage square inch | 


‘miniature accelerometers, could have been used with the short a 
CALIBRATION 
“ Scmenndiin (9) reported on an application of the University of Florida 
calibration chamber to study the effects of the state of shearing stress on the 
_ dynamic bulk modulus of sand. The senior writer participated in that study 
_ The calibration chamber is essentially a large triaxial cell- like chamber about 
1.25 m in diameter and 1.25 m deep (Fig. 6). The state of stress in the sand 
can be controlled by applying programmed vertical and horizontal stresses. _ 
_ Compression and shear wave velocities v were measured © on n planes at Selected | 


“confining conditions. These wave velocities used to compute the 

" low- strain, bulk modulus of the sand. To measure elastic wave velocities in 
the sand in the chamber, Schmertmann used miniature accelerometers embedded 
at selected locations during deposition of the sand, (Fig. 7). The accelerometer 

locations were Selected based on a the angry) of a access ports i in the top ‘Platen 


the ports in the top platen of the chamber. Using four accelerometers and _ 

three source locations, wave propagatiun was studied along four planes through rz 

velocities were determined by time f a to 


7 = of the chamber. Seismic energy : sources were operated through _ 


| 
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Shear wave velocity, 
(meters per second) 
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FIG. 7.—Location of Accelerometers and Sources in ) Calibration 
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FIG. 8.—Shear and Compression Wave Velocities in Calibration Chamber, Octahedral — 
~ Both compression and shear wave velocities were determined in planes at — 
_ Various orientations through the sand. An example of the results is plotted — 
in Fig. 8 in polar coordinates. In Fig. 8, velocity is measured radially from _ 
angle from the horizontal, 6, represents the orientation of - 


origin and the 


| _- between two transducers a known distance apart. The output from a set 
| _ transducers was displayed on a dual trace oscilloscope to determine the inte ; 
- - time. The accelerometers were located in the planes of interest and orie 
perpendicular to the direction of wave | 
 £ 
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The seismic om used in hese tests were essentially the interval velocity 
 eross- hole technique described by Stokoe and Hoar (12), expect that the receiver _ 
and source were at fixed locations rather than being moved from one elevation 

to another. ‘The major ‘differences between this: work and field cross- “hole 


Careful placement of the accelerometers ¥ was s also necessary t to assure me accuracy 
in spacing. Miniature input sources consisted of: (1) A 100-mm diameter plate 

_to which sand had been cemented to the bottom surface; and (2) a scissor-jack 4 

= which was expanded against the walls of a bored hole. Both torsion | 


Mope Founpation 1 Faciuity 


ilit 


m high and 2.3 m in diameter, containing a fine, poorly graded, _ 
‘round grained sand. The unique feature of this facility is that 
qj its contents may be reconstituted to any desired density by an appropriate — 
combination of liquefaction and vibratory densification. upper sand 
. Test bed preparation is initiated by creating a ‘ ae in the sand 
_ by an upward flow of water. Sand in a loose state is achieved by gradual 
‘ flow reduction after a high flow rate, and sand in a denser condition can be — 
_ achieved using an electric spud concrete vibrator at preselected depths and 
_ To establish material properties for use in the analysis of pile studies, Chon 
performed seismic cross-hole tests in the facility to determine the variation 
of dynamic shear modulus with depth. The crosshole apparatus used is shown 
in Fig. 10. The state of the sand as by a and = 
‘process was relatively dense \ with D, = 16% and e= = 0.6. Spacing between a 


generated with the scissor-jack packer. 
6 was described by Chon (4) and is shown 
q schematically in Fig. 9 This facility consists ofa vertically oriented, corrugated _ 
4 
wil 
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FIG. 10.—Schematic ot Hole Apparatus 
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1 Sheer Modulus Versus Depth in Model Pile Facility (4) 


| 


ase - rafts Mingle: 


4 the impulse and receiver t rods \ was 1.22 m and data were obtained from. depths -e 
7 of 0.61 m, 0.91 m, 1.22 m, and 1.52 m. The seismic source consisted of an sw 

Pa impulse rod coupled to an open-ended pipe. At each depth, the direction of — £ 
& the impulse was reversed to facilitate identification of the shear wave oll 
_ Shear wave velocities from th the cr cross-hole tests were converted to shear moduli . 


shear wave Fig. shear modulus versus "depth fi ‘from these 


_ The determination of soil properties near the surface of a cohesionless material 
by the cross-hole technique is very difficult because shear waves are not readily — 
_ transmitted through the lightly confined regions of the soil near the surface. 

Consequently, shear moduli were not assigned to the first 2ft(0.6m), 


; ae experimental moduli were compared to moduli computed from resonar 


by ‘Haupt excellent agreement was s achieved between the 
- cross-hole test, the resonant column tests, and the empirical equation. These 
results support the use of the cross-hole test as a means of establishing in-place — 
‘dynamic soil modulus in laboratory testing facilities. 


surf ace was unconfined and apparent cohesion did not provide adequate confine- 
ment. The junior writer undertook the task of correcting this deficiency | by 
adapting the facility so that an air pressure surcharge could be applied to the 
sand. The system is shown schematically (Fig. 12). Pure gum rubber air bags 

a were constructed to cover the entire surface area of the sand. The bags | 

Simultaneously reacted against the soil surface and a composite plywood and - 
4 expanded steel sheet backed by an I-beam grid attached to the culvert by stout 
: brackets at four points. Building supply compressed air was supplied through 


a aie leading to an manifold that ee air to each of the three air 


| 
int 
i 
Richart | used. As ff 


. manometer. It a found that at least 35 5 kPa 1 may be safely used 1, resulting a 
in a force of approx. 217 MN distributed uniformly over the area. 16.0% — 
7 new set of seismic cross-hole tests was performed by the junior writer 
with improved cross- -hole apparatus and techniques to establish shear = 
velocity as a function of both applied surface pressure and depth. The objective 
was to determine the surcharge pressure that produced the most favorable elastic 
property distribution for model footing tests. A schematic of om the apparatus 
its arrangement is presented (Fig. 13), 

The impulse source for the crosshole tests (Fig. 14) consists a rotating 
hammer supported by wooden bearings in a 64-mm diameter pipe. Attached 
to the upper end of the pipe is a target arm. The lower end couples to a 
rod, which in turn attaches to a Shelby tube. The rod is 25 mm in diameter 
and 2.44 m long while the Shelby tube is 355 mm long and 51 mm in diameter. 

. -_ Impact of the hammer Lecyttvay target arm arm applies a pure torque to the tube. © 

Wied 


SURFACE PRESSURE 


TRANSDUCER TRANSDUCER 
13 —Schematic of Torsional Cross Hole ‘Apparatus ai 


_ In an actual test, the Shelby tube is embedded 152 mm into > the soil otow 


: 4 Two horizontal ‘velocity transducers, placed at the same depth as the Shelby | 

; “tube and along a straight line from the source, were used as receivers. These 
geophones were oriented perpendicular to the array line for maximum response a 


(13). The hammer and the target ¢ arm n of the 1 torsional device are separated by a 
a wooden bearing which, in addition to its structural function, serves as an 
electric insulator between the two. As in Chon’s (4) tests, hammer-target impact 
initiates triggering of the oscilloscope sweep by closing an RC trigger ciecelt: 
_ The torsional device was excited in both the clockwise and counterclockwise 
directions, and clear shear-wave reversals were noted. No distinct trend, however, 


was ‘detected for the | compression wave because. of its low ar a and the 
Test were conducted at depths of 0.15 m, 0.31 m, 0.61 m, 0.91 m, 1. - 
 m, and 1.52 m. For each of these depths t the tank water level was vas slowly : 
- ‘ raised and lowered prior to testing in order to maintain a consistent water content. 

The source and transducers were then placed and the pressure was raised a 

& kPa to seat the air bags. Testing at a given depth consisted of obtaining 

_cross-hole data at each surface pressure of the following unloading-loading — 

sequence: 35 kPa, 28 kPa, 21 kPa, 14 kPa, 7 kPa, 14 kPa, 21 kPa, 28 kPa, a 
' q and 35 kPa. A 15-min delay was adhered to after every cm a change to 
_ allow for the minor time effects characteristic of cohesionless soils. 


3 


calculate wave the path distance receivers 
was divided by the time taken for a wave to travel from the near to the far 
transducer. Trigger error was found to be significant for the computation of 
velocity based on the first arrival at the near transducer; consequently, velocities 
were based upon interval timing only. Several points on the wave train were 

considered for wave arrival identification, and the results for a given impulse . 
averaged to arrive at the shear-wave velocity for that test. 2 

_ Fig. 15(@) shows the results of these cross-hole tests as a family of shear 

; wave moe versus depth curves with surface pressure as a third —,. 


| 


inst and t the variation of wave form with 


‘decreasing shear wave velocity with decreasing pressure may readily be observed 

by comparison among the records. Also apparent with decreasing pressure is 
‘the degradation of wave form believed to be associated with the inability 

the material to transmit shear disturbances at low levels of confinement. _ cs 

Bn of interest is the quality of the reproducibility indicated in Fig. 15(6) 
and the preservation of the wave form shape between receivers, especially s : 

at the larger pressures where the medium is most receptive to shear wave 
transmission and transducer Inspection of the | records 
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FIG. 15.—(a) Shear Wave Velocity Versus Depth for Model Facility; 


reveals a wave-form which is | in ‘each of the far transducer 

traces by a low amplitude disturbance, which is not present in the near transducer i 
- records and precedes the major vibration. This initial segment is believed to 

- be composed of compression waves that travel at a greater rate than shear — = 

waves, and reflections and refraction from a high wave speed surface or boundary. 

A final characteristic noted was a slight decrease in wave velocity with later — 

points on the wave train implying effects of either material nonlinearity, wave : 

_ dispersion, or reflection-refraction interference. This trend, however, did not %j ’ 
“ee — with depth. The greatest velocity variation among individual computations ai 
wave for one record was 10. of the average 


| 
in 


velocity for the record and occurred for the abet. test depth of 0.15 m 


7 at which influence of boundary phenomena may be expected to be large. At 
each depth, | the ‘maximum velocity variation determined and the smallest, 
-. 4.6%, occurred at the level of 1.22 m, 1, well removed from disruptive boundaries. ty 
- bh To compare with these experimental results, shear-wave velocity profiles as 
a function of surface pressure were calculated by use of the Hardin equation. _ 
; = void ratio was taken as 0.69. It was assumed that the air bags produced — 7 
a uniform vertical pressure equal to the applied pressure throughout the facility ia 
and K, was as assumed to be 0. 4. The confining stress distribution due to _ : 


- surcharge was ‘superimposed on that caused by geostatic stress. The results 
of calculations for surcharges of 7 kPa and 35 kPa and for geostatic confining — 
— only are also plotted on Fig. 1S(a), 4 
_ From Fig. 15(a) it is clear that an increase in the surcharge pressure results 
in an increase in shear-wave velocity at all depths. For each pressure the theory 4 
underestimates the wave velocity | near the > surface while overestimating wave 
induced stress in the vicinity of the surface, implying the existence of an arching - 

effect created by the tank walls. Studies which tend to confirm this idea are ab 

Am 1 currently being conducted in an effort to explain the experimentally observed 

| Additionally, these cr cross-hole tests results are now being used as a “basis | 7 
for redesign of the facility to avoid arching under the application of surcharge _ 
and further tests may be used to study the effectiveness of any future alterations. an 
_ ‘The demonstrated ability of the cross-hole test in the detection of small changes | 
in stress, as well as its potential for identification of a regional stress distribution 
which may have otherwise been difficult to obtain, is noteworthy j and suggests — 
‘solution of "the inverse problem, the estimation of the confining stress or a 
change in confining stress that produces an existing or a change in an existing 
= wave velocity distribution. Possible applications include the evaluation 


= pore water to soil during consolidation following construction, and the 
study of any process, which is associated with a change in confining stress. — 
_ By virtue of its scale, the execution of the laboratory « cross-hole test _—— 
a from its field counterpart. For short source to receiver spacings such — : 
as those required in laboratory testing, a given triggering delay may represent _ 
a significant fraction of the disturbance travel time, leading to large errors in 
shear-wave velocity computations. To avoid this difficulty, t the use of multiple — 
transducers permitting transducer-to- to-transducer- based she shear-wave velocities ‘is 
—— essential. Because receiver- to-r -receiver distances are are small, wave form 
spreading, which has been observed in the field by Stokoe and Hoar (12), is — 
te particularly significant in the laboratory. This permits shear-wave velocities 
to be based upon the average of velocities obtained from different wave points — 
appearing within the initial portion of the wave train. The advantage of this — 


procedure is that random shifting of peaks, due i in part | to boundary reflections — 


and refractions and in part to compression ‘wave separation, “may be averaged | 


Seismic are used in the field can 


_af damace ta earth macces fallawine earthanakes monitoring of the stress t 
| 
| 


a ~~ (1) Size and energy content of s seismic source; > (2) receiver size and sensitivity; ; 


are advantageous in laboratory, namely; the tests are nondestructive 
can be performed on materials in situ. Sow, Brow 

The seismic cross-hole tests proved to be very effective means for sre-wenall 

the spacial vi variation n of elastic soil Properties ir in 1 large scale laboratory | facilities. 


confining pressure spr voting a characteristic of considerable practical value for 
estimating effects of facility modification, and for the preparation of the most 
favorable material property distribution for foundation 
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DRILLED-SHAFT INTEGRITY BY WAVE 
PRopAGATION METHOD* 
‘By Thomas M. Hearne, Jr.,' Kenneth H. Stokoe, II,? Associate Members, 
ASCE, and Lyndon C. Reese, F.ASCE = 
use se of cast-in- -place concrete shafts and the difficult soil 
conditions under which “many of these deep foundations are constructed, a : 
uncertainty often arises concerning the integrity of the cast members (3,19). 
Foundation performance may be greatly impaired if defects or irregularities _ 
exist in these shafts. Defects and irregularities such as ; voids, soil inclusions, 
"separations, enlarged re regions, or necked regions may occur during construction _ 
(6,13,14). Therefore, definitive and economical test procedures are needed that 
can be used shortly after construction to evaluate the structural integrity of 
drilled shafts. Then construction can proceed with So or the eo 


a recent of the kinds of that come to the 

; - have settled about 4 in. (100 mm), even though the full load had not been 
rs applied to the structure. While there were some unusual features in the construc- 2 
4 tion procedures that led to the likelihood that there was excessive settlement 


the possible lack of integrity of the drilled ‘shafts at the | site. Tt is not 
7 ‘Man proposed, of course, that every drilled shaft be instrumented, but it is ’ 
- suggested that strong consideration be given to instrumenting routinely a certain 
d fraction of the drilled shafts on many construction projects. Had instrumentation 


been used at the bridge s site, and instructive have 
been gained at a small cost. 


“Presented at the October 22-26, 1979, ASCE Annual Convention and Exposition, held 
my Research Asst., ‘Dept. of Civ. Engrg., Univ. of Texas, Austin, Tex. 78712. oie soars 7 
Assoc. Prof., Dept. of Civ. Engrg., Univ. of Texas, Austin, Tex. 78712. 
Prof., Dept. of Civ. Engrg., Univ. of Texas, Austin, Tex. 78712, = 
_ Note. Discussion open until March 1, 1982. Separate discussions should be submitted © 
an for the individual papers in this symposium. To extend the closing date one month, _ 
a written request must be filed with the Manager of Technical and Professional Publications, 
ASCE. Manuscript was submitted for review for possible publication on February 17, 
1981. This paper is part of the Journal of the Geotechnical Division, 
the American Society of Civil Engineers, 
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in ¢nd Was no way to prove the settiement not OCCur 


the structural integrity of drilled shaft foundations. These methods include nuclear — 
radiation, sonic, stress wave propagation, , and vibration methods (1,4,5,7,11,12, 
& 17,18,20). Their application has met with varying degrees of success. During, 
the last 2 yr, the writers have s successfully adapted and employed the stress — 
wave propagation method to evaluate drilled-shaft integrity. In theory, this method 7 
is quite simple; that is, one just measures the time for compression =. 
to propagate vertically through the shaft from a source to one or more receivers. © 
method a viable field which is described 
Of the many different source-receiver sills that can be employed | 
with the compression wave propagation method, two basic configurations have 
been used in this study. In the simpler of the two configurations [which has | 
been used with limited success in the past (1,11,12,17,18)], both the source 


. and receiver are e installed on the top of the shaft. The time required for the | 


m 
and then return to ‘the top again is measured. _ This test configuration i is refered * 
to herein as compression wave propagation with surface receiver,  — ; 


In the second source-receiver configuration, which was developed in this & 


: study, the source is installed at the top of the shaft, and the receiver is mounted 


within the shaft. The receiver is usually mounted on the steel reinforcing cage 


ata predetermined elevation. With the receiver position known, the source-to- -re- 


7 ceiver travel distance is also known. Therefore, P- -wave pro opagation velocity 

can be determined directly once the source- to-r -receiver travel time is measured. ; 

__ This test procedure is referred to herein é as compression wave propagation with © 


Regardless of the source-receiver configuration employed, successful use of 


_ the method requires identification of reflected wave arrivals so that travel times" 


a of reflected waves can be used with propagation velocity to determine « elevations 
7 of reflection points. Reflection ‘points for sound, straight-sided shafts are the 
- concrete-air interface at the top of the shaft and the concrete-soil interface — 


at the bottom of the shaft. For shafts with discontinuities and irregularities, 


_ additional reflection points occur, and it is the monitoring of reflections from 
ae additional points which ‘permits identification of defective or irregular 

sKafts. The amplitudes of any monitored reflections depend upon many factors” 
” including the input signal amplitude, concrete quality, surrounding soil conditions, __ 
.~ size and nature of the discontinuity or irregularity, and the location of ¥ 

the receiver relative to the elevation of the discontinuity or irregularity. = 
me: _ The following formulae are used in the compression wave propagation method 


to calculate the P- ~-wave © velocity, Ve, and and Young’ s modulus in a rod, E, for 
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DRILLED SHAFT INTEGRITY. 

= travel distance; ¢ = travel time; y = . unit weight of concrete; “ee 


i at low strain levels. Strain levels | ‘resulting from P-waves propagating down 
the shaft are typically less than 10~* mm / mm. Asa s modulus 


at higher strain levels, such those determined in compression 


The most successful application of the compression wave propagation method _ 
in this has ‘esulted oom the use of embedded receivers. The 1 receive 


= 
FIG. 1.—Compression | Wave rie FIG. 2.—Determination of ‘Direct and 
are cast in shaft t preselected depths during the construction 
operation. Each receiver consists of a vertical velocity transducer sealed 
an epoxy-coated aluminum case. An electrical cable extends from the aluminum Ka 
case to about 10 ft (3 m) beyond | the top of the shaft. The receiver cable ‘, 
ay usually taped to the reinforcing steel along the shaft length to permit | ease a 
a placement during construction. The receiver case i is joined to a metal bracket 
which is in turn attached to the reinforcing steel with steel hose clamps. A s 
silicon rubber pad is used to isolate the bracket from the reinforcing steel. A | 
In these tests, isolation of the receiver from the reinforcing steelis animportant | 
woud onsi ion because the compression wave velocity in steel is greater than J 
concrete. If significant levels of compression wave” energy are 
_ within the reinforcing steel and isolation is not provided, then monitored travel 
may correspond to the compression wave velocity in steel rather than ff 


= 
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hand- held hammer and chisel); (2) an electrical trigger; er; (3) a vertical receiver(s); 
@ a storage oscilloscope (or digital oscilloscope); and (5) a Polaroid camera. 4 
_ Compression waves are generated at the top of the shaft by releasing the a 
_ drop hammer from a predetermined level. Hammers used in this study weighted © , 
from 10 lb-32 lb (4.5 kg-14.5 kg) and drop heights ranged from 6 in.-30 in. — 
(150 mm-760 mm). An alternative approach for generating P-waves is to replace 
the drop hammer arrangement with a hand-held hammer and pointed chisel. 
However, the drop hammer technique provides a constant input amplitude which _ 

_ is desirable when analyzing the amplitude of monitored P-waves. 
When the falling hammer contacts the steel plate, the oscilloscope is triggered, 
“and recording of the receiver output commences. The time for the compression ; 
_ wave to travel from the top of the shaft to the receiver is then recorded. ‘This 
a is called the direct travel time, and, with the travel distance known, one all 
can compute the direct P-wave velocity. 

For accurate timing measurements, it is essential that triggering of f the 
oscilloscope. occurs at a known time after the hammer contacts the steel plate. 
If triggering is delayed some unknown amount, then the direct velocity measure- _ 
ment will be in error on the high side, with the magnitude of the error dependent 
upon the delay time and the travel distance. a eco, 
__ One means of eliminating triggering problems is to install two or more receivers _ 
with 2 vertical spacing of 10 ft (3 m) or more. Measurement of an interval 
travel time can then be performed. The interval travel time is the time for 
he wave to Propagate from one receiver to another. Such a measurement is 
_ Shown in Fig. 2 and is based on the initial arrivals of the P- wave, ‘denoted a 
by the P’s in the figure. The output of both receivers for the same impulse - 
_ has been captured on the oscilloscope record, and any time delay in triggering : 
_ did not affect the interval travel time measurement illustrated. Use of a vertical 7 
spacing of 10 ft (3 m) or more between receivers minimizes the influence of 
 f slight variations in interpretations o of the initial P- -wave arrival times. I ‘However, — 
any time delay or variations in interpretations of P- -wave affect 


"which can complicate interpretation of P-wave arrivals. Rayleigh waves (surface 
waves) are also generated by the hammer impact. These surface waves produce | 
- very n noisy environment in the upper portion of the shaft [most severe in _ 
the u upper 10 ft (3 m) in this study]. Therefore, embedded receivers are placed — sal ; 
at a depth of 10 ft (3 m) or more to minimize any unwanted monitoring of _ 
_ Rayleigh waves. The noise generated by surface waves constitutes a major 
problem when using a surface receiver(8).  j§.§ 
- To monitor wave arrivals accurately, receivers must be used which track 7 
closely wave particle motion. Frequencies generated in concrete by mechanical 
Hi ‘sources vary from several hundred to several thousand Hz (cycles per second). & 
‘Receiver packages utilizing small velocity transducers with natural fi requencies ~ 
between 4.5 Hz and 40 Hz are economical (about $40 each, including cable) 
and can be used satisfactorily inconcrete. = j= 
Vertical transducers are utilized to e the vertical particle motion 


Used in recordin -wa 
“ie 
4 
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q these receivers, one must ensure a good bond between the. concrete and receiver. 
Also, isolation of the receiver from reinforcing steel | must be provided 
. _ A storage oscilloscope and Polaroid camera are well suited | for monitoring — 4 
receiver output and obtaining a permanent record at a reasonable cost. More 
precision and greater ease in data analysis can be obtained by using a digital : 
oscilloscope with magnetic recording capabilities. However, the digital oscillo- 
i ; scope is about two to three times more expensive. An electrical trigger is used 
to initiate timing with either oscilloscope. The trigger is essentially instantaneous : 
when properly used and typically delays the oscilloscope sweep by less than 
0.01 ms. A complete description of the trigger and possible triggering postions 
is given Hoar and Stokoe (9). € 


(DA compression wave source “pe drop hammer with steel plate or hand- held 
hammer and chisel); (2) an electrical trigger; (3) an accelerometer; (4) a charge _ 
conditioner (used to empty the accelerometer signal); (5) a filter with 


& uration. The procedure involved in generating and recording the compression 


Pa wave is the same as that described for compression wave propagation with 
De The monitoring equipment and testing configuration in this test procedure _ 

te from that when testing \ with embedded receivers in several respects: (1) * 
The receiver is bonded to the top of the shaft with epoxy rather than being 
embedded at some depth within the shaft: (2) a small accelerometer is used } 
_ the surface receiver instead of a velocity transducer because reflections are 
‘more easily monitored with the accelerometers (velocity are used 


output is typically filtered to ‘reduce n noise ¢ generated at the top of the shaft. = an 
__ In terms of convenience, economy, and ease of testing, the procedure using 
a surface receiver has obvious advantages over testing with embedded receivers. 
One receiver can be used to test all shafts, and no instrumentation need be 

before or during shaft construction. However, these advantages are 
overshadowed by several severe limitations. 

RR The success of the method requires that the reflection from a bottom _ 
of the shaft (or a defect or irregularity) be monitored. If a reflection cannot — 

be monitored or identified, the results of the test are inconclusive, and no : 
accurate information can be obtained to judge the concrete quality or integrity. — 

2. P-wave velocity cannot be measured directly. A method exists for estimating - 
_ the P-wave velocity by monitoring the frequency of the surface wave reflecting 
a _ back and forth across the top of the shaft (18). By rememarie this frequency, q 


= 


|. 
7 an estimate tne compression wave velocity in the shaft. However, at dest 


i is a crude estimate. Drilled shafts encountered i in often have highly 

surface shapes, thus rendering this procedure questionable and many 
3. The most important limitation of testing with a surface receiver is that | { 
any impact at the top of the shaft generates surface waves which reflect back __ 
and forth across the top of the shaft for a time period that extends well beyond - 
P-wave reflection arrivals for most shafts. These surface waves act as noise, — 
which complicates and many times precludes identification of the arrival of 

_ the reflected P-wave. This noise problem is further compounded ea - reinforcing 
cage or w steel dowels which extend above the shaft surface. 


HAMMER 


TRIGGER 


FIG. 3—Compression Wave Propagation Method with Surface 
Filtering may be used to minimize surface noise. ‘However, filtering may 
a - cause additional problems, many times unrecognized problems, by altering “ 
amplitudes and travel times of reflections, and one must be very careful when 
. filtering receiver output. If the generated P-wave frequency is close to the 
_ _ surface wave frequency, filtering will not be advantageous because the P-wave 
‘ will be attenuated along with the surface noise, resulting in no =e 


in the monitored receiver output. As the diameter of the shaft decreases, the 
‘frequency of noise created “ae the surface waves increases . Filtering is most 
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“DRILLED SHAFT INTEGRITY : 


sources on this or smaller, the generated compression = 

é frequency differs from the frequency of noise generated by the surface wave, 
and a filter can be used effectively. 
In this study, P-wave velocities in 3,000 psi-4, 000 psi (20,700 kPa-27, 600 
kPa) concrete have ranged from 12,000 fps to almost 15,000 fps (3,600 m/s—4, 570 
m/s) for measurements made 2 or more days after casting. Compression wave sai 
velocity cannot be directly related to concrete strength due to mix design — 

variations. However, for a given mix measured ~wave velocity 


measurements as made with embedded receivers are very worthwhile. * 
oul A suggested rating of concrete quality based on P-wave velocity is lle 
in Table 1 (10). This rating is based on ultrasonic pulse tests (high- frequency 
4 tests) which were used to measure P-wave velocity through a concrete medium. a 
In wave propagation testing, the compression wave velocity” in a rod differs 


Compression wave velocity, ‘i 
(meters per second) condition 
above 15,000 (4, 570) 
12 {000-15,000 (3,660-4,570) good 
-“4 10,000-12,000 (3,050-3 660) } 
(2, 130-3,050) poor 


rom the compression wave velocity ina continuous medium of the same ‘materia ' 


a rodlike structure, the wave propagation velocity as determined herein will 
be less than the velocity r measured i in a concrete medium at ultrasonic frequencies — 
(15). This velocity reduction is ; approximately 10% for a Poisson ratio of 0.25. bs 
- Rae, velocities shown i in Table 1 should be reduced by 10% if — 
ay) 
Shaft lengths can be computed from reflection arrival times for measurements _ 
with embedded receivers. These lengths should correspond closely to the actual = 
, shaft lengths for sound shafts without irregularities. The success of this technique — 
is, however, based upon the ability to measure or estimate P-wave — 
and to identify at least one reflection arrival. wend 
When the P-wave encounters the concrete-soil interface | at the bottom of E] 
the shaft, some incident wave energy will be reflected and some of ia. 


4 
i 


‘eer wave energy will be transmitted into the soil nenostie’ the shaft. If a sighifiens 
portion of the energy is reflected, the P-wave will be detectable at the embedded 


receiver. Under these circumstances, the possibility also exists of identifying — 

_ an additional reflection from the top of the shaft where essentially total reflection of 
“4 of the energy occurs. In this case, the 2L travel time can be determined as 
- _ shown in Fig. 4, and an estimate of the effective shaft length, L, can be made ‘oad 


A typical record | showing the propagation of the compression wave as outlined ; 


t 


with three receivers embedded in a large [48-in. . (1, 120-mm) diam] shaft at 
; the depths indicated. The 2L travel time for each of the receivers is the time _ 
- interval from the initial compression wave arrival to the second reflection arrival, 


Sweep 
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Rb-Refiection From Bottom of Shof 
Reflection From Top of Shoft 


Rt, which is from the top of the shaft. In this record, the 2L travel, time could 
be determined from successive reflection from the base, 
_ It is interesting to note in Fig. 5 that the first reflection has the same polarity — 
as the initial arrival, indicating that the compression wave | is reflected as < 
“tensile wave from the base of the shaft. This tensile wave is then reflected 
aS a compression wave from the top of the shaft, and the process repeats. — 
‘This represents a typical reflection pattern for sound drilled shafts. 5 
Typical shaft length computations based on reflection travel times for 4 
te _ with embedded receivers are shown in Table 2. The listed shafts are representative 
_ examples from three different sites in Texas. An example of a record from 
4 a shaft with a known defect is shown in Fig. 6(a). During construction, a 2- # 
©. 6-m) thick clay seam was added near the middle of the 39-ft (11.9-m) a 
\ shaft to simulate possible defects such asa soil inclusion | or a eam 


F 
| 
FIG. 4.—Reflection Analysis with prd for 92-ft (28- 
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‘ben 2 —Calculated § Shaft Li Lengths Based. on on Reflection Travel Times 


ghaftlength, —onreflections, 


number 


in feet (meters) 
tess 0 


9.9 28.0) 


98.1(29.9) 


28.8) 


11.9 @G4.1) 


(7. 


in feet (meters) 


103.6 31.6) 


102.0 (31. 


113.5 (34. 6) 
44.7 (13. 
(14.1) 


25.0 (7.6) 


Percent — 


difference 
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@ Receiver Location 
Direct P-Wave Arrival 
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Reflection From vem Yep of 


«Shaft Length=39 (119m) 
Clay Defect = 18.5 ft (5.6m) 
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covering ng the entire cross section of the shaft (1). The shaft i ii insnnemiented “a 
with three embedded vertical velocity transducers, as shown in Fig. 6(b). 
_ Several important points that highlight inherent strengths and limitations < of 
P-wave propagation testing are shown in this record. The receiver located at 
a depth of 5.2 ft (1.5 m) monitored two distinct reflection arrivals much earlier 
than would be anticipated for a sound, 39-ft (11.9-m) long shaft. By using the _ 
2L travel time determined from this receiver and a P-wave velocity of 12,200 
ft/sec (3,719 m/s) determined from the interval travel time between the top * 
- two receivers, the effective shaft length, L, is 18.3 ft (5.6 m), indicating a a 


Receiver 


Rate= 
ms/div. 


div. 


Reflection From Bottom of men 
From Top ot Shalt 


Direct P-wave Arrival 
® Reflection from Top of Shaft ; 
Reflection from Bottom of Shaft 


Shaft Length =112 (341m) 


Shaft Diameter = 36 in. cm 


7 7.—Reflection Reseed for 112- (34-1 m) Long Shaft with Defect 


_ The record from the receiver at a depth of 18.2 ft (5.5 m) can be interpreted : 
two | ) ways if considered by itself. First, the 2L travel time (as | shown in Fig. . 
4) may be interpreted to be 6.14 ms, | which gives a shaft length of 37.5 ft 
(11.4 m). This length is close enough to the actual length so that it is possible 
a 0 misinterpret the record as being representative of a sound shaft. This occurs 
because the receiver is adjacent to the defect and is located near the middle = 
of the shaft length. However, the amplitude of the first reflection arrival monitored 


e this re receiver is es than ‘normally encountered for similar sound shafts, 
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a reflection point, there is only one reflection arrival (from the top of = 
shaft) for each 2L_ travel distance of the P-wave . Therefore, th the second 
: interpretation is that the 2L travel time is 3.00 ms, which is represented by 
the time interval between successive troughs on the record. This interpretation 
gives a length of 18.3 ft (5.6 m), which correctly identifies the location of 
7 the defect. This second interpretation would automatically occur when the output — 
from this receiver is analyzed with the other embedded receivers. 


The bottom receiver located at an elevation of 38.5 ft (11.7 m) recorded e 
avery weak P-wave arrival, and it is obvious that some type of defect covering 


a large percentage of the cross section is present. (The sensitivity ‘setting for 
: the bottom trace is 50 times greater than that of the upper two traces.) Additionally, 
_ the measured interval P-wave velocity from 18.2 ft-38.5 ft (5.5 m-11.7 m) is 
reduced by almost 2,000 ft/sec (610 m/s) from the interval velocity measurement 
between 5.2 ft (1.6 m) and 18.2 ft (5.5 m). This observation further indicates 
a defective zone exists within the shaft. However, ‘essentially no information 


is present in the bottom trace, which can be used to determine ‘exactly where 


Anexample of a record froma shaft with a possible defect which was discovered | 
during routine field testing is shown in Fig. 7(a). This 112-ft (34.1-m) long 
shaft was Cast by the slurry displacement method and was one of many which 

were tested on a large construction ‘project. The shaft was not cored, however, 
because it was not going to be used in the structure. Fig. 7(b) shows the possible 
— paths for the P-wave record shown in Fig. 7(a). A distinct reflection 


Ri occurs in addition to the reflection monitored from the bottom of the shaft 


(Rb). The reflections are most easily seen on the bottom trace of the record. — 
This output was obtained from a receiver which was located near the middle 
of the shaft length, and reflection arrival times should be approximately — 
spaced across the record. However, two distinct reflections are monitored within 
avery narrow time frame. Because the direct wave amplitude is not significantly — 
reduced from the monitored amplitude above the region in question, the —s 


However, the ‘Teflection from “the bottom of the shaft is not recorded | by the 
: top two 1 receivers, thus n minimizing this possibility. The most likely pees 
were determined to be a severely necked region or a weakened zone of concrete 
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ty 


Attenuation of the P-wave as it propagates down a shaft is influenced by 


ee several factors including concrete quality, the presence of irregularities or 
discontinuities, surrounding soil conditions, and method of concrete placement. 
(Olin concrete may result in segregation which decreases strength and 
i increases attenuation. Drilled “shafts s surrounded by hard material have - 
- found to exhibit higher P-wave attenuation than similar shafts in soft soil. _ 
’ ‘Fig. 8 shows P-wave attenuation in a 92-ft long (28-m), 48-in. (1, 120-mm) - 


; diameter drilled shaft founded in clay and cast by the slurry displacement metho 


d. 
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“The shaft was in 8 ft Q. 4 m) of water with the top 10 ft (3. 0 m) of the shaft 
_ above the mean water line. At the bottom of this shaft a portion of the compression pe 
_ wave energy was transmitted into the underlying soil. Theoretically, for a plane — 


wave, the ratio of reflected-to-incident compression wave amplitude depends 
on the density, and P-wave velocity contrasts between the concrete and soil _ 
at the base of the shaft (15). This ratio is expressed as follows: Hiaenhig 


‘in which ¢ C= seflocted compression wave amplitude; _A = incident compression 

wave amplitude; p, = mass density of concrete; p, = mass density of f underlying — 


soil; Va= = P-wave velocity in concrete; and V., = P-wave velocity in nunderlying — 
at these properties of the soil at the base of the shaft are similar t to those _ 
concrete, the reflected wave will be small. The greater 


Yet 


300 
Bottom Retiection 


in. (1, 1,120-mm) Diam Diameter Shaft P-Wave -Wave 


the contrast in the properties of the two materials at the reflection interface, 
the gre ‘greater will be the ratio of the “reflected- -to- incident compression wave 
7 ; amplitude. Therefore, almost all of the incident compression wave amplitude 
is reflected from the top of the shaft at the concrete-air interface. The theoretical 
relationship i is shown in Fig. 9. « ode 
To monitor accurately P-wave attenuation and reflection amplitudes, several 1 
d embedded receivers are required in each shaft, an approach usually not economi-— 

: cally feasible. However, the reflected-to-incident amplitude ratio at the bottom 
of a shaft could be closely estimated by monitoring the output from a single 
4 receiver embedded in the lower portion of the shaft. With this —— 

= gross indication of fics competence of the ‘material immediately underlying 


4 oneas is no single “best” location for placement of a receiver for all conditions. 
.. The best Teceiver location or locations ‘Should be Selected | d by considering the the 


‘ 
_ 
q 
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following i items: (1) Stiuctentien importance of shaft to be instrumented; (2) receiver _ 
costs; (3) existing soil conditions; (4) construction method; and (5) anticipated — 


very Unfiltered Output From, Surtoce “Accelerometer ols ay ws 


Rb-Reflection From Bottom of Shoft 
6 Filtered Output From Surface Accelerometer 
Reflection From Bottom of 


Untiitered Output From Embedded Velocity Tr 


FIG, 10—Surface and Embedded Receiver Records for 25.2-ft (7.7- -m) Long, 18-in. 

be Likely problem r¢ regions where embedded receivers should not be placed include z 
the following (based on in wavelengths ‘generated in this study): : (I) Withia a 10 
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ft G.0 0 m) of the top of the shaft a noise is a major problem; (2) near 
the elevation of the bottom of the nese if the casing method of constructio ed 


: (3) of on short reinforcing « cages s [approx 10 ft GB. 0 m) long or r less} because of 7 
noise problems; and (4) at the bottom of the shaft or in the proximity of possible F 
locations of defects or irregularities because wave arrivals and reflections will 
overlap, which makes interpretation difficult if not impossible. = | 
_ If concrete quality is the major concern at a job site, it is recommended : 
- that a receiver be positioned about 1 wavelength [approx 15 ft (4.6 m)] from 
the bottom of the shaft. It is important that the receiver elevation be recorded 
as accurately as possible. Direct velocity measurements can then be made to 
the depth of the receiver. A relationship between P-wave velocity and concrete a 
strength can also be developed for a given mix design. Then, by maintaining — 

a a log « of velocities, shafts $ with problem concrete can be easily ide identified _ based % 

7 _ Shafts naa than 2 ft (0.6 m) in diameter that are not founded i in hard material — 
can be tested by using a surface receiver. However, even under favorable — 


Recewer — Receive 
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«FIG. “11.—Reflection Record for 10-ft 12. —Retlection for 58-ft 
(3-m) Long, 12-in. (300-mm) Diameter, (17.7-m) Long, -in. Diame- 


- conditions, some of the tests are likely to be ate, On the other r hand, 
even under the most adverse conditions, testing with an embedded receiver ‘ 
near the bottom of a shaft will give the direct P-wave velocity as a minimum. ‘ 
If the bottom reflection is monitored, the shaft length can be determined. 
_ Additionally, P-wave amplitudes and other reflections may be monitored to 
aid in determining the locations of irregularities and discontinuities. Ages 
If defects such as soil cave- -ins and entrapped ‘drilling mud are + likely to be | 
- encountered, it is most advantageous for a ra receiver to be placed n near the middle 7 
of the shaft length if soil and construction conditions permit. Reflections within 
sound shafts without irregularities are most easily monitored in this region. —__ 
Three oscilloscope records obtained from a 25.2-ft (7.7-m) long, 18-in. (460-mm) 
- shaft tested with both a surface receiver and an embedded receiver are a 
shown in Fig. 10 to show the difference between the two methods. The unfiltered ¥ 
_ oscilloscope record obtained from the surface receiver (acclerometer) is shown 
7 in Fig. 10(a). This record is typical for drilled shafts with irregular surface — 
shapes and with reinforcing steel extending above the top of the shaft. As 


Fig. 10(a) shows, surface noise masks any reflection arrivals. The filtered 
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in the shaft pa be determined by dividing | twice the shaft t length by the time 
_ The unfiltered receiver output obtained from an embedded receiver at a depth 

; “of 16.5 ft (5.0 m) is shown in Fig. 10(c). Some noise is present, but two reflection — 
j arrivals can be identified from this record. With the direct travel time measure- 
ment , the compression wave velocity in the shaft can also be accurately 
determined, and the shaft length can easily be 


J xd ‘Carreatly, sev several J testing limiations exist for the compression wave propagation 

; 3 procedures described herein. For a shaft with several defects or ‘irregularities, , 
it is very likely that by using a surface receiver, only the elevation of the 

_ uppermost defect will be determined and defects at lower elevations will not 

be identified. Also, even when embedded receivers are used, the extent and — 7 
severity of defects cannot be accurately determined in most cases. For example, 

thin soil seam across ‘the entire shaft cross section might have the same 

reflection characteristics as a soil inclusion co’ covering about half the shaft cross 

section. However, if an embedded receiver is located beneath the soil sien 

an indication of the seriousness of the defect might be determined from reductions 

_ in the monitored P-wave amplitude and direct P-wave velocity. ae a 

0 The success ss of both ‘methods requires reflections to be monitored and i identified. 


Often testing techniques are developed under conditions which do not 
field conditions, such as on model piles Of precast, concrete piles lying in 
construction yards for the techniques presented herein. Results under these — 


unrealistic conditions are ¢ falsely encouraging because reflection ai arrivals are re easily 


surface shapes and protruding reinforcing ‘steel are absent. Fig. ll shows a 
record obtained from a receiver embedded in a 10-ft (3.0-m) long, 12-in. (300-mm) : 
diam free-standing concrete shaft. Reflection arrival amnpleates are high and» 
damping is low which results in an incorrect prototype test. Fig. 12 shows 
a more representative record. This record was obtained from a receiver embedded 
in a 58-ft (17.7-m) long, 24-in. . (610-mm) diam drilled shaft founded in a a clay a 
deposit. Only two reflection arrivals are monitored, and reflection amplitudes 
are strongly attenuated and are easily masked by noise. Of the many factors | 
affecting compression wave attenuation, the surrounding soil is a major contribu- 7 

im Using the methods described herein, generated wavelengths are approx 10 © 
-ft-25 ft (3.0 m-7.6 m). With these wavelengths, defects on the order of a small 

a4 ‘portion of the wavelength will be undetected. One improvement in the technique 
~~ be to generate higher frequency P-waves (thus shorter wavelengths) to 


( 
sample effectively smaller regions of concrete. However, wavelengths cannot ~ ! 
be allowed to decrease to much less than the shaft diameter or the shaft will ,. 


no longer represent a rod-like structure. Instead the shaft will behave as an © 


elastic medium, and P-wave reflections will originate at all shaft boundaries. i 


is shown in Fig. 10(b). Filtering eliminated the surface noise, and two —_ 
| 
| 
| 


$0} 


techniques herein, the tests ‘geovide an “economical and rapid means means 
_ of assessing the overall structural integrity of drilled shafts. When embedded - 
peceivets are used, sound shafts without irregularities are any identified 
are shafts with discontinuities such as complete separations and large voids” 4 
that span a major portion of the cross-sectional area of the shaft. However, ’ : 
testing procedures and interpretation methods need to be improved so that the 
full potential of the methods presented herein can be utilized. “ae We rede 
4 _ Of the two compression wave propagation techniques reviewed herein, the | 
= most successful application has resulted from the use of embedded receivers. 
; The advantages of using one or more embedded receivers instead of a surface : 
_ receiver include the following: (1) Direct and interval P- “wave velocity measure- 
ments are re performed; (2) two reflection arrivals occur for each 2L travel distance; 
@ lower noise levels are encountered; and (4) attenuation is measured directly. — 
_ The major disadvantages are costs of the embedded receivers and the installation — 
required prior to and during construction, 
_ Advantages of using a surface receiver instead of embedded receivers include wd 
and ease of placement since only one receiver "is required for 
phase. The disadvantages of ‘using a surface receiver include the following: 
(1) Direct P-wave velocity measurements are not possible; (2) testing is inconclu- — 
= if reflections are not meme and only one reflection arrival i is mneqpared ra 


by The value of propagation ‘testing asa been clearly 
- demonstrated to the writers during the performance of this study. With only 
one exception, drilled shafts that have exhibited irregularities in testing with is 
: embedded receivers have given no Mitteation by visual inspection during con- 
_ struction that any problems existed. In addition, many shafts which were cast 
under difficult conditions in which delays and other construction problems 
_ This work represents part of the first writer’s doctoral dissertation. The work 7 
was supported by the Bureau of Engineering Research at The University of 
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AppeNvix ll—Notation or fos fet | 
on The following gymbols are used in this pa paper: va. fs 
= incident compression wave amplitude; — 
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reflection from bottom of shaft; 
= reflection from irregularity; AAT 


_ feflection from top of shaft; rime 

compression wave vi velocity i in concrete; 


OV 


= compression wave Ni in underlying soil; 
weight of concrete; 
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OF SPT-N AND Revanive D DENSITY 
Michael N N. Fardis,' A. M. ASCE, and Daniele 


“The | Standard Penetration Test (SPT) hes | ‘for a number of 
4 years as amethod of indirect evaluation of subsurface soils. The test is particularly — 

~ useful for granular cohesionless soils, e.g., sands, for which it provides a measure 
of denseness. In fact, laboratory experiments indicate that for a given sand 
and vertical effective stress o,, the blowcount number N is highly correlated 
» 
with relative density D,; N is a measure of “strength and for a given ‘sand, 
_ Strength increases with D, and o,. The fact that measuring N in the field + 
_ less expensive than measuring D, has encouraged development of relationships 
N, D,, and G, (3,4,5,8,10). These relationships are dependable under — 
laboratory-controlled conditions but their r validity in the ; field is reduced due 4 
to the effect of factors, such as soil structure (fabric), consolidation history, 
inference of D, from N under both | laboratory and field conditions. The model 
in the laboratory (‘‘reference model’’) is somewhat more sophisticated than 
those presently available i in that it accounts | for errors that have been previously 
- neglected, specifically, those in the determination of D,; it is also based on 
a more comprehensive set of data. Finally, the log-linear form of the proposed 
& relationship offers computational advantages, especially when one considers 
_N and D, as random variables and wants to derive probabilistic characteristics 
Only some of the features of the laboratory model are assumed to hold in 
the field; other properties, such as trends and probabilistic correlation in space, 
are site-specific and must be inferred from site data. 
" Once a site-specific relationship has been constructed by appropriately modify- 
ing the reference laboratory model, one can infer properties of in-situ D, from 
i field measurements of N. One such relationship has been used in Ref. 2 to 
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infer probabilistic characteristics f Dasa random field i in space from charac- 
_ teristics of the waton field, N, and to evaluate the hazard of a 


actual values of N, D,, ¢ etc., and measurements or estimates in the Laboratory 
(N.,D,, _ or in the field (N,,D,, ...). Some of these quantities, ¢.g., N, 
<q do not exist prior to the experiment “that measures them; others, e.g., D_, are 
Se a in terms of a mathematical idealization of the soil material. In both 


the field be 6,, .. .. (here 0, stands for either 6,, or 6, with subscripts, 
L and F, denoting ‘ "mabermaty and ‘‘field,”’ ” respectively) and let 0, by a point 
of @ from @,, .. 8, = =, 6,. Then by 4...) 
of we mean the limit 
lim 6, 
provided t this limit exists in some probabilistic sense. As. defined by ae: a 
- the actual value of 6 depends on the type of experiment and possibly on the 
form of the estimator 6,. Of course, indeterminacy disappears as reference 
is made to a specific experimental procedure and a specific estimator. 
Retationsuips Between N, D, ano G, 


Regression Models.—For normally consolidated sands, laboratory information 


on N, D, and a, suggests a  log- -linear regression model of the type 


in 1 which N, and denote the actual values of the associated parameters. 
of-fit term measures departure of In N from a, + a, ln D, + 

InG,. It is modeled here as an uncertain function of D, and 6 G, with marginal in 
7 distribution NO, a?) (normal with mean 0 and variance o?). Eq. 2 fits data 

in the literature (3, 4 ,5) at least as well as the Gibbs and Holz family of curves 7 
‘ (3) or the recently developed WES nonlinear relationships (4,5), and offers — 

_ Unfortunately, standard regression techniques cannot be i to estimate 
“the parameters a,, a, and a, in Eq. 2 from laboratory data on (G,, D,, N) because 
of agew in the measurement of the last two quantities. | Uncertainty on N, 


given the actual value N is easy to treat, especially if one can assume that 


_ estimation cred are multiplicative. In this case one can write il 


that from Eq. 2 
L 


For of N at given and @ 
the random terms « and €,, VN, may be considered independent and normal with 
zero mean. It is convenient to assume that the same terms remain independent — 
in different measurements of N, although the error from lack- of-fit, €, is likely 


| 
| 
| 
q 
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Uncertainty on given D, is caused by of ‘the of 
the limiting dry densit densities Ymax and and Y in the unt poses more 


D, in terms of y,, dry unit weight, and of) 


D Vmax 
D, 


laboratories 


performed m maximum n and minimum am dry density for the 


_ cohesionless soil. Using data in this last reference one can separately infer 2a 
ro the variance of limiting density values reported by the same lab, s’, a measure - 
of within-lab variability, and th the variance of the value determined 


sity determined non-ASTM’| bop! lbs s? = 0.54? pef? 


s? 


ie) ASTM' | | lab: 


sity determined non-ASTM? ‘ lab: = 0.467 he, 

vector of parameters in reference N .06 r 60.034 0. 
— D, relationship in equation2 | 0. 00092 

lack of fit term, equation 

residual term, equation4 _ 


measurement error in the ‘ 0280 to 
measurement error in the field 0.0 


‘American Society for Testing Materials. 

different labs, a measure of interlab variability. Values of s’ and s?, are 
Ss in Table I. They are _ for a particular sand but they can be considered 


representative of all uniform medium sands. The distributions of + Tome © iven 
aNd Of Yuin, BIVEN Y min» both within and across laboratories, are nearly 


*. normal. Although s and s,, are small compared to bas average values of the 


-limitin densities, the sensitivity of D, to y and mis such that these errors 
& y max 
Cit Sia 


could not be neglected when regressing N against D, pale 

- Inaccuracy of relative density data can be accounted for in the regression — 

- of N, on D, by substituting ‘for D, in Eq. 4 in terms of yp Ymax and Yin > 


This gives the nonlinear regression model ’ 


| 
am 
fa] 
| | 
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3, are ‘ase control and a,, a3, Y min aFe 
regression parameters. An additional complication of the model is that data 
are from sands of different types for which y,,,, and y,,,, are different and 7 
the parameters a,, a,, and a, may or may not be the same. Due to these 
features, the following special procedure has been devised for parameter 


— Analysis of Regression Model in deci 6.—Let the index k(k = 


n,) identify the type of sand and suppose that for the kth sand type, M ae 
provide limiting density data. If (y>,, ..., y,,) denote the n, values of yi... 
from the rth laboratory, then before regression “analysis, Youn; may be considered 1 


to be a random va variable w with mean mean value 


= 
The above expressions are obtained by specializing results of the linear Gaussian Ss 
- model with noninformative flat prior (12). If n, = n for all laboratories, then _ 
7 and 8 take on the simple 


Similarly for m 200 the course of regression analysis, the distributions 
Of ys ‘ax and yi. are mee in the same way as the distributions of the parameters 


ze 12, and ay. - This i is achieved by including the the quantities 


_among the components of the parameter vector B, in 
Bayesian analysis | can be used for the ‘estimation of the vector B, provided 
one knows its prior distribution. This distribution can be taken to be multivariate 


normal, B ~ N(m,,%,). The inverse of the prior covariance matrix, @a) 4 


is diagonal with the firet om a, coefficients we to zero, to express 


‘et 
= estimation procedure is much simplified if the right-hand side of on 
6 is replaced with a linear function of B, e.g., with the linear Taylor t Ea. 
expansion about the peter mean values Of , and 1d min and the posterior mean 


— 

| 

| 
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ESTIMATION OF SPT- 


a of a,,m (Notice the more common of about 
the prior mean * of B does not apply here, because the prior mean of a, is 
not ¢ d.) One can om write, for the ith measurement of sand ndk, 


Is in the prior mean value of In D 


2,49 


aun 


Ynax 


4 One may ‘rewrite Eq. 14 as 


== 
AY ant wen 
and €, is a vector of independent normal errors with zero mean and unknown — 


Bayesian regression with unknown residual variance is feasible and reasonably _ 
oe only in the case when B and a7 have conjugate prior distribution (6). 
- Unfortunately, our prior information is not of this form and we have to resort 
to an approximate procedure. (Modification of standard Bayesian regression — 
would have been needed even in the case of conjugate prior, because the matrix a 
. H depends on the a priori unknown quantity mm ,)- The procedure is iterative. 
First one assumes some values for _ and With these parameters given, 
the matrix H is known and the | prior ‘normal distribution of B is in in conjugate aa 
a form. Under the same condition, the dare ataal distribution of B i is s normal, with 
q a mean valns (1) 


” 


and covariance matrix the 


— 


1 
L 


An estimator of the variance, , is 


N 
a a 
In the 
- _ _in which vector z contains all the observations In N, ,. H is a matrix with -_ 
| 
veges 
4 


7 from Eq. 20 und the estimate of a? from Eq 'n are not close to the a 
‘ initially assumed, thea these values are modified and Eqs. 20, 21, and 22 are 


value i is attained. Convergence i is usually fast. oa 
100 


‘The number neor each 
isa, in psi 


— 
ao 


= 


RBMS 


FIG. 1.—Reference Model between SPT-N and D.: (a) Laboratory Data and Fitted 
Relationship; (b) Initial and Revised D, psi = 6. 89 ) 


__ in which n, and n, are the dimensions of z and B, respectively, =. | 


_ 


— n Ref. 3, , the Reid Bedford model sand (RBMS) in | Ref. 4, and the : standard 
concrete sand (SCS) and the Plate River sand (PRS) in Ref. 5. Analysis was 
a in two stages. At first, the four sands were assumed to have different : 

coefficients a, ,a,, and a,. Because the posterior mean values of these coefficients 

were nearly the same for all sands (in consideration of their posterior standard — 
_ deviations), it was assumed, in the second stage of analysis, t that the regression © 

of Eq. 2 is identical for all sand types. Under this condition, the best ontimates — 

4 of the limiting densities changed from their prior values by only small amounts. | 
by The regression variance from the final analysis, ¢ a7, = 0.094, is not aoa 


estimation errors by pooling together information for all sands. The mean 
: vector and the covariance matrix of [a,,a,,a,]” are given in Table 1 for & 
_ expressed as a fraction and G, in psi. The covariance matrix includes the | 

of uncertainty on the limiting dry densities of the four sands. The eae z 


_ of the residual terms, e¢, in Eq. 2 and, iin in Eqs. 4 and 6 are given in ol 
Eq. 2 with parameters a, and residual distributed according to 
is called here the reference relationship between N and D, for normally consolidated — 
4 sands. This relationship applies for G, in the range of 1-50 psi (6.9-344 kPa) . 
and refers to actual values of N and D, under laboratory conditions, i.e. o 
and D, determined from a large” ‘number of laboratory measurements (D, 
is calculated from Eq. 5 with average limiting densities from a large number - 
of laboratory measurements according to the American Society for ‘Testing 
Materials procedure). Eq. 2 is shown in Fig. l(a) after setting the a, equal 
to their posterior mean values, 4, = m*, 2 = = 2. 67, 4, = m*. = 2.06, 4, = 
, = 0.442 and for « = 0. Also shown in Fig. 1 are data points used in > 
the analysis. The horizontal axis Bives values of D, calculated from Eq. 5 using 
the posterior mean values of and Fig. ‘shows the difference 
between posterior and prior mean values of D, for each of the four sands, 
as a function of the posterior mean of relative density. For example, consider 
the point with coordinates D, = 27% and N = 2G, = 9 psi) in Fig. \(a). 
This point corresponds to a measurement of fine sand. The associated difference 
between posterior and and prior mean of D D, is, from Fig. 71%. This means 
‘that a priori the mean value of D, was 27% + 7% = 34%. Updating of D, 
is a consequence of updating the estimates of the limiting densities. - 
Be. The reference relationship between N and D, cannot be used directly in the 
“field, because of factors, such as fabric, fines content, and consolidation history, 
and because of different ways in which the Standard Penetration Test is conducted 
(9). The next section shows how one c: can generate site-specific models — 


Correction of Equation 2 ror ConomionS 
_ Visual display and formal analysis of data from 17 different sites confirms — 
existence of a log-linear trend of N versus dD. . In none of the 17 ad cases was 


we — 


— 
_ The regression procedur 
at ciean San wi Stain Size The Tine sand 
| means that the final model fits the data almost as well as the four separate 


it ated to wade the slope parameter a,, but significant and sometimes very si 


_ important differences were found in the dependence of In N on depth, h. The | 
_ physical explanation is that factors such as fines content, fabric, consolidation A 
* history, etc., may vary considerably in-situ, thus distorting the dependence oe 
co E(in N) on h as predicted from Eq. 2) with o, a function of depth. For this 
reason it is mers ‘fealistic, in the field, to the parameter, with a 
ide 


4 


ve 


= 


Elevation (f 


! 
d 


af, 
5 2.—Relationships between N and D_ at Four Different ft = 305 


ric function of depth, s.(h)' er from site to site. F 
= g,(h) + a, In D, + 
he residual term €p has a structure similar to & in Eq. 4 1.€., €p ' - 
‘in w, in which €,, y,i8 measurement error with the same probability distribution — “ 
in Ea. 3, and €, represents lack-of-fit at the site. 


he 
| 
&g 
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TIMATI N OF 


_ Eq. 23 assumes that the deposit has nent stratification. If percent f fines, 
fabric, degree of consolidation, etc., vary -y also horizontally, then clearly g,(h g 
should be replaced with a more general function of location. ah il 


‘9 Estimation of Parameters in Eq. 23.—We take as our objective nthe: field 2 ay 


od to estimate parameters of In D, (D, = actual relative density in situ) as a random 


function of spatial location. Data consists of field measurements of N and D_. 

_ Because measuring D, is more “costly than measuring N, it is. typical that 
measurements of only blow-count ‘number are taken at ‘many more locations — 
than joint measurements of (N,D,) and that the latter are at discrete points — 

along one or very few vertical borings, 
- As a first step, the differences In N, — G, In D, — 4, In G, are calculated 
at the locations in which N and D, are stimultaneously measured and plotted 
versus depth, h, to estimate the vertical trend g,(h) at the site. Plots of this 
type and fitted trends at four sites are shown in Fig. 2. Data from the first 
site (these were used in Ref. 7 to conclude that “ no usable correlation exists 
between density and SPT values at the site’’) show a strong dependence of > 
— «g, On A, at least below a depth of 11.5 ft (3.5 m) . The other three cases, 

all drawn from ‘nuclear power plant foundation studies , show no significant” 7 
variation of g, with depth. Although in these cases one can take g,(h) as a 
constant, the fitted values differ significantly from the laboratory estimate, a, 

= 2.67. Analysis of data from 13 other sites gave similar results: considerable — 
Site-to-site variability of g,(h) and sometimes apices variation of g, with 
depth. As an indication of this site-to-site variability, if one considers g, (h) 


to be independent of depth at at each | site, then its variance over 17 sites equals» 


0.777 "with an associated mean value of 3.01 (2). Analyses of this type. a all 
~~ the importance of obtaining a few joint measurements of N and D, at the site, ad 
i order to estimate g,(h) and remove bias from the relationship between N © 

a... The function g, provides information about in situ fabric, overconsolidation, 
etc., as a function of depth. For example, a value of g,(h) much larger than 
the laboratory estimate a,= 2.67 indicates overconsolidation or in situ structure 


with considerable stability u under undrained cyclic loading, such as that induced y 


G,, the higher the liquefaction resistance of in situ soil compared to a sample 


in the laboratory to the same D,. 


of « oF, can n be calculated a as the average squared of In 
in D, — G,\nG, from g,(h). As examplified in Fig. 2(a), fitting a depth- -dependent - 
trend g, reduces | the variance of the lack-of-fit term €,,0; =o; - 


Spariat VaRIATION oF N AND oils ig yd 


' _ Estimation of the mean value, variance and correlation function | of In D, 
is pursued conditionally on given g,(h), a,, and a,. If this function and parameters — 
cannot be assumed known, (e.g., equal to 2,(h), @,, and 4,), then a viable — 
_ procedure is to simulate them numerically and each time estimate the parameters — 
of the random field In g & _ This simulation approach | was used in Ref. 2 to 
probabilistically ‘characterize liquefaction resistance of a soil deposit. 
_ Given g,(h), a, and a,, one can calculate, for each omnes N; of Ni Bd 


4353 
q 
| 
| 
| 
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group oup these according to’ depth. The sample mean and variance he 
of the residuals at discrete depth, h,, are estimators of the mean m4,,y,,'(h,) 
ye and the variance Cnn, ‘(hy ) of (in N r) The same quantities are estimators 
of a,E {in D,(h,)] and a? (h,) + Notice that the lack-of-fit term «, 
_ does not contribute significantly to the sample mean of the residuals because 
“= €, varies on each horizontal plane around its zero mean. Thus, one may estimate _ 
first two moments of In D |.(h,) as 
Pn p,(h,) — (h,) 


“and p,(h) = max 40, — [o2,y, th, - 03, 
In some cases it is useful to know the degree of correlation in space of 
es blowcount number and of relative density. A realistic and relatively simple 


_ model for the covariance function of (In N,)’ is one of the following type. 

Bi(ry,hy, hy B, (r 


in which r, = [(x, — + - '/? is horizontal distance, 
— Ah,| is vertical distance, and | B,, B,, B, are the covariance functions of a 
In D,, €, and Tespectively. The additive form of Eq. 27 holds under 


“covariance functions B, and Bi are further assumed t to © have similar ‘ “separable” if 


=o 
whereas is iven b 
expr ess the fact that ei y, is not correlated. The term equals the 
ance of the error in the laboratory (0.028) and = 
with a7 estimated by the procedure at the end of the last section. 
4 "AS a special case, one may set r, = 0 to obtain the vertical covariance 
hash )= [a3 Finv, (hy ) on, (h,) + er.) + I(r, ). BD 


in which Ir,) ) = 1 if r= = 0 and /(r,) = 0 otherwise. Similarly, by setting 


= one obtains the horizontal covariance function bg 


J 
‘ 
{ 
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—ESTIMATIONOF SPT-N 
llieais observations of N on a regular (ec. Su ; , Tectangular) network of points = 
it is possible to form estimators of the covariance functions in Eqs. 29 and 
; in which the summation is over all measurement pairs at depths h, and hh * 
and at the same horizontal location, and n h,r, 1S the number of eoch pairs. — 


Walk 


«FIG. 3.—First Two Moments of (in N,)’ and D, for the Site Analyzed sos Ref. 2 
Ti. 


Similarly, for the covariance function on the horizontal at depth, 


1/2 


x,)? ] = ry and depth, h. The number | of 


From knowledge of By, and it is possible to estimate the vertical 


and horizontal correlation functions, p, and Py: Convenient estimators are 


] 
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_ Summations are over the values of h, for w which ow + r,) exists 
,, terms) and over the values of h at which h) exists terms). 
Notice that the functions in Eqs. 35 and 36 do not depend on the form 
of g,(h) or on the parameters a, and a,. ‘Therefore, even if g,(h) and a,, 
a, are simulated numerically, these fun functions need to be calculated only once. a 4 
oe , exemplify the analysis, data on N has been processed from an actual am 
sand deposit. In situ measurements were made on a 6 x 6 horizontal square 
- grid of points, with regular spacing of 200 ft (61 m) [hence over a horizontal 
SS of 1,000 ft x 1,000 ft (305 m x 305 m)], and at several depths up 
te 105 ft (32 m). This data was first processed to obtain the probabilistic 


of the random field (In N )’. The mean value and the 
F 


“horizont 


ncti 


Correlatio 


. 


pk 200 300 40 00 
(ft.) Horizontal distonce, 


FIG. 4.—Correlation Function of In D, tor the Site Analyzed in Ref. 2 (a) V Vert ical 


Correlation; (6) Horizontal Correlation (if ft = 0.305 a ‘s) 


«* the field are shown in Fig. 3 as functions of depth. The same — shows * 


- the median of D . [equal to exp (1/a, m n igs )] and the variance of In D,, 

Cum The standard deviation of In D,, Fino,» equals in approximation the — 

- coefficient of variation of D,. Fig. 4(a) shows the sample horizontal correlation 

- function 6,(r,) and a fitted parametric model of the simple exponential type. 
As shown in Fig. 4(5), data on horizontal correlation is scanty . Assuming that 
the correlation function is exponential - also on horizontal planes, the model 

in Fig. 4b) can be fitted. It is worth ‘mentioning , that very similar correlation 
results were obtained at the same deposit for other soil properties, e.g., grain a 


size distribution This fact gives support to the correlation 

‘consistent of is for inference of D, from N under 


| 
a 
| 
| 


~ 
both field conditions. The model in the laboratory (reference 
a model) is based on a large data set, which consists of essentially all available — 
x laboratory information on joint measurements of N, D, and G,. It takes explicitly — 
into account errors in the laboratory determination of D, and is of a form — 
makes its use particularly simple. The reference model is not directly 


' applicable in situ, because of the effects of factors” such as in situ structure =. 
i consolidation history, fines content and variations in the way the SPT is performed 7 7 
However, the model can be applied in the field after a simple, but sometimes _ 
depth dependent, modification, which requires a few joint measurements of ithe 
NW and D, at the site. Once constructed, the site-specific model can be used Eon 
- to infer D, from N at individual locations and more importantly to determine es 
a the: mean value and the variance of D,asa function c of depth from measurements 7 q 
_of N along a large number of borings. The procedure is a numerically | 7 4 
by using data from an actual site. _ 
The ‘National ‘Science Foundation supported this study under Grant es. 
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= parameter in N-D, relationship, Eq. 2 (i = 1, 
= covariance function of random field X; dy 
¢ ), B, 2(1), B,(-) = covariance functions of a, In D,, «, and _ (Eas. 
-depth- and site- “dependent intercept in  N- D, relationship 


observation matrix with | ith row given by Eq. |! 19; — 


_h,h, = depth, with continuous or discrete values; _ 
Peli 8) ei) = function of x, equal to | for x = 0 and equal to 0 
index for limiting density measurement from rth labora- _ 


subscript indicating a quantity measured in the lab; 


= quantity defined by Eq. 15, for | ith measur ae 


as 2. number of laboratories providing for 
m, = mean value of variable x or mean value vector of 


mim, prior and posterior mean values of x or 


blowcount number of the SPT; 
of measurements; 
= number of limiting density measurements from rth — 


= number of sand types in the data; 
dimension of vector z orB, respectively; 
ny, = number of terms in summation of Eqs. 33, 34, 35, 36; 
index for laboratory providing limiting density data (r 
= horizontal and vertical distance between | two points; 


= within lab a variance of of 


| 
: 
&g 
le 
%y = horizontal coordinates; 
+= vector containing observations of In N, , 
= coefficient defined by Eq. 16 for sand k; 
uh aa = coefficient defined by Eq. 17 for observation i, on sand i 


a * = vector of parameters to to be. estimated, » Eq. 13; oP 
= maximum and minimum dry densities; 

= = jth measurement of limiting density of sand, k, from 
in =. deviations of , from their prior mean values; 
€€, 5 = lack-of-fit terms in N — D, relationships in the lab or 
= residual term in Eqs. 4 and 23; 
vector of values of €, for observations of In. N, &§ 
ae = error in measurement of In N in the field or in the . a 
measured d parameter and the outcome (Of its i ith measure- - 7 
and vertical | correlation function of Inf D, and 


covariance matrix of random vector x; 
and posterior covariance matrix of of x; comparing 


variance of random variable, 
27 


a B kitty: cal - 


= vertical effective stress; and 
“mean value or sample average as an estimator ‘ofa : random = 


vo 


z= 
Bea 
4 
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STATISTICAL ANALYSIS oF SAND LIQUEFACTION 


i Michael N. Fardis' and Daniele Veneziano,” Associate te Members, ASCE. small 


Because of the large damage - potential of ween failures, the | profession 


_ has made a major effort to understand and design against this phenomenon, 
parila under earthquake conditions. Although explanation of the physical 
‘mechanism and _ extrapolation from laboratory experiments t to field conditions — : 


still p pose se controversial problems, a ‘number of useful design yn methods are now 
available. The simpler procedures are empirical; they consist of comparing — 
- ‘Standard Penetration Test results at the site and a design earthquake acceleration 
_ with similar data from sites of known performance during historical earthquakes. 
More sophisticated ‘‘fundamental’’ approaches require calculation of the shear 
‘stresses at various depths in the soil deposit and comparison v with | laboratory ; 
test data on undisturbed or reconstituted samples. oe 
Most design methods are deterministic, but some empirical procedures have 
been extended to estimate the liquefaction probability from given, Standard — 


a: or other earthquake intensity parameter (5,34). Probabilistic versions of the — 


fundamental approach have been proposed by Donovan (6), Faccioli 
and Tang (13), and McGuire, et al. (20). 


Res - cornerstone of any probabilistic fundamental approach is a stochastic 


q 
Penetration Test, SPT, blowcounts number, N, and peak ground — : 


tn early studie: studies (18, 19 ,29), relative density was considered as the main sand 


_ has emphasized the importance of sample preparation method and soil strecture 
- so that present design practice relies on tests of ‘‘undisturbed’’ samples. Validity 
« this procedure is seriously diminished by the inability of sampling techniques 

to preserve in-place density of cohesionless material (except for medium-dense 
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sands) (3), in situ structure, and the effects of previous stress history (23) * 
Ps As a result, laboratory testing of undisturbed samples seems to offer little o 7 
= no advantage over testing of reconstituted samples, provided that in situ D, <7 
can be reliably estimated and appropriate corrections made de for the effect of 
— probabilistic model of liquefaction in the lab for undrained simple shear — : 


tests is developed first. The model gives the probability distribution of the 
A number of cycles to liquefaction as a function of uniform shear stress amplitude, | ¢ 
mean initial effective stress, and relative density. It also accounts for sample — 
preparation method, system compliance, stress nonuniformities, and errors _ 
_ inherent in the laboratory determination of relative density. Parameters of the 7 


model are estimated from a large number of test -Tesults obtained by several 


of undisturbed medium dense sands and reconstituted samples and on the effects 
4 of motion bidirectionality are used to convert the laboratory liquefaction model 
into one applicable in the field. Finally, the model is generalized to load ool 


q 

For given cyclic shear stress and initial ef fective stresses (vertical and lateral) 
“the i in situ liquefaction behavior of a sand depends on ‘grain size characteristics, > 
} telative density, soil structure, and previous static and dynamic stress history a 
(28). Because of the complexity of these effects, a liquefaction model based 
on tests of undisturbed samples would be most desirable. Unfortunately, the _ 
imitation of current sampling techniques and of available data forces one into - 
using indirect approaches, such as that proposed here. First, a probabilistic ; 
‘model i is developed for liquefaction | in a the laboratory as a function of relative — 
P density, D,; initial mean effective stress, ¢°; and sample preparation method. 
‘Then the model is transformed for in situ conditions by incorporating the effect 
of in place structure and multidirectionality of motion. 
Analysis of Simple Shear Test Results.—Undrained cyclic simple-shear tests _ 
Provide the closest laboratory approximation to in ‘situ stress conditions 
vertically propagating shear waves. A relatively large data set (192 tests) can 
be assembled using published results (1,9-12,23,26,32,35-37) and classified into 
- eight groups pr ee on testing devices. Characteristics of each group are 
given in Table 1, including sand type (all sands are uniform, medium and clean), 
initial vertical effective stress So,, and sample ‘Preparation n method. The 


in which N, = the number of uniform stress cycles that produce an excess _ 
pore pressure equal to &? (a condition that i is associated here with “liquefaction”; 
+ uniform shear stress amplitude; o>, = initial mean effective stress, ed 
also the confining stress: ‘under simple shear conditions (16); ond 
D, 


| 
Ca uc Uala O a we @ UCA | 
— Tang (13) the ratio t/ D, a, is expressed as a linear function of In N,. i, ' 


SAND ANALYSIS 


Preliminary analysis with a, and a, in Eq. 1 indexed by k 
parameters 2 are re practically the same for all groups. The zero-mean Gaussian 
for | given sample | preparation method and has the same variance for all groups. 
‘Thus, group dependence in Eq. | is limited to the intercept parameter a path a 
Standard regression techniques cannot be used to estimate the cosfficients 
(a‘,a,,a,) from observed values of (N,,1/&°,,D,) because of unreliability of 
_ laboratory determinations of D, . This problem was previously encountered by 


pounds 


Testing square = Preparation = 
device me = = method 
(2) 


Ring- -torsional with Bandaijima, Niigata .22 pluviation through 
Ring-torsional with Ottawa, American 


rigid walls § | Society for Test- 
ing Materials, 


volume Ottawa, American | dry vibrations 


dry vibrations, 
Society for Test- 
ing Materials, “+4 
shaking table |Monterey, Number} pluviation through 
Ottawa, 28. Pluviation 
Society for Test-| 
Monterey, Number pluviation 
pluviation through 
Roscoe-type i pluviation through 


the following: (1) Replacing in 
limiting dry densities, y,,,, and y,,,,; (2) calculating the variances of within- lab 
_ across-lab errors in the determination of the lnniting densities and using — 
these variances to form a prior distribution of iia? and +y,,;, for each sand; — 
(3) linearizing the model with respect to y,,,, and y,,,,; and (4) using on 
—— in an iterative mode to revise the initial (laboratory) estimates of 
max and ¥,,,,, While also estimating the parameters a‘, a,, and a,. The prior a. 
_ distribution of these last parameters is assumed to be noninformative flat. 
‘The 192 data points and the fitted regressions are shown in Fig.  1a)-1(h), 


for k= = 1, 8. For | each experimental point, two values of D, are given: aA 


authors i the ahalysis OF regres Telauonsmips between SFi-N ands 
fy D,; it can be solved here by the same technique. In short, the technique consists ae 
4 


- the estimate from laboratory determinations of the limiting densities (value in 

parenthesis) and the estimate using the updated values of y,,,, and y,,;, for iy 
the appropriate group after regression analysis. The latter estimate is, in 
approximation, the posterior mean of D,. Due to large sensitivity of D, to | 
the limiting densities, prior and posterior values are sometimes different. For 
groups that include tests at a single value of relative density, updating of — 
and Y min is not possible (due to non- -informativeness of the prior ~ eee 
of a‘) so that in this case prior and posterior mean value of D, coincide. parental vat 
~The posterior mean (4 ) and variance of a‘ are also given (Fig. 1). ‘The 


55% (41%) 


48% (33%) 


(6) Group 4 (1,25) 869 191%) 
Vor (o* ) #0.648 Dy =85% (83%) 

(55%) 


Number of cycles to liquefaction, 


FIG. 1 —Labora atory Liquaaction Data and Fitted Relonahips 


‘limiting densities. The latter variances are re especially. large when the 
densities are not updated (groups 7 and 8) or when one starts with large prior — 
_ Variances because information on y,,,, and y,,,, is from one laboratory only _ 
a (groups 1, 7 and 8) or is obtained from non-American Society for Testing and 
- Materials procedures. First and second moments of a,, a,, and € are in Table 
In the past, it has often been found convenient to assume that the shear _ 


fferences 


SAND LIQUEFACTION ANALYSIS 
stress ratio esis that causes liquefaction in a given number of cycles varies ee 
~ proportionally to D, (15). This assumption has been based on early cyclic triaxial — 7 
tests by Seed and Lee (19,29). In terms of the present model, Eq. 1,theimplication _ 
a is a, = —a,. Based on more recent simple shear testing, Seed and his co-workers — 
4, 28) have presented in a graphical form a relationship between shear stress 
Tatio and relative density (for given N) which suggests that varies 
zz with nm approx 1.35. The model proposed herein, which is based on a more — 
i comprehensive set of data and accounts for uncertainty in D,, suggests a similar 
_ relationship, but with a slightly modified value of n; n = —a, /a, ~ 1.5. @ ieay” 


10; 


20% Vor = 65-70% (56-63%) 
= 60-65% (50-56%) 
45-50% (28-36%) 
* 40-45% (20-2 8%) 


a) 


ti * Dy "52% (50%) 
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remainder of this section aims at erences 
intercepts a‘. These differences are due to many factors: (1) Preparation method; 
m7) initial mean effective stress &°,; (3) system compliance and stress nonunifor- 

; - mity; (4) sand type; and (5) testing devices. After calculating the effect of | 7 
&g — the first three factors, it will be found that the corrected intercepts differ little 
_ from one another, thus, that a general model for liquefaction in the lab needs» 
. ‘not explicitly recognize § sand type (within the present class of sands) or testing | 


results suggest that preparation method, compliance and 


wih 
| 
| F 
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nonuniformity ‘be accounted for ‘through constants ants whereas : 
dependence on &°, should be included by considering mean effective stress A 
a an additional regression variable. After doing so, the intercepts a} will be 7 


a 
ta 


rows 

‘in which A‘ A‘ are corrections for the ‘method of 


Som 


group k (relative to a reference method), system compliance, and ‘platten” 
TABLE 2. -—Probability Distribution of | Uncertain Parameters 
a4 


Uncertain 


oe 


Difference in in- 0.038 0. 
ne pluvia- by 
tion through ee 
water and plu- 
Nation 
through air' or 
ay vibration’ 
Difference in 
tercept due to 
situ struc- 
Difference in in- 
tercept due to 
motion bi- 
directionality 
Vector of 
parameters in 
for li- 
quefaction in 
the laboratory, 


9, or (in “a 
parenthesis) in 
situ, Eq. 13 


im 


* Values in parentheses refer to a} 


testing (when all the factors are zero and 
= o,.,). These reference conditions are taken to be: sample preparation - 
pluviation through water, testing at constant volume (no system compliance), — 


uniform Stress, and = on = 4. 8 psi (33 kN/m* This value of 


Fs 
1 
Symbot bution | Mean | Second moments, 
= 
#9 
ia 
4 a,,(a™*™)! 0 
} 
roughness, respectively; o,., is a reference value of initial mean effective stress, =f 
| 
a Effect of ¢°,.—There is enough experimental evidence that liquefaction behavior } 


does not t depend only « on 1/8 but rather ont and o o° 


term that depends on Available experimental suggests the 

In N, + a,ln— + a,In D, + a,1n 


‘ With only one exception (26), tests within ¢ each | group were made using a ‘single 
q value of vertical effective stress °. Therefore, simple shear test data cannot 
i be used to estimate the coefficient a,. A viable alternative is to estimate a 
from cyclic triaxial tests, emuning that this parameter is the same for <oll ; 
Some of the triaxial tests reported in the literature (4,11,19,24,26) we 
performed on samples of the same sand, prepared by the same method, with _ 
the same D,, but different o”,. From regression analysis using this data (8) 
one finds that the parameter a, “Giffers significantly from sand to sand; therefore, 
a separate estimate of a, was obtained for each of the five sands in the quoted — _ 
references. Considered as a statistical sample, this set of five estimates has 
empirical mean m a = 72 and vi variance = 077 Because i it is unclear 
which sand characteristic is responsible for this “dependence, a, will have to 
= be treated as a random variable with the moments given previously, = oo 
In a recent state-of-the-art paper (28) Seed applied a multiplicative factor 
a of 0.9 to the shear stress ratio causing liquefaction in the lab under G? = ‘ 
7 8 psi (55 kN/m* ) in order to “allow for increased pressures, up to about 1 : 
_tons/sq ft, which would be applicable | in the range of interest in the field. 


_ According to the present model and for | tsf (2,000 psf or 95. 7 kN/m? ) this 


correction factor equals [2000/(8 x 144)] aoe Te 92, in good agreement with 
Effect of Sample Preparation.—The method by which samples are prepared 


- influences their liquefaction behavior through modification of the soil structure. 
. Experimental results suggest that ‘sample preparation method affects the regres- PS 
sion in Eq. 7 only through the coefficient a |. Objectives here re are to quantify has 
: effect of sample preparation and to correct the parameters a‘ by ‘referring 


them to a single preparation method, which is taken to be siovlasion through 


: and the ith Preparation method depends only on i (and not also on sand type 


Assume ne that the difference in intercept A’ ‘between pluviation through as 
“testing: procedure, etc.). the parameter A) can be e: estimated = i by 


in which I, is - a (0,1) indicator variable with value | if preparation method 


iis used, and the intercept a = a, + a, ln D, + a, In & ° /Grep iS @ Constant, 
if preparation method is all that is being changed (in addition to +/%, ) —— 
= experiment. Refs. 12 and 24 provide the data one needs for this analysis. — 


" The resulting mean value and covariance matrix of the vector with components _ 


| 
4 
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the next section, a correction of the type a‘, be introduced to convert 
the lab model with samples prepared by pluviation through water, to a model 
applicable in situ, i.e., to ‘“‘undisturbed’’ samples. 
og Effect of System Compliance and Stress Nonuniformity. —Membrane stretching 
in the corners, expansion of the confining frame and compliance of the pore- pres- 

_ sure measurement system (these phenomena are collectively called — 
compliance’’) have an effect on the liquefaction of undrained samples tested 

by simple shear devices with relatively rigid walls (21). Such devices were used 
in test groups 1, 2, 5, 6 and 8. Quantification of the compliance effect is oe oo 
if one assumes that only the coefficient a‘ in Eq. | is affected. For analysis, - 
one can use data in Ref. 12 from experiments under constant volume (uncompliant) 
and undrained (compliant) conditions. In this set of experiments, sand type 5 
and all other testing parameters except 1/6, were kept the same; ; therefore, _ 


. 


k the values 6. ‘Analysis gives A‘, as. as a 
Dy mean -0. 99 and variance 0.5”. This means that the parameters a‘ of groups — 


, 2, 5, 6 and 8 should be decreased by 0.99 in order to remove bias from 


_ Simple shear tests of specimens covered in part by flexible membranes, such _ 
as those i in group 7, suffer from n complaince due to reduced membrane penetration | 
é into o peripheral voids. Following the procedure suggested in Ref. 21, it is found 
that in this case ~ Nii. 3,0.67 ), thus, that a‘ for test group 7 should 

be decreased by 1.3 to remove this source of bies. - Results in group 4 4 are % 

already corrected for compliance, whereas tests in group 3 are free of compliance ~ 


For group 6 the effect of stress non-uniformities due to 
roughness can be corrected by adding 0.61 to a,; the correction A rough ~ 
- N(0.61,0.87) has been found through analysis of a model of the type in ~ 
5, using data for rough and smooth plattensin Ref. 11. 0 2s” 
Model of Liquefaction in Laboratory.—Previous results allow one to correct _ 
the coefficients a’ for system compliance and pilatten roughness and to refer 
them to the same preparation method (pluviation through water) and to a given 


value o,., of initial mean effective stress. The corrected coefficients, denoted . 


Values of all the terms in the right hand side of Eq. 6 are listed in Table 
3. Uncorrected and corrected coefficients are compared in the same Table. 

_ Residual differences among the corrected values reflect variations in sand type _ 
and testing procedure, as well as possible statistical error in the determination 

of the correction terms. (See variances of the quantities for A given previously.) 


_ The variance of the corrected parameters i. (0.747) is three times smaller 


| 
| 
| 
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than the variance of the uncorrected estimates (1.317), and i is of the same order 

as the estimation variances of a‘ in Fig. 1. We shall take advantage of this 
fact to replace the eight with a single uncertain parameter a, 
“ (this is the first term in the righthand side of Eq. 2), which can be written uo a 

Since it is not clear v which of ‘the testing procedures represents best the in 
- stress conditions, and since the model applies to the generic sand, the 

variation of a, due to sand type and testing procedure, modeled by e«, 7 wil 

“be included in the final in situ model. ll 

oi objective here is to estimate m,, and the variance of «, in Eq. - The 


0.95 


74°. For estimation of n m one may use the average age of the quantities ay 


or more efficient ‘Bayesian or linear-minimum- -variance which give 
different weights to the quantities a‘, depending on their estimation variance. 


From Eq. 8, m, = = -0. 64, whereas from Bayesian estimation, which is the 
procedure used in ‘Ref. 8, one finds that t m,, is isa ‘normal: variable with mean 
a 0. 66 and variance 0.677. This last variance accounts for statistical variability 
of the estimator in Eq. 6, and, therefore, for uncertainty on the correction 
of the corrected coefficients, 2“. By adding the estimation variance of m, 
to the variance of €,, which is due to sand type and testing procedure, one 
btains a total variance of 0.747 + 0.67 0 for a,. 
_ For o,., = 4.8 psi, D, portnccedly asa fraction, sample preparation by p pluviation 
through water, and an uncompliant testing device, the final model of —- 


in a the laboratory has the form, 


terms A. It is large in comparison toa? because of the strong positive correlation j 


pik. 
=— 
a 
| 
— 


+ a, 10D, 


in which the vector of {a, a? has normal distribution with 
te 


The residual term ¢ has normal distribution with zero mean and variance 0.4 457. 
te 


_ The previous model for liquefaction in the laboratory, Eq. 9, refers to 


sample of uniform medium clean sand prepared by pluviation through water 

_ and subjected to unidirectional simple shear cycles of uniform mg 

_ order to obtain a model applicable in situ to a small volume of the same e sand, 

— with the s same D,, subjected to a known time history of horizontally propagating 
seismic S-waves, one must account for three factors: (1) The effect of difference 
in soil structure and age between in-place material and laboratory samples; 
(2) the effect of bidirectionality of the stress time history; and (3) the nonuniformity 
of stress cycle amplitudes in the field. The first two factors are consider : 
next to produce a model for liquefaction i in the field under uniform shear stres ? 
cycling. A to incorporate cycle nonuniformity be given in the 
Effect of In ‘Situ ‘Soil ‘Serecture. —On the basis of experimental evidence in 
Ref. 25, it appears reasonable to confine the effect of soil structure and aging 


_ to the parameter a, in Eq. 9. If a;" is the corrected value of a,, then the © 


difference Aj. .i. = aye", will vary from site to site depending on local soil 


structure and age. One can obviously write 

in which vee is the difference between the values of a, for undisturbed samples — 


samples remolded to the same D, by the ith procedure, and is the 
correction in the lab calculated in the | previous section. Eq. 11 is the instrument 7 
by which A,,,,, can be evaluated using information in Refs. 24, 33 on the 

liquefaction behavior of undisturbed samples and samples reconstituted by five ~ 
different procedures. Nineteen values of A,,,;,. have been calculated in this 7 


way, from 19 different natural deposits. The variance of these values is — 
to site- -to-site differences in ‘soil structure and age and equals 0. 63°. BAIS) Ont . 


‘~*~ An estimate of the mean value o of Aj, vi, in Eq. 11 can be obtained by \ various _ | 
rocedures. Perhaps the simplest estimatoris j= 
in ie the index, 4 is for preparation method and depends on the deposit, | 
- . 


4 | 
7 


efficient « estimators (Bayesian, linear-minimum- recegnine ‘that 
A‘, have different estimation variances for different i, and, therefore, should 
“receive different weighting in a formula of the type in Eq. ‘12. The Bayesian 
; _ approach was pursued in Ref. 8, with the result that the mean value of Ase 
has distribution N (2.58,0.177). By adding the estimation variance, 0.177, 
- the variance for within-group fluctuation one obtains a total variance of 0. ‘5°. 


P. Thus, the distribution of the correction term in Eq. 11 can be taken to . 


_ In Ref. 28, Seed has applied a multiplicative factor of 1.8 to the shear stress _ “4 
- “ratio which causes liquefaction of a sample prepared by pluviation Gaon 

air, in order to correct for in situ structure and long-term loading effects. According ; 
to the present model, this correction factor equals exp [— with 
| are A‘, = 2.58 + 0.43 = 3.01, for i that corresponds to pluviation 7 
through air. correction factor according to the present model 


-isexp [- GB. 01)/(-S. 5.17)] = 79, in excellent agreement with the value lue suggested 


Although VaTies across sites to the distribution above, indirect. 
information on its value at a site can be obtained from a few joint measurements 
, - of SPT-blowcounts number N and D,, on the basis of the fact that factors — 
‘such as soil- -structure, time under pressure, stress history, me . affect in the 
same direction both the liquefaction resistance (through a and the value 
(of N for given D, and (28). Ifa high positive correlation coefficient i is assumed 


measurements of (N, D,,&,) to limit uncertainty on atthe site(8), 


Maltidirectionality of Motion.—In the of seismic response it is 


is selected with frequency content and intensity (e.g., peak acceleration) equal 

to the average f: requency content and intensity along two orthogonal horizontal 

’ directions at rock outcrop. The effect of bidirectionality of motion (more 
7 _ Pronounced rearrangement of the grains | and smaller liquefaction | ewe 


shall assume ‘that ‘this ‘model is s used and ‘that the unidirectional input a We a 


Let RD be the factor by which the densification rate of dry sand as 7 
vertically propagating multidirectional S-waves exceeds the same rate under 
a unidirectional input motion of average intensity. The value of RD depends 
on details of the ground motion. Its probabilistic characteristics can be estimated 
7 from experimental results in Ref. 27 on the densification caused by artificial 
4 acceleration time-histories applied to dry samples, separately or simultaneously 
q in two orthogonal directions. The latter condition realistically simulates actual . 
4b bidirectional ground motions: when averaged over a sample of 8 simulation 5 
a pairs, the ratio of the smaller to the larger peak acceleration is 0.87, which 7 
is close to the value 0.82 estimated from actual records (22). From the same 
8 motion pairs, the sample values of the mean and variance of RD are 1.7 


ha and 0. 37, respectively. These s statistics are taken as the true ‘mean and the 
pe true variance of RD and denoted by RD and Var (RD). is he Bays 


4 
7 | 
| 4 
only the intercept a, in Eq. 9. From the same reference one can estimate the q 
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‘They 
(RD = 2) = 63 and Var [4,,(RD = 2)) = If is linear 
in RD, then A,, has unconditional mean A., = —0.63 x 1. 7/2.0 = —0.44, 
and if Var [A,,(RD)] is quadratic in RD, then Var (A,,) = 0.247. Hoecrnye 
a. suggested to apply a multiplicative factor of 0.9 to the shear stress” 
ratio which causes liquefaction i in unidirectional simple shear testing, in order 
to correct for the effect of bidirectionality of the actual motion (28). According — 
to the model proposed herein this correction factor equals exp (A,, »/a@,) = 
exp (—0.44/5.17) = 0.92, in good agreement with the value suggested by Seed. 
In summary, multidirectionality of shear waves reduces liquefaction resistance. 
This reduction can be accounted for by adding a motion-dependent random 
quantity 4,, with mean —0.44 and variance 0. 24” to * right hand side of 
Field Liquefaction Model.—For bidirectional, constant amplitude shear waves, 
‘the he number of f cycles to liquefaction in the fic in the field satisfies at 
InN, = 4 a, + a,\n D, + a, in 
in which = 4, + + Aap has N(- ~0.66 + 2.58 — 
| In this final model for field conditions and for given t, D, anda om , statistical 
uncertainty on the coefficients a, and a, and variability of a, due to 8 
type ‘contribute little to. “uncertainty o1 on N,. ‘Uncertainty on which is due 
to details of the ground motion, is also relatively “unimportant. The major 
- contributions to uncertainty on N, come from site- to-site variations of i in place 
soil structure, the effect of sand type, uncertainty on y,,,, and yi, for 
= that provided data for the laboratory model, and uncertainty on the - 
- ‘The term €, which represents the effect of minor differences i in soil structure 
for the same sample preparation method, is retained in the model for liquefaction : 
in ) the field as an uncorrelated random function of location. _ Eq. 2B is used 
n-uniform stress 
Pore unoer Nonuniform Simpce SHear Cycies b 


- model of Ref. 32 predicts that after. tiles n of uniform shear stress 
_ cycles, the undrained excess pore pressure u, is related toa® and N, as ae 


with 6 as a parameter. Information i in Refs. = » 10, 15, 26, : 35, 36 on simple 
shear testing of different sands with different valnes of D, provides 22 22 sample 


q 


SAND LIQUEFACTION ANALYSIS | 
(With D, as a. fraction, €, a log normal scatter term, and In 


order to obtain the undrained evolaticn. of u, _ under nonuniform | 
cycles, the assumption is made that, given &°, 06, and the parameters a, 
in Eq. 13, the quantity u, is a state variable of the evolution of undrained — 
- pore pressure, i.e., future values of u, depend only on its present value and 
on future shear stress amplitudes. Then, if the effect of any individual stress 
cycle on u, is small so that one can consider n as a continuous “‘time’’ variable, 
the firstorder differential equation that governs the ved evolution of u, in time i - ’ 
 Ofthetype ASCE, Vet, 100 No Peper 


in which +, is stress amplitude at at time in. The function ae should be such | 


Z that Eq. 14 is the solution of the differential equation, Eq. 16, in the oe j 


case ¢ of uniform stress cycles. Differentiation of Eq. 1 with respect to n gives: 


in which +, appears implicitly through N,. Eq. 17 can be integrated analytically 


from On to give ; 


In essence, Eq. 18 gives an expression of the ere uniform — concept 
described by Seed (28). In Ref. 8 it is shown that Eg. 18 18 corresponds to the 


quantification of this concept proposed by Annaki and Lee Ga 

3 The p present paper develops a a consistent set of stochastic models of liquefaction | 
resistance of a homogeneous mass of sand, such as a laboratory sample or 

-a small element of in situ material. A probabilistic model for liquefaction in 
the laboratory for undrained simple shear tests is developed first. This model 
gives the probability distribution of the number of cycles to liquefaction as 
s function of uniform shear stress amplitude, mean initial effective stress and o 
‘relative density and accounts for errors inherent in the laboratory determination a 
‘of relative density. Corrections for sample preparation method, system compliance 

and stress nonuniformities are calculated and incorporated into the model. 

7 Parameters are estimated from a large number of simple shear test results obesned : .: 
eight different groups of researchers using s several sands and testing techniques. 

Because of the type of sands in the data set, the model is applicable o only ae 

to uniform, medium clean sand, which is the most liquefiable class of sand. 
Information in the literature on the difference between liquefaction — 
of undisturbed and reconstituted samples and on the effects of motion bidirec-— 


= —arc sin | D> 


| 
q 
q 
tionality 1s used to convert the laboratory liquefaction model into one applicable 


in the field. Finally, the model is genccatined to load | cycles \ with nonuniform 


The major contributions to uncertainty on the number of cycles to hesdiectiiie, 
for given stress amplitude, relative density and confining stress come from 
site-to-site variations of in-place soil structure, the effect of sand type, uncertainty — 
~@ the limiting densities of the sands that provided data for the laboratory 
model, and uncertainty on the corrections of the laboratory data for preparation © 
method, system compliance, etc. The model ee here is used to predict 
the liquefaction potential of sand deposits (8). 
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The following symbols are used in this paper: | to 


= parameters of laboratory liquefaction model, Eq. 


ay intercept parameter for group k in Eq. 1; 


vie = value of a‘ after correction for preparation method, — 


msitue 


ith en) ated = - ©, 1) indicator variable equal to | for preparation method, 
=) index for sample preparation method; 
K, = coef ficient © of lateral stress at 
index for test group (k = 1,8); bee , 
m, = mean value of random variable 


og? 
= -blow- -count number of the Standard Penetration Test; 
= number of uniform shear stress cycles to liquefaction; 
haamalacneedin RD = factor by which densification rate of dry sand under 


shaking exceeds same rate under unidirec- 
Yous Vans maximum and minimum dry density; 
difference in intercept of liquefaction model for in 


soil and for samples prepared by pluviation through 


™ comp’ ™ row 
4 


it 
| 
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reference value of , equal to 4.8 psi 
= Var(x) = variance of random variable Xx; 

a i = shear stress amplitude, constant or dependent on on oe 


mean value of pa average as an estimator of arandom 


hon 


a The covers: 


-_ €, = variation of a, due to sand type and testing device; 
| 
a 
q | 
i 
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 Time-DEPENDENT BEHAVIOR OF SOLUTION 


Jan By Jamshid Ghaboussi,' M. ASCE, Randall 
creeping of salt surrounding ‘solution caverns over long of t time 
causes continuing changes in the state of stress and deformation around the 
¥ caverns. The long-term behavior of solution caverns is dependent on many 
factors. For caverns which are located at sufficient distances from the boundaries — 
of salt domes" and other caverns, the main factors affecting their long term 
behavior are cavern shape, cavern depth and in-situ stress, material pi properties 
of rock salt, and internal pressure. The emphasis in this study is on the influence 


7 i of the cavern shape and the internal pressure history on the 8 behavior — 


8) the influence of creep deformations in a cavern 1 tightness pressure test: a 
and (4) the effect on analysis results of the / way in which the leaching oe, «i 


is simulated in the analysis—instantaneous versus gradual simulation of the 


reep are discussed briefly, followed by the description of the caverns 

nd the method « of analysis. In the remainder of | the paper, various. aspects a 

of results of the analyses are presented and ee eee & 


a 
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one set of creep material properties was used throughout this study. _ doidw oi 1 
7 a number of creep finite element analyses were performed on three caverns. 
_ The cavern geometries were selected so as to cover a wide range of height-to- _ 
diameter ratios, varying from a disc shaped cavern to a long cylindrical cavern. 
_ These caverns were analyzed under several different internal pressures. In _ 
_ addition, analyses were performed to study the following: (1) The effect on § 
7 cavern behavior of a sudden drop in internal pressure; (2) the response of _ : 
verene thet af cavern in an array of enared caverns: 
& q 
| 
— 


: Material Model. —The strains are assumed to consist of elastic s 


=€ 


4 elastic material parameters: Shear modulus, G, and the bulk modulus, B: 


which 8, the Kronecker delta; ‘summation convention over the 


_ The volumetric component t of the creep strains are assumed to be zero. Thus, a 


volumetrically the material behaves elastically: 

=0 


~The creep strain rates are assumed to be proportional to deviatoric stresses. ; 


in the the time rate; and é° = the rate of the 4 
creep strain which is given by the following power law: A 
which k, m, and n = the creep constants; = the second 


invariant of deviatoric stress 

—5,5, 


y _ The power law is know n to reasonably well approximate the creep: of rock 

: salt, based on field measurements of pillar shorting (2), laboratory tests of 
model pillars (7), and laboratory testing of rock samples (1,4). Power law has © 

also been used for analytical investigation of openings in salt (1,5). It should g 
be pointed out that most of laboratory and field measurements are for time | 
ranges of short to moderate term. Reliable long-term creep data for rock salt : 


_ behavior is lacking at the | present In the absence of long-term creep data, 

_ the power law is used. However, the long-term portions of the analysis results 

_ presented i in this paper must be considered qualitatively and must only be used 

& comparative purposes. Improved material models for rock salt are ; a 

a Only one set of creep | ) material parameters was one throughout this study. % 

‘The material parameters are k = 1.36 x 10-", m= 1.45, n= 0.50. These 

‘parameters were chosen from the creep test results ‘penaalil a in Ref. 1. The 
_ creep strain from Eq. 5 with the above parameters is compared with typical | 


Bes extension test result from Ref. 1 in Fig. 1. Included in this figure is 


The material is ‘assumed to be 4 ; 

} 

| 
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0. 25% elastic strain which was estimated from the test results. It can be seen a 
from Fig. _ 1 that the power law reasonably well approximates the test results, -a 
= at a very early time where the creep Strains ar are underestimated. 5 stele ef 
Finite Element Solution Scheme.— —The creep finite element analysis is carried 


a in a time marching scheme and the solution is obtained at discrete times, 
,i=1,..., with time steps At; alia? 
rab 
The stress vector, o,, at time /,, is written as the sum of the known wd 7 
G,_,, at time t,_, plus the incremental stress vector, 
_ The incremental stress vector, Ao,, is related to the increment of the total 


strain vector, be, and the increment of the creep strain vector, through 


(oy "2500 psi bs 
| 


Elastic Strain» 0.0025 


FIG. 1 Vereus of Test and Creep Used 
in Study (1 psi=6.89kPa) wre 
the elastic material property matrix, 
Ao, = C(Ae, - 


The creep strain increment Ae‘ at time /, is computed by the foliowing equation 


= nk 


which S = the vector of deviatoric and = = 
finite- element equation can be written as follows: 
in which - = element Boolean incidence matrix; B = element strain-displacem ent - 
matrix: = BAu; Au = incremental vector of nodal displacements; P,= J 
nodal hey vector at time, ¢,; kK = superscript denoting the element number; ia 
_ and T = superscript denoting the transpose of a matrix. In Eq. 11 the integration __ 


is over the volume of each element and the summation is over all the elements © 


| 


1382, 


‘K4u,=P,—R, . (12) 
a in which K = the elastic’ alttasss matrix; and R, = an equivalent load vector _ 


4 The time response of the ‘system is acaasinie wad successive solution of Eq. 


‘The three cavern geometries considered in this investigation are shown in | 


Fig. 2. Cavern dimensions and depths are summarized in Table 1. All three . 
caverns are approximately at the same midheight depth. The cavern dimensions © 
are chosen so as to allow consideration of a wide range of possible height 
to ‘diameter ratios, h/D. Although these caverns are hypothetical, the cavern 


Depth Below 
Ground Surface 


FIG. —Cavern Geometries Analyzed (1 ft = 0. 305 | m; 1 = 02 kg/m 


Z and dimensions are similar to many in sak domes. 
_ Axisymmetric finite element meshes for each of the three caverns are shown 
in Figs. 3, 4, and 5. The meshes given in Figs. 3, 4, and 5 are for single 
caverns; i.e., it was assumed that there were no other caverns nearby that < 
* would influence the behavior of the cavern being analyzed. Cavern 2, however, 
was also analyzed f for the multiple cavern case. As shown i in Fig. 4, a 
array - of identical caverns was assumed. For such an array, the vertical planes 
a of symmetry form, in plane view, a hexagon around each cavern. Thus, due ah 
te. to symmetry, the multiple cavern case can be analyzed by considering just me. 
<3 one of the caverns and that portion of the surrounding medium that is located | 
a within the hexagon. An axisymmetric finite element mesh was generated for 
the cavern 2 multiple | cavern case, with a radius corresponding to that of the | 
circle in Fig. 4 (area of circle = area of hexagon). This type of analysis is ae 
reasonable approximation to the much more expensive three-dimensional 
1 the Caverns Studied 


Thickness of salt cover, in feet 7 90 


Diameter, in feet (ms 


Height, h, in in feet 
Initial volume, in 


(cubic meters) (1,980, ag a: mi 920,000) 
> 


(1,920, 
analysis that would be to analyze | the 
_ It should be noted that this approach to the analysis of the multiple cavern 
se yields results that are applicable only to the caverns located in the res a 
ee beats areata of ‘identically shaped caverns that are all leached at the same 
Bach anelysis can be divided into three phases as the basis of the various — ; 
solution steps. These phases are the in-situ phase, the leaching phase, and the 
operation phase. The in-situ phase corresponds to solution step | in which 
* the in-situ stress state in the salt mass is simulated. The left boundary, mesh a 
.s center line, and the bottom boundary, are assigned displacement boundary 
Bion (no displacements normal to boundary). The top boundary is assigned 
a force boundary condition corresponding to the weight of the overburden and 
cap rock. The right boundary is assigned a force boundary condition corresponding 
to vertical | of in-situ stresses. The vertical stresses 


| 
= 


resulting from the self-weight of the salt a are introduced forces 
_ - The leaching phase consists of one or more solution steps during which the 
a. leaching out of the cavern is simulated. The majority of analyses performed _ 
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Ht = 


cavern One 
assumed the cavern to be formed instantaneously at time, ¢ oe Thus, in = 
~ solution step 2 the force boundary condition at the right boundary is replaced _ ke 
= by a displacement boundary condition, the salt elements within the perimeter — 
- the cavern are deactivated to simulate creation of the cavern, and forces 
are applied to the roof, floor, and wall of the newly created cavern to Fam onal 


pressure exerted by the fluid within the cavern. _In one 


4 
ss 


of. cavern the gra leaching out of the over a of time was: 


600ft 


_ The operation phase cor consists of all solution on steps after completion 1 of the - 
leaching Se Because the mesh geometry remains unchanged aan this nil 


| 
+. 
a 


| 
: abtilia that after a cade length of time the brine is replaced by oil, or = 
_ that a pressure test is performed. These changes were simulated by altering 
the e force distributions } applied to the roof, f floor, and wall of the simulated — 
cavern, Basically, three force distributions were considered: Brine-filled 
cavern, column of brine to ground surface; (2) oil-filled cavern, column of F 
oil to ground surface; and (3) brine-filled cavern, column of brine to ground 2 pa 


: surface plus well- head pressure. It was always assumed that the internal pressure = us 


| Simulation| 
pressure leaching si 
P one-step 
| multistage 
one-step | 
one-step | no 
one-step 


column of brine, y = 75 pef a a 
= column of oil, y = = 50 "ee (800 ae » 


thet internal pressure did not increase. All the analyses are 


The cause creep strain in salt is the octahedral as can be 
seen from Eq. 5. In the in-situ phase, the octahedral shear stresses are negligible 
since in all the analyses an in-situ coefficient of lateral stress of K, = 1.02 
was used. When the cavern is created, high octahedral shear stresses occur 
around the cavern. . This -Tesults in high creep. “strain rates. As time elapses, 
_ the octahedral shear stresses around the caverns generally decrease, resulting — 
if decreasing creep strain rates. If this process is allowed to continue for a 
ee long period of time, the octahedral shear stresses will vanish in the limit, | . nie 
the of-stress around the cavern will a hydrostatic stress 
of the individual stress compon nents, the one acting normal to the cavern 


-% 10° psi (27.6 x < 10” kPa), v = 0.25. All cases: k = 1.36 x 0% m = 


AY inward creep of the surrounding salt, fluid was bled off at the surface sich 
TABLE 2.—Summary of Finite-Element Analyses Performed _ 
[Sina 
= @€=—ss Case Cavern multip 
one single 
| two single 
4 


‘surfece: must remain constant and equal to internal fluid pressure. ‘The magnitude 
* of the other two normal stresses generally approach the value the internal pressure, al 


and the shear stresses approach zero. This pattern of stress changes is the 
_ dominant aspect of the creep behavior of the medium around the cavern. Any 
-method which does not take into account such creep-induced stress an 
can lead to grossly misleading results. 

_ The above arguments strictly apply to single caverns within very large masses 
of salt rock. During the process of creep, arching takes place, and the — 


Cavern One 1 

— Insitu 4 


- 


FIG. 6.—Radial Distribution of Stresses in Roof and Floor of Cavern 1 at ¢ = a 


- will continue if the extending zone of creeping medium does not intersect cuniler 7 
zones from other caverns. Thus, it can be > expected that the behavior of a , 
cavern within an array of caverns will differ somewhat from that of a single - 
cavern. In a multiple cavern situation, arching is not possible and the average . 
- vertical stress in the salt between the caverns must remain constant and equal 
to the overburden stress. Therefore, the limiting value of the octahedral shear 
= is determined by ‘the: difference the average and 


| 
| 
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a To illustrate the main aspects of the c creep-induced stress changes, the stresses » 
around the three caverns plus multiple case of cavern 2, under constant fluid — 
pressure due to column of brine, are explained (cases 1, 6, 8, and 9 of Table 2). 

The distribution of stresses in the roof and floor of cavern one are shown © 
in Fig. 6. There is little variation of stress magnitude with radial distance from 

the center line except around the corners. It can be seen that, upon instantaneous 7 
creation of the brine-filled cavern, the vertical stress, o,, immediately drops _ 
to a value equal to the internal pressure. In response to this stress reduction, - # 
the radial and circumferential stresses also drop to values less than their in-situ _ 
"magnitudes. There is a general stress reduction with the passage of time. However, ey 
immediately above the roof and below the floor, the vertical stress remains a 


unchanged, being controlled by the constant value of the internal pressure. 


> 
Line 


— 


Stresses in Wall of Cavern 2, Single - Stresses in Wall of Cavern 2, Multiple * 
Cavern Case (Case 6, Table 2) (1 psi = Cavern Case (Case 8, Table 2) (1 psi ms : 


With time, the magnitudes of the radial and circumferential stresses tend to 
approach that of the internal pressure, and VJ, approaches zero. 

The vertical distributions of stresses in the salt just beyond the 
are given in Figs. 7 and 8, for cavern 2 and the cavern 2-multiple cavern case, 
respectively. Upon the instantaneous creation of these brine-filled caverns, the 
radial stress drops from its in-situ magnitude to that of the internal pressure. 
_ For a single cavern, the vertical stresses near the wall essentially remain at 

4 their in-situ magnitudes because the “effect of arching is concentrated away 

_ from the cavern wall and over a large area. However, for the multiple cavern 
case, the vertical stress initially carried by the salt removed to form the covers 

7 must now be carried by the salt remaining between the caverns, thus resulting J 
in an increase of vertical stress magnitude from the in-situ conditions . With 


14 
(| | VW 
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the passage of time, the radial and stresses tend to appecach 


zz magnitude of the internal pressure. The vertical stress in the salt adjacent 


gon 


ae 


16,380 hrs = 187 


FIG. 9.—Contours of Square Root of Second Invariant of Deviatoric Stress" Around 


1 1, Table 2) (1 psi - = 6. 89 kPa) ob ib 


FIG. 10. —Contours of f Square Root of Second sci of Deviatoric Stress for —— 
2 (Cases 6 and 8, Table 2) (1 psi = 6. 89 we) met = un 

toa single cavern also tends to internal pressure. Additional 
of vertical stress to the salt farther away from the cavern Are pay i 


390 


stress at t any pestioutoe point in the ‘“‘pillar’’ between caverns may vary vith 
Contours of Vd; values in the salt around the caverns are given in Fig. 
9 for cavern | and Fig. 10 for cavern 2 (single and r multiple cavern cases). 
_ The second invariant of the deviatoric stress tensor, J,, is a convienent measure 
of the shear stress, and is related to the octahedral shear stress as follows: 
4 Fig. 9 shows that at ¢ = 0 the very high VJ; ends are confined to relatively _ 
$ small zones around the “‘corners’’ of the cavern. Along | the center line above 
f —s and _- below the cavern, the largest vi values occur some distance away from 
the cavern, instead of right at the roof and floor as is the case for cavern ~ 
2 (Fig. 10). As shown in Fig. 9, after_a period of approximately 2 yr — 
has been a significant reduction of VI, , magnitudes around the cavern. The 
highest stress differences remain around the cavern wall, the V J, values there 
being almost twice the magnitude of those in the roof and _— at the e center a 


Fig. 10 shows that high shear stresses for cavern 2 occur at it the cavern walls, 


“and the zone of high shear stresses extend some distance beyond the walls. a 
As shown in Fig. 10, the Vd, contours for the multiple cavern case do not _ 
- differ at all that much from those for the single cavern, _ although it is obvious 


_ that the differences do exist. Stresses around the caverns for the multiple cavern” : 


In general, creep of the surrounding salt results in the dissipation of the 
“high J, concentrations in the immediate vicinity of a cavern, as can be seen ; 
in Figs. 9 and 10. However, concurrently there is an expansion of the zone 
of resulting i in a greater volume of salt that is subjected t to 


7 _ case, of course, are quite dependent on the spacing between caverns. ‘=e 


‘ 


examination of the volume change, vc, versus time curves obtained 
from these calculations, it was observed that if these curves were divided into 

_ increments, and if the volume change rates were calculated for each increment, 

_ the resulting volume change rate (VCR) versus time (¢) relationship, when plotted 
ona log VCR-log (t) graph, closely approximated a straight line. The log (VCR)-log 


Cavern VoLume 


; mG) plots obtained for several cases are shown in Figs. 11 and 12. Similar relations - 
_ were obtained for other cases. Therefore, it is reasonable to assume that the | 
- equation for time rate of volume change of caverns, VCR, will have the following — 


in which (VCR), = volume change rate at ¢ = 1 hi; ie slope of the log 
(VCR) versus log (¢) line; and ¢ = time, in hours. The integration of of Eq. 14 : 
with respect to time yields an equation for the f percentage of | volume « change 4 


ncase. For this 

cas etween cavers 

; es remain constant with time and equal to the weight of all material above the al - 

| 

and 
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Analyses Of Aa 


Cavern Two 


Internol Pressure Due To 
* —Cavern Volume Change Rate Versus. Time 


+ 845 psi 


ll FIG. 13, —Cavern Volume Change Versus 


| | 
{ 
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q Fig. 13 is ’ ye log plot of volume change versus fins In this figure, the - 
results from the different finite element _ analyses are represented by the data i in 
points. The curves shown in Fig. 13 were obtained from Eq. 15. It can be 
seen that in each case the finite- —, analysis results and the equation are ai 
For each case considered in mss 13, the output from the corresponding finit 
element analysis was used to obtain the values of (VCR),, S, and VC, fo 
§ in Eq. 15. Values of of the parameters (VCR),, S, and | (VC), are given in ie 
- Table 3 for the six cases. These parameters are functions of cavern size and 
shape, cavern depth, magnitude of the in-situ stresses, magnitude of the cavern 


*S ‘TABLE 3.—Values of Parameters for Eq. 15 as Determined from Finite- Element Results 


Volume © 


percentage 
per hour 
0.02720 
| 0.05830 
0.02285 
0.00565 | 0. 0. 05910 
So . 0. 10320 
“Case numbers correspond with those in Table | 


- P, = pressure due to a column of brine. sma a a 


= pressure due toacolumnofoil. 
~ Py = wellhead pressure = 845 psi (5,822 kPa). 
7 7 internal pressure, the creep properties of the salt, and the presence of surrounding. 
_ _ cavers. On the basis of the relatively small number of analyses performed 
in this study, it is not possible to determine with any degree of certainty the 
3 relationships between each of these three parameters and the individual factors ed 
q the above list. ‘However, , given the data that is available, it does seem oa 
few comments are in order. 
“g It appears that (VCR), is a function of the difference between a magsiodes 
of the internal fluid pressure and the stresses in the surrounding salt. It can 
be seen in Table 3 that the highest value of (VCR), occurs when the internal 
pressure is lowest and the smallest value occurs when the internal pressure _ ‘ 
is highest. Also, comparison of values for single caverns 2 and 3 show that 7 
the higher value of (VCR), occurs for cavern 2, which extends to the greater 
~ depth and therefore is subjected to the larger o,,, — P, stress difference. 
Comparison of the for the single cavern and cavern 


4 


cases of cavern 2 indicates that the presence « of adjacent caverns wei, 
the magnitude of the parameter. Note that there is a considerable difference _ 
s between the (VCR), values of the brine-filled caverns | and 3, which are both — J 
located at the same midheight depth. It would seem that differences of — 
shape | are primarily responsible for the difference between the (VCR), values 
of these two caverns. For cavern 1, the majority of volume change is due _ 
to inward displacement of roof and floor, whereas, for the other caverns, wall y 
_ displacements account for most of the volume change. = | bese 
_ The S values given in Table 3 seem to indicate that the parameter, Sa 
a insensitive to most of the factors that influence (VCR) sik It would & 
4 the presence © or absence of adjacent caverns. The value of S may be linked 
to the redistribution of stresses around the cavern during creep. There are no . 
7 obstacles to stress redistribution for single caverns, and, therefore, the creep 
‘Tate (and thus the volume ' change rate) is allowed to reduce Telatively rapidly. 
TABLE 4. —Velume | Changes (Closures), asa Percentage = 4. 


bos 
Cavern multiple 
single 
‘Single 
single 
‘multiple 
single | 
“Case numbers correspond to those in Table 2. on 
P, = pressure due to column of brine. 
= pressure due to column of oil. she ing 
= wellhead pressure = 845 psi (5,822 kPa) > 


& wen, for the multiple cavern case, the amount of salt surrounding any 

one cavern is rather limited. Thus, stress redistribution is hindered, and large _ 
Ste te differences remain in the salt for a longer time. This weute cause the 
the 


_ The percentages of volume change (volume reduction or closure), as obtained — 
_ from Eq. 15 for the six cases considered, are given in Table 4 for five aiterat 
times from 7.5 yr-40 yr. As shown in Table 4, the smallest volume change 
occurs for the single cavern 2 with an added wellhead pressure. These values _ 


this same cavern is surrounded by an array of identical caverns. At 7.5 yr, 
(A vi V)% for the ee cavern case is about 140% larger than that for = 


| 
i 
‘i - are less than 25% of those for the same cavern without this additional internal | ; 
q 
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the difference is 200%. Comparison of cases 9 and 6 dente that (A V/V)% = i 
for cavern 3 is approximately 50% of that for cavern 2 which h has almost three 
times the height, of cavern 3 and therefore extends to. much greater 

Cavern Response to Suopen Decrease oF INTERNAL PressuRE 
A sudden decrease in the internal pressure of a f a cavern would « cause a : 

neous elastic « changes in the state of stress around the cavern. Generally, the 
octahedral shear stresses are suddenly increased and as a result the rate of oe 

creep increases. Presented in this section are the results of analysis of simulation _ : 


of pressure decrease in cavern | (case 3 of Table 2). In this analysis the cavern 
was subjected to the following history of internal pressure. 
1. t = 0: Leaching of the cavern (one-step ‘simulation-a -assumes leaching . 
_ 2.05 ¢ = 7.5 yr: The surrounding salt is allowed to creep, the cavern 
‘internal pressure is maintained constant, equal to the cavern filled with brine 
a column of being to ground surface. 
t= : Instan e from 
a in 2 to that for the | cavern 1 filled with oil, plus a - column of oil beds une 
- 4. 7.5 yr = t = 15 yr: The surrounding salt is allowed to creep in = 


After, as before, the pressure change all stresses are compressive. _ 
” In the salt comprising the central portions of the cavern roof al leer i i. e., 
a away from the ‘‘corners,”’ the vertical stresses, @,, undergo a decrease of 
magnitude that is identical to the internal al pressure change. In these same locations — A 
the radial and circumferential stresses, and o,, also decrease i in magnitude, 
but not to the extent that o, does. Although there is a general stress reduction 
7 in these locations in response to the pressure drop, because o, decreases more - 
than do o, and o,, there is a resultant increase of the stress differences. Thus, 
the magnitude of VJ, immediately after the pressure drop is greater | than that 
_ just before the pressure drop. The magnitudes of VJ, in the cavern n wall also 
_ increase immediately following the pressure drop. Here, however, this is to 
a great extent due to the large increase of the shear stress, t,,, rather than L 
a the changes of the normal stresses. In general, it can be said that in response e 
to a reduction of cavern internal pressure the magnitude of the octahedral shear > 
_ stress, Vid, is increased everywhere in the salt surrounding the cavern. As — 
; a consequence, the rate of creep is also everywhere increased. Saar a 
: 4 a Results from this analysis are shown, in terms of cavern volume change 
Closure), by the lower curve in the log (AV/V) versus log (t) plot of Fig. ~ 
‘See 0 <1 = 7.5 yr, the volume change values are identical to those obtained 
from the analysis of cavern one with brine pressure which was mentioned earlier. _ 
7 Extension of the solid-line curve beyond ¢ } anes 5 yr shows the cavern ie 
changes that would have resulted if the cavern aie had — 


4 
q 
| 
| | 
i : Surrounding Sa t change So as to be in equiliDrium with the new interna ressure : 
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_ maintained at that on ite The upper curve in Fig. 14 represents the cave ~— 
pressure corresponding to that for oil throughout the entire time span illustrated - 
‘(case 2, Table 2). It can be seen in Fig. 14 that, simultaneously with the pressure a 
oe drop, the salt around the cavern experiences a certain amount of elastic inward 


= changes that would result if cavern 1 were subjected to an internal 


— 


‘displacement which results in an elastic component of volume change of 
approximately AV / V = 0.35%. Following the pressure drop, the increased stress _ Pf 
differences in the salt result in a higher rate of creep and, thus, a higher rate 
of cavern | volume work sgial 


oat | 
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At 1=75 Years interno! 

Pressure is Reduced From 
Brine To Oi! 

: 


Cavern Volume Change (VC), AV/V_ 


The qualitative observations can be made the cavern 
behavior following a sudden internal pressure drop. __ 


l. rt gr the ‘internal pressure decrease, a sudden it increase in the creep 


: - approaches that of a cavern which has been under the lower internal pressure 
‘ . from the beginning. In effect, the “‘memory”’ of earlier higher internal pressure 
> m gradually disappears. Tis behavior is evident from cavern volume changes — 
An important practical implication | of these observations is that in the analysis — 
of cavern behavior, the recent past history plays an important role and must 
be simulated. However, the distant past history need not be considered since, | 
for all practical purposes, their influence has disappeared. 
 Mutristace Versus Sincue-Step Simutation oF Leacuinc 
The actual leaching out of caverns takes a relatively long period of —_ 
before a final configuration is reached. During this period of time the salt creeps, — 
‘Tesulting: ‘in redistribution ‘Stresses: Ina 


4 


‘single- step excavation analysis, the creation of the cavern is simulated in one — 
step, representing in: instantaneous formation of the cavern. Thus, in this: type 


is final configuration. In order to allow creep during the leaching period, it 
is necessary to simulate the leaching process. In this method of analysis, the a 
leaching of the cavern is simulated in a number of stages with real time intervals od 


leaching, large octahedral shear : stresses ca can be e expected to occur around the 
_ cavern, which may suggest yielding and possible failure of the rock salt in 

these zones of high octahedral shear stress. However, in reality, as the leaching 

“out of the cavern proceeds, high octahedral shear stresses are dissipated due - 
the the creep of th the s salt. 

Cavern 
Remains Unchanged 


Within These Time 
Stages Of Leaching 


Intervals 


First 


Fourth Stage Final Configuration 
Att =24Months Of Covern 


FIG. 15.—Detail of Finite Element Mesh Around Cavern 1 at Various 


to the ‘simulation of the leaching process in cavern 1 was performed in| ‘ 
four stages as shown in Fig. 15. The fourth stage leaching, which completed © 

E the formation of the final configuration of the cavern, , occurred two years after 
the leaching of the first stage. Between each of the leaching stages, 8 months ry: 


- of elapsed time was allowed. The results of this analysis are compared with 

7 the results of the single-step analysis in Figs. 16, 17, and 18. The vertical — 

_ displacements of the center of the cavern roof as shown in Fig. 16. The second 
invariant of the deviatoric stress at the same point versus time is shown in 
Fig. 17 for the two types of analyses. The sudden jumps in the displacements 
and stresses at the beginning of each stage of excavation are due to elastic 

_ instantaneous response. In a - leaching out of the cavern is gradual 
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- and continuous, _ to very small steps of excavation. In such a case, at 
each instant of time the deformations will consist of an elastic and a — 
component, and thus the portions of the curves in Figs. 16 and 17 which consist — 
of elastic jumps followed by creep should be re] replaced by continuous curves. 
However, during the leaching-out period, the rate of deformation will be higher 


the rate resulting from the single. “step This should 


Vertical Displacement 8 (inches) 


—Root Center- Line Displacement of Cavern 1 One- and 


Point Above 
Roof Centerline 


Dever 1 for One-Step and Muttiotage Simulation of Leaching (Cases 1 and 4, 
Table 2) (1 psi = 6.89 kPa) 


‘4s be expected, : since in a single- step analysis all the elastic deformation occur 
at time = 0 when the formation of the cavern is simulated. 
After the leaching phase, the displacements and stresses from the multistage 
Lou analysis approach those of the single-step analysis. This indicates that for the — ; 
_ study of the long-term >10 yr effects the he single-step ex excavation — is will 
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Be 
FIG. 16 
Rs ‘Simuta 
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quare Root of Invariant of Deviatoric Stress Around Cavern 1 during Multistage Simulation of ib 


's of Squ 
si = 6.89 kPa 


Contours ol 
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18, respectively. For the : single ‘step analysis, VI, for times 1 = 
_ (instantaneous formation of the cavern) and V J, for t = 1.87 yr are shown 
in Fig. 9. For the multistage leaching analysis, the contours of VJ, for t = 
0, after the first stage of leaching, before and after each of the leaching stages _ 
2, 3, and 4, and at time ¢ = 5. 7 yr, are shown in Fig. Bw ; ar 4 
al It can be seen from Fig. ds that in single- “step analysis there are very large 
; — zones of high VI, around the cavern at t = 0. By comparison, in multistage 
7: analysis, at the end of the fourth stage of excavation the values of V J, around» 
the cavern are much smaller. This seems to indicate that in real caverns which 
p are leached out over a period of time, very high values es of V Jz, and therefore 
and failure are unlikely to occur. he — 
_ Another observation which can be. made from Figs. 9 and 18 is that the 
| ‘contours of VI, which are initially at some distance from the cavern, 1, generally 
tend to move away from the cavern with time and their shape becomes more A 
spherical. This indicates that, with elapsed time, the effective zone of creep 
Ss the cavern becomes larger. Also, it can be noted that after a sufficient 


ont of time and at sufficient distances from the cavern, the medium behaves - a 


as though the shape o of the cavern is spherical. bluow 


Simu ATION OF TiGHTNESS Pressure Test 


ae a cavern is sufficiently free of fractures that the cavern can ‘ee 
_ used for the purpose of storage. The tightness pressure test is conducted by | 
the wellhead pressure by pumping in brine and determining whether 

the wellhead pressure can be maintained for a period of time, usually a few 
days. Usually, the wellhead pressure is increased e 90% of the in- situ stress 
- at the depth above which the pipe is grouted. nd a 


During the tightness pressure test, the vale. of the cavern is also changing 

- due to creep by the increase in the internal pressure. The role of the creep 

- volume change in a pressure test has not been thoroughly studied. The amount 

_ of brine which is pumped into a cavern during a pressure test not only represents 
“the brine which may be lost through fractures or other causes, but “ahead also 


the volume changes i in a cavern due | to creep deformations. 


_ Theincreased internal pressure caused an instantaneous elastic volume increase 
a of the cavern of about 55,000 bbl (8,800 m’). This was followed by a gradual 
. outward displacement of cavern roof, floor, and walls. Due to creep of the 


salt, the cavern roof and floor continue to move outward for about 10 days, — 


GTIO 
__-The contours of the square root OF the second invariant Of deviatoric stresses | 
q 
= ——_— was simulated in cavern 1. In this analysis, first the leaching out of the caverns 
= 6—Ssé«ta: Simulated in four stages over 2 yr, as described in the previous section, a 
i and then the cavern was allowed to creep for five and one half — ; 
with an internal pressure resulting from a column of brine. Att = 7.5 J 
— 7 yr, the tightness pressure test was simulated by instantaneously increasing the Jf 
&§ 
| 


an and then they re reverse direction and start moving inward, as would under 


wellhead pressure increase = 713 psi(4,913 kPa) 


volume decrease of the brine i in the cavern = 3,900 bbl (624 m *) ae 
¢ elastic volume i increase = 55, 000 bbl 600 m 
maximum creep volume increase = 960 bbl (154 m 


creep volume increase in | h = 800 bbl (128 m = 


a In reality, to increase the wellhead pressure by 13, psi 4, 913 kPa), ‘sufficient - 
rine must be pumped in to overcome a 3,900 bbl (624 m ny volume decrease se 
of the brine and 55,000 bbl (8,800 m’) instantaneous elastic volume increase 
of the cavern. _ Obviously, it will take more than 1 h to pump in that much > 
brine. Thus, at least the first hour creep volume change must be included in 2 
the volume of brine segues to raise the wellhead —— thy 4 713 psi (4, 913 


160 bbl (26 m*) 


ae lastic Volume Increase = 55, 000 bb! 


(instantaneous) 
tone Volume Decrease Of Brine = = 3,900 bbl 


FIG. 19. —Volume Change of Cavern 1 Pressure Test (1 barrel = 0. 16m 


ad 
(19. It can seen that Maximum increase iS { 
“However 800 _bbi (126 of this volume increase occurs within | h. The 
ac 


‘Strain of 160/1: 12,400, 000 = 13x 10° in the brine. Assuming the bulk modulus 
4 of brine to be 4 x 10° psi (2.8 x 10° kPa), the | maximum pressure change 
s during a 10 day y pressure test due to « creep would be 1.3 x 10° x 400,000 
=5. 2 psi (35.8 kPa), or less than 1% of the test pressure. This analysis seems r= 
. to indicate that creep deformations have no significant influence on the results — 
of the tightness pressure test for the cavern analyzed. 
ee _ However, it must be emphasized that the above conclusion is only valid 
for the cavern which was analyzed. In general, the following factors will have = 


a significant i influence on the results of the simulation analyses of the tightness — 


Cavern shape; the height- -to-diameter ratio. 


= 3 Recent past pressure history; whether the | cavern has been subjected to 
significant pressure changes i in the recent past. 

B order to arrive at a more general conclusion hajeoioa the infleence of the 


“creep deformations on the tightness peeurs test results, these factors must 


ve deformations was of wack 
‘The fin financial su support of the Solution a Mining Research Institute, Inc., Floss- 


AppeNDIX.—REFERENCES _ 


> 3 . Boresi, A. P. and De Deere, D D. i” “Creep Closure of a | Spherical Cavity in an Infinite 
_Medium,”’ Report to Homes and Narver, Inc., May, 1963. | 
2. Bradshaw, R. L., Boegly, W. J., and Empson, F. M., “Correlation of ene 
fe Measurements in Salt Mines with Laboratory Creep-Test Data,’ ered 
__ Sixth Symposium on Rock Mechanics, Rolla, Mo., Aug., 1964. —— 
3 Hermann, W., Wawersik, OW. R., and Lauson, H. S., “Creep Curves ‘aa Fitting 
- Parameters: for Southeastern New Mexico Bedded Salt,”” Report SA ND-80-0087, Sandia 
Laboratories, Albuquerque, N.M., Mar., 1980. aed 
. Lomenick, T. F., and Bradshaw, R. | “Deformation of Rock Salt in Openings Mined © 
_ for the Disposal of Radioactive Wastes,’’ Rock Mechanics, Vol. 1, No. 3, 1969, PP. 
‘ 5. Maxwell, D. E., , Wahi, K. K., Reaugh, J. E., and Hofmann, R., “A Time- Dependent — 
_ Creep Model for Salt,’’ The Fifth International Symposium on Salt, Hamburg, Germany, 
6. Nickell, R. E., ‘“‘Rock Salt Constitutive Modelling,” Report UCRL-15207, Lawrence ; 
Livermore National Laboratory, Livermore, Calif., Mar., 1980. 
s 7. Obert, L. E., ‘‘Deformational Behavior of Model Pillars Made from Salt, Trona and 
Potash Ose, Proceedings of on mock Mechanics, Rolla, Mo., 
Aug., 1964. 


i 
q 
= 
— 4 q 
| 
| 


WW 


4 

ink 
nO 


the pr 


rd 
o ON PEAT: OBSERVATIONS AND DEsicn 


Edward T. T. Hanrahan,’ and Martin G. — 
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a peat bog of average depth 25 ‘ft 6 m) ¥ was reconstructed “during 
period between August 1953 and March, 1954. The thickness of gravel base, 
originally 12 in. (305 mm) was to be increased to 2.5 ft (0.8 m). The control 

of the thickness of gravel proved to be difficult, and excessive quantities were 
used in certain locations resulting it in serious consequences in the f form of 

"settlement, upheaval, a and lateral spread. three lengths in particular, 
3,630-3,820 ft (1,107-1,165 m); 3,905-4,090 ft (1,191-1,247 m), and 4,630-4,783 

: ft (1,412-1,459 m), the rate of deformation was of such magnitude that it was 
decided to adopt the expedient of YN iy (5) i.e., to remove part of = 


The 3.28 2. 19 x 1.64 ni m x m x 0.5 m) 
placed i in one, two, or three layers respectively, as shown in Figs. | and 2 
Remedial measures were complete before November 1954. Road levels were 
~ taken during the ensuing period of nearly 26 years, and there is now an eet a 


An investigation by means Oo} of ten trial holes was undertaken in ‘July, 1962 — 
to determine the thickness of gravel base and the condition of the peat bales — 
after eight years of submergence. At the same time, a comprehensive laboratory — x 
investigation was carried out to make a general study of the behavior of = 

“under load. A simple method of predicting settlements was formulated and ‘J 
a comparison was made with the settlements which had occurred baad to that 
date. This study was published (6). 
_ was observed that in the baled sections which, during construction, had . 
; ee to be the weakest sections, the road was stable and in good shape. 
Elsewhere, it was badly distorted, and has since required continuous addition 
- of material to maintain its surface. A further investigation was made in March, 
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1980, mainly for the purpose of determining ‘the thickness of of ec construction at 


each of the ten previous locations. 


settlement to be separately evaluated nema three headings and the i 
combined to give total settlements, 


a Elastic Settlement.— Various formulas and charts based on n Young’ s modulus. 


| 


son’s 


and Poisson are 
~~ settlement was made | during the 1964 investigation using plane ne asin 


deformation in a model tank. If it is assumed that the elastic properties ar 
_ proportional to the shear strength, the elastic settlement (stiff strip) may then 


be from the formula 


n which Ve: = elastic settlement (m); B = width of strip (m); ¢ 


strip (kPa); and Ss = shear s strength of soil (kPa). 
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—Consolidation settlement is measured with the 


- oedometer apparatus and requires that a number of similar specimens be tested - 
under a range of loadings. Each load is applied as a single increment. The — 
"wellknown time-scaling method, based on the Terzaghi equation | of one- 
dimensional consolidation, proposed by Hanrahan (5,6), is used for forecasting _ 

- field settlements. The method is especially useful for soils such as peat, the 
properties of which change significantly under load and time. a> ae 

e Firstly, the laboratory curve relating vertical strain and time is established — 

(Fig. 3). For forecasting field behavior, the time scale is then altered in the : 
ratio (D/d)? in which D = = maximum: length of drainage path in field; and 

. d = maximum length of drainage in laboratory. It is generally understood that 

4 this method is not applicable in the so-called secondary phase of consolidation. — 
This does not correlate with this author’s experience. 
If the loading does not remain constant, a simple method of curve- hopping \ 
is used, as illustrated in the following example. Suppose it is required to foreca: 
the field settlement of a stratum (D) under a unit load x for a time t, subsequent - © 
to which the loading is to be increased to y for a time interval At. Two — 
single-increment oedometer curves (d) ABC and DEF, (Fig. 4), are soaked. 
The field-time scale is established by correcting the laboratory scale by the 4 
ratio (D/d). The strain produced by the load x at time ¢ is then shown at 


_ the point B. The total strain produced at time t+ At by the loading sequence — 
is given by the point F. The method of curve-hopping is valid both for increasing * 


Creep Settlement. —This important type of settlement is taken into account 


— in which y, a creep settlement (m); B = width “ loaded strip | (m); q = unit 
_ load on strip (kPa); and s = undrained shear strength (kPa). ee ees 
_ This equation is the result of a regression analysis of a family of cl curves 
©). Eq. 2 may also be used when either the loading or the time is a variable. 
_For example, suppose that the creep settlement is known at the start of a 
_ time increment Af, during which a load x is in operation. Clearly, this settlement 
: may have developed under a succession of previous loads of different values. a 
Eq. 2 is then used to find an equivalent time, t., . This is the time it would 
_ take for the known creep settlement to develop under the single load x. The © 
required settlement at the end of the time interval, Af, is nd 
t=t,, + At, and again bq . 2 for a new value of y,. 


Design data are used al the application o of the simple method to the - County 


_ In-Situ Characteristics.- _ were previously described (6). The average a) 


~ undrained shear strength was 0.75 psi (5.81 kPa) and the moisture content 1,150%. _ 


A relationship between undrained shear eng aioe — content was also 


| 
| | 
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“ROAD. (ON PEAT PEAT 


the: simple assumption ¥ was ‘made that there was a uniform vertical strain, a 

Fy reduction of moisture content, and an increase of strength throughout the full 

rf -. of the stratum. This is the basis of Fig. 5. 

_ Oedometer Settlement Curves.—Fig. 6 shows three curves. relating vertical 
Strain in the oedometer with time, obtained using | a standard 76 mm diam pe a 


loads may be as shown 6). 


| 


‘STRENGTH 


hand 
Zoe 


6.—Oedometer Curves Effect of Varying Loads 


> 


7, 

Thickness of Peat Layer. —= 25 ft (7.6 m). The ‘stratum was underlain 

_ by impermeable glacial till, and the maximum drainage path was, therefore, 
“taken equal to the thickness of the pent stratum. Thus, | min laboratory time 
= yn 3,500 and 9,500 days. Any change in n loading ‘was assumed to *~ 
be caused by periodic maintenance, and observed increases in thickness lol ae 


shed O24 . Strenotheninge which tabec nlace on 
|: 
IG. 
Tye 


pe 
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i 
he 


18.00 
— fie: 


The changes in loading at 100 Pirelli the aaa: lightweight 
_ bales for the gravel. Reference to baled sections, such as at trial holes 2, 3, _ 
a and 6, confirms the stability achieved by the lightweight fill. Unbaled sections, 
eB, g., 4 and 9, exhibit large and persistent def deformations (Fig. 8). ee 7 
Comparison Between Computeo ano Osserveo Serriements 
Figs. 7 and 8 show, respectively the computed and ¢ observed settlements of 
the subgrade at each of the 10 locations. Satisfactory agreement is seen to 
exist between both sets of curves, _ thus confirming the validity of the simple — 
_ Method of e, and e, Parameters. —The simple method of computation outlined — 
“4 previously, although yielding satisfactory results, may be criticized on a number — 
s counts, chief of which is that it relies for the estimation of consolidation 
settlement on the oedometer apparatus. The oedometer cannot measure lateral is 
strains, and does not give | essential information on the lateral stresses. These — 
were shown (10,11) to vary with time, size of specimen, and location of — 
in a stratum. The importance of varying stresses was — demonstrated by 


Biot (1) in his equation _ = 


in which = excess pore pressure; c coefficient of consolidation; 

La Place operator; and @ = o, + + 
a Unfortunately, Biot’s equation is also unsatisfactory, for the following reasons: 


i 
@ cays 100-3,500 days 5,500-8.000 days 
bunds} = |pounds}| =  |pounds 
= lopascals | foot | kilopascals | foot | kilopascals | foot ri  & 
2069] 10.99 | 229 109 «229 
4 6.91- 144 144 10.56 =| 220 
460 
106 
1. Solutions are available for only a few cases, Gibson and McNamee 1987, 


ROAD PEAT 


days | 6,500-8,000days | 


In 
pounds 


square 


26.06. 
10.08 


7 De Jong 1957, Davis and Poulos 1972, Smith and Hobbs 1976. a 4 
a _ The theory is confined to ideally elastic materials of constant properties. 7 


a ° ’ _ The equations do not permit a ready conversion of pore pressure necessary 
for estimates of settlement to the components of principal strain. = 
_ 4, The theory presupposes that strain can be defined solely in terms of effective 
‘Stress. It was shown (9) that this is apes only with an equal all-round stress 
_ 5. The theory is defective in its ignorance of time-dependent distortional strain. 
Most of these defects are eliminated by the n method of the e,, e, parameters 
i a. 9,11). The following is a very brief outline of this method: Total principal 
: Strains e,, e,, e, are subdivided into two categories. These are termed the 
and e, strains (or parameters) and are measured i in the 
_ laboratory. Thus 


the e, ‘strains are those v which c can develop ‘only under drained conditions. Their 
= is determined by the consolidation (diffusion) process. The e, strains can 
occur in both drained and undrained conditions. Only in saturated and non-dilatant | ’ 4 
soils are the e, strains synonymous with distortional strains. In other soils they © 
include a component of volumetric strain, although in soft compressible soils — 
the latter is usually negligible. The rate at which the e, strains develop i is governed 
the theological properties of the Unlike: ‘the strains, this rate is 


4420 
= 
4 
0,UUU-9,9UU GayS 4 
| pounds | ounds | lla 

| square | | square 
kilopascals | _—foot kilopascals oot 

(10.99 10.99 | 229° 

648 | 135 6.86 7.20 
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independent of the the size or nines path of the specimen or element. 
9 shows ys a diagrammatic -Tepresentation of both types of strain for 
4 system of axisymmetric principal stresses. To illustrate the important continuous 
interaction of both types of strain, the special condition of — 


7. at Ton s Using Simple 


A Se 


TIME (ays), 


= consolidation | may be examined. TI This i is expressed 

o,>0,= 
The e, = Ois described (Fig. %) 


| 
|. 
) 
a 


Q 


onreat 

n which ¢ = time; and z = dep th. ) 

” Ea. 7 implies a continuous equality between two types of strain of opposite a 


a But, under a constant stress system, the two types are known to develop _ 
at different rates. Therefore, in circumstances where the strains cannot develop 
_ in full freedom, the stress System is subject to continuous change. This has Pp | 
been demonstrated experimentally (10). The lateral stress has also been shown > + 
to vary with size of specimen and location of element in a stratum. The special 
-_ savaage of the method is that it permits a continuous evaluation of the lateral si 
stress, o, (11). For one-dimensional consolidation, essentially two laboratory- 
"measured relationships, F, are required. — ra 


7 9. of e.,e, Strains 


_ The first, F,, is based, using an undrained triaxial specimen, on observation | 


iS _ with time of the e,, strain for a range of values of principal stress- difference. 2 

The second relationship, F,, i is based on the | _ observation, with time of the 
volume change of a radially drained triaxial specimen under a range of values 
of effective stress. Equal all-round stresses are employed, thereby avoiding ¥ 
the complication of distortional strain. For applying the measurements, if the 
stress is equal in the three principal directions, the e, Strain may be taken <. 

to be one- -third the volumetric strain. If unequal, correction factors o,/o,., 
o;/0,., are applied for each of the three principal directions, respectively 

-_ Field rates of strain ‘may be computed from the wns rates by applying 


1950 
yo? (AWA! 
LA 


Ina typical example of one- -dimensional consclidation, the following are usually J 
- hewn or specified quantities: (1) Vertical pressure, o,; (2) time, ¢; and (3) 4 
pore pressure, u. The latter i isa specified time. Hence, 


the lateral ral stress, sy I may be solved at all values of time and depth below — 


_ Having thus fully defined the stress, it is possible to solve all the remaining — 


and strain components by using relationships similar to Eqs. 8 and 9. 
Finally, using Eq. 4, the total strains may be determined and the settlement A 


i Application of Method to Plane Strain. —Deformation of soil under a —_— 


s+ ers sum) . (12) 


» is better "Stated | by the equalit 


| + (34) 

As: a consequence of the — Eq. 13a, it follows that the intermediate 

- Principal stress, o,, in the longitudinal direction, is subject to continuous — 

As seen in the Biot Eq. 3, this change in stress has a fundamental influence © 

on the rate of is of pore pressure, and, in turn, on the rate of consolidation — 

— 

Computer .—Full details of a ¢ computer program (PLANE) for 

application of thee, and e, parameters to the problem of forecasting ¢ the ae 
ci a road are available (12). A few comments may be made: | 


‘equation by Scott, (13) (Eqs. 6-62, p. ). 272), | 
& 
_ 2. The time is subdivided into intervals, during each of which the; pore pressure 


is assumed to dissipate while the applied stresses remain constant. The term 


3. At the of interval, the in total stress which | is 
to satisfy Eq. 1 


pore pressure in with the equation of Scott, 13). ved 


‘The following are the data used for the application of the program (PLANE) 


= Data.—A regression analysis of a family of curves, obtained by testing 
ji undrained triaxial specimens under an axisymmetric stress system, yielded the — 


” 
= 14.85 $2.35 
ty) 


tty 


felationship 


It will be noted that the principal (o, is expressed 
in the dimensionless form, t/t, (8), in which 21 = (0, — o;); and 21, = (o, 4 
03) a at failure. The assumption of constant volume during undrained shear 
reasonably valid for peat. Hence bar ore 


oy, (15) 


= 

variation in ins moisture- content is “expressed t 
5 


ave 


in “which shearing strength (kPa); ad in in cc 


~ 
FIG. 10.—Comparison Between Observed and Computed Settlements Based on _ 
Method of e, and e, Parameters a Four Locations 


: 3 triaxial specimen initially 76 mm by 38 mm diameter, at initial moisture-content, — 
+ 150%. For convenience, the e, parameters are measured in the laboratory _ 
using an axisymmetric system of stresses. For use with a nonequal system — 
they must be corrected. The factor to be used is obtained by comparing the 
extended oa of both the axisymmetric and three-dimensional stress tensors. — 


which = required value for nonequal three- stress system; 
and e,,,,, = value of e,, from Eq. 15. 


A 
Data. —Data can be expressed by the 


| 
r 
| 
18) 
i 


in which Av = loss in ~_— (cc) ot triaxial specimen, ‘initially 76 mm by ; 
; 8 um diameter at initial moisture- content, 1,150%. And, c, = 38 in which 


_ Comparison Between Observed and Computed Settlements.—Fig. 10 shows a 
os “comparison between the observed values" of settlement and values obtained | 
g the e, and e, parameters, from four locations. The level of agreement | 
is not noticeably better than that shown in Figs. 7 and 8, which is based on 
the simple method. Nevertheless, the agreement is satisfactory. The availability | 


4 the | ‘computer ‘program PLANE makes the application of this method a 
Two methods of evaluating settlement are described. The first is a simple — 
ia procedure, requiring the addition of settlements under the three headings, elastic, 
consolidation, and creep. The second is based on a general time- eer 7 
constitutive model, relating stress and strain in three dimensions. This model __ 
permits the important change in stress, required by “‘forced strain” = 
Zi such as plane strain, to be continuously evaluated. The complete definition 
of stress, thus afforded, in turn allows the strains to be solved from basic. 
- elationships. Records of settlement of a road, extending over a period of 
_ twenty-six years, indicated good agreement for both methods of computation. 4 
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= width of strip load; = 
= coefficient of consolidation; L 
maximum length of drainage path j in — 
= maximum length of drainage aii in oratory, 
= principal strains; 
principal strain parameters; 
unit load on strip; 
-undrained shear strength of soil 
excess pore pressure; 


= elastic settlement; 
creep settlement; 
change i in volume; 
+0, 
principal stresses; 
= octahedral stress; 
1/2(0, — 93), ‘with reference to triaxial test; 
1/2(o, — ie. princi pal stress difference at at shear failure; 
Lal Place ee. 
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TECHNICAL NOT NOTES 

aided by the feegatact, 41 freqochdy iacressts, TR 

_ INSTRUMENTATION CORRECTIONS TO ' 

By John R. Hall, 'M. ASCE, John T. Christian,’ F. ASCE 
A. M. ASCE, and James L. Wolfgang* j 

The detection of the arrival of seismic stress waves for a field 1 measurement 
of wave velocity requires that the incoming signal be recorded and displayed — 
: by an instrument such as a seismograph. F ilters, introduced into recording circuits — 
a to limit noise or to prevent aliasing in digital recording, can have significant 
effects on the observed arrival times and on the computation of wave velocity. — 

- Such filtration is often not taken into account in processing the data, but the 
_ example that follows indicates how the effects can be evaluated and the data 

corrected by simple techniques. The properties of electronic filters are, of course, 
well known to electrical engineers, and the following comments are an 
with an application in which their importance has not always been apprecintes 
by those trained i in other disciplines. — as ; 7 


Ficter CHARACTERISTICS “AND Errects 


the filter of which i is 


filter i is 


: in = output amplitude; E, = voltage f= 
signal feocaher in Hertz, tz, fi = filter cutoff frequency in Hertz. The phase 
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: One widely used recorder is the Bison Model 1580 
- composed of two first-order, low-pass, Butterworth filters Cascaded in series, 
ies of 160 Hz and 89 Hz, respectively. The ratio of 
he input voltage for a unit gain first-order, low-pass, . 
a 


of the pore Se for the individual filters. For the Bison instrument with | i= 


and a combined phase 6, in 


‘The phase a 8, causes a time lag, At, which is equal to 6/(27/). The 
amplitude ratio from Eq. 3 and the lag from Eq. 4 are plotted + mies —_ 
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FIG. 1 —Response of Filter to Steady State Sinusoidal input 
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TECHNICAL NOTES 


will be seen that with cai frequency | there is a decrease of amplitude — 

and an increase of phase angle. Because the time lag is equal to the phase = 

_ angle divided by the frequency, At decreases as frequency increases. The waves % 
measured in the example given below had a predominant frequency from 100 
Hz to 200 Hz. For f = 100 Hz, the phase angle is 80.3°, and the time lag — 
_. Transient Signals.—The signals tr transmitted in wave velocity measurements © 
a are not simple sinusoids. However, since the system is linear, it is possible 

iA to find the solution to transient signals by means of the Duhamel integral. _ 


Output Dh. 
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(msec) 
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ponse « of Filter to haan T 


for transient signals 
wares from computer 
1000 
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A typical result is s shown in Fig. 2, where the time lag ‘can be seen clearly. 
The shape of the wave form and its dominant frequency were varied to obtain: ‘ 
the values of time lag for a range of transient signals as plotted in Fig. 3, 
which also shows the steady state time lag for comparison. The results between 
- 100 Hz and 200 Hz indicate that time lags from 0.8-1.2 msec can be expected : 
i _ Example.—A recent field determination of shear wave velocity by cross-hole — 
measurements in a competent rock was performed using three borings to " 
by. _ of over 300 ft. The energy was imparted to the walls of the holes by 


identified 
J wave FIG. 
fires motiog 
oe 
| 
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Bison Model 1465 ‘hammer. The Bison 1580 seismograph was to record. 
2 all signals. The surface configuration of the holes is shown i in Fig. 4, but evaluation - 
- of verticality indicated considerable horizontal drift in the | positions of the holes | 2 
hs depth. Corrected horizontal distances were used in all calculations. _ eas 
.: _ The shear wave velocity was originally computed by assuming no time lags — 

in the seismograph and no delay in the trigger timing. The velocity appeared — 

= increase for longer distances (Fig. 4), a result that indicated rene time > 

a a linear regression analysis v was made of the travel ange versus ‘distance from ie 


“all three the The results showed that the mean of 


BORING LAYOUT 


AtoB Regression using data 


A to plus mean time intercept 


B toc raw data 


FIG. 4. and Corrected Shear Wave, Velocities for Rock 
intercepts computed at all elevations was 1.0 msec with a 
deviation of 0. 3 msec, These values are consistent with those predicted from 


g With the mean of 1.0 msec as an additional data point for each elevation, 
ta new linear regression analysis was performed to compute the shear. wave 
_ velocity profile plotted on the right in Fig. 4. Due to the high velocity of nl a 
materials, the difference between the corrected and uncorrected velocities is . 
large. For soil sites, the influence of a 1 msec timing delay would be ewe 
*. The time ines calculated in ‘the regression analysis of the field data also contain me 
other effects such errors: in picking wave arrivals, 


| 


induced by the filters lends confidence to the corrected shear wave velocities | 
= Fig. 4. The computed shear wave velocities are affected by errors in the - 
Ee by deviations of the boreholes from the vertical, especially at depths - 
‘greater than 200 ft (61 m). However, since the values of corrected shear wave _ 
yelocity show the anticipated gradual increase ‘ai | depths below 125 ft (38. 
m), the errors in deviation surveys at large depths would seem to be less a 


a Two Nails conclusions can be drawn from this work: 
1. The characteristics of instrumentation can have important effects on the 
; sestihe of field geophysical tests. In particular, the effects of filtration should 
- Oe Regression analysis of travel time versus distance is effective in determining — 
timing error and correcting for it. A isto: oil vo Of 


with reversible polaitey, shear wave arrivals in cross-hole seismic tests can 
be identified unambiguously because the direction of first motion in the shear — 
wave reverses when the direction of the impact is reversed. The direction of 
first motion in the P- -wave i is assumed not to reverse. - Unfortunately, shear 
: wave sources such as the Bison hammer also produce P-waves of opposite 
polarity when the direction of impact changes. 
on The reason that P-waves reverse polarity is that a downward excitation of 
_ the hammer causes initial compression below the level of the hammer and initial 


tension above it. An excitation reverses this pattern so the P- waves 
must reverse polarity. A theoretical argument is that in a linearly elastic material, 
the superposition of equal and opposite upward and downward excitations gives 


“not cancel, and there will be energy propagated with no input.  —™S 
. _ Non-reversing P-waves may be caused by several factors, including late arrivals 
@ re) 


a case of no net excitation. If the P-waves do not reverse polarity, they will ‘ 
of reflected or refracted waves and secondary dilatational effects near the source. : 


However, this misconception probably arises because in the down-hole test, 

_ shear wave energy is often created by striking a horizontal beam or plate on 
g surface in opposite horizontal directions. There is almost certainly some 
—_ component to the blow, which does not reverse when the beam 
= s struck ir in the opposite h horizontal direction. ° This unreversed vertical load causes 


an ‘wareversing P-wave, which can be observed in field data. 


Mullin, and in of Stone & Webster Engineering Corporation. aur 


4 
. energy to the soil, and errors in the trigger; but in this case, filtration is the J 
a 
| | 
a 
1 
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IntropucTion 


ro to a point load, line load, and strip load is being performed using the modified : 
4 psy of the Boussinesq’ s equations. For the strip load type of suschenpe, only 


: Se user to evaluate the total lateral surcharge pressure by: mechanical integration 
procedures which obviously are time-consuming. 
: _ This paper is presented to show that a direct solution for the ‘total — _ 
surcharge pressure can be derived by mathematical manipulation of the variables — af 
in the given equation for the unit lateral surcharge pressure. Simplified expressions a 
for the location of the centroid of the total lateral surcharge pressure and the os 
‘point of maximum unit lateral pressure are also derived to complete the analysis. 
si Retaining wall structures supporting continuous wall footing, highway, and — 
railroad loadings are practical examples in which strip load type of surcharge — 
— is applicable. The calculated total lateral surcharge pressure is then added to 
‘other lateral pressures such as earth pressure, water pressure, etc., , for for stability se 4 
and structural analysis of the retaining wall structure. 
Elastic analysis 0 of pressure a against unyielding anfin from surcharge los load (3) 
_ This value is doubled to fit ex imental data (4). Eq a 
experimen al data (4) Eq. 1 is written as 


o, = — (B — sinB cos 2a) 
‘Civ. Engr., Planning Div., New York City Transit Jay St., 

7 Note. —Discussion open until March 1, 1982. To extend the closing date one month, 
a written request must be filed with the Manager of Technical and Professional Publications, + 
ASCE. “Manuscript was submitted for review for possible publication on January 21, 
- 1981. This paper is part of the Journal of the Geotechnical Engineering Division, Proceedings 
of the American Society of Civil Engineers, ©ASCE, Vol. 107, No. GT10, — 


rons 
+ 
= 


Cartesian ‘coordinates. Fig. 2 shows ‘the notations introduced for t this Purpose, ~@ 
e., make the substitutions thus: 


when 


- Eqs. 3, 4, and 5 in Eq. 2 to yield oe 

x=— -[@, —9,) (sin 0. cos 0, — sin 6, ‘cos 


the notations and “ “g are hencef orth discarded. Introduce 


= 


trigonometric fancti ons in terms of dimensions in Eq. 6 to obtain 


a + 
7 can bei integrated using standard formulas Therefore 


_ Fig. 1 shows the geometric representation of the variables involved in Eq. 2. _ 
g Transform Eq. 2 to a form suitable for integration. To do this, it is necessary _ 
il — —— : 
SAR 
4 


“od 


Failure Plone (0 


2. —Moditied Notations 


Titten 


in 6, = 0, at subgrade level = arctan - 


= 0, at subgrade level = arctan — (13) 

and dimensions bs = uld be included within the a wedge defined by 

For active pressure condition 
6 =(« 


When 6, and 6, are expressed 10 j 


TECHNICAL NOTES 


err, a 

Foenpaniony- 
pei: Problem No. 1 —In Fig. 2 calculate the total lateral alia pore 
against the vertical wall when the soil behind it is in the active pepeeane mode 
& failure given the data: a = 6 ft (1.83 m); b = 4 ft (1.22 m); h = 20 ft 

— (6.10 m); g = 1,000 psf (4,882 kg/sq m); and > = 30°. Solution: foonite 
+ b) < 20 tan (45° - — 15°); tan 6, = - 10/20 = 0.50; 6, = 36.57°; tan 6, 
= 0. 20; 6, = 11.31°. Then by Eq. 11 P = 1,000/90 ‘[20(26. 57 131) 
and P = 3,390 lb/ft (5,044 kg/m). > 


_ The equivalent uniformly distributed surcharge pressure against the wall = 


3,280/20 = 170 pat (630 (he fretion and at (he 
_ The location of the the centroid of the total lateral surcharge pressure is determined 5 


(0, — 6) + (R - 
2h(8, -8 
Problem No. 2.—Determine the centroid of the total lateral 


y wal 
0,) = 26.57° — 11.31° = vera 
‘Substituting in in Eq, 15 


26) + 5,084 — 57.30 


was 
‘The of maximum unit t lateral sur charge is evaluated dif: feren- 


in which = = a’ + 3b(a 


it lateral surcharge pressure and calculate its =r: Solution: ean 
= = (6° = (3)(4)(6 + 4) = ga 2) 
0 = (4) 6 +4) = 6,4 6,400 


.. 
= 3.53 ft (1. 
ened this value in Eq. 7 to get the maximum unit lateral surcharge pressure, 
21,00) [ ( 353 3.53 
— arccot —— — arccot —— 


= 


1 psf 1,762 k 


The use of the derived formulas will save the user time in computation of 
surcharge pressure involving the strip load type of surface load. oh an oe ee 


idere 
applicable to yielding o1 or retaining wall structures with the soil 
_ behind it either in the active or passive mode of failure depending on the given 
oe -- of the problem. (The lacie derivation may be furnished by the 
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tiating Eq. 7 and equating the ro and solving for Y ihe, 
Y,, = | 

] ; 

q 
i 
: 
| | 
— 
= 
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FOUNDATIONS UNDER VERTICAL wine 


Foundations, especially those for machines, are usually embedded. For 


embedded footings, damping exists from the friction and adhesion at the side - 
interface of the footing and the soil, in addition t to that due to radiation. During © 


4 


7 ee damping and hence it can be neglected. Exact solution for steady-state 
vibration of a system damped by a combination of dry and viscous friction 
Fig. 1) has been given by Dem Hartog (2) for frequency independent excitation. 
_Anandakrishnan and Krishnaswamy (1), following the same procedure, gave 
: a solution for the case of frequency dependent excitation. The exact solution 
_ given in the previously cited sources is very complicated and lengthy. aa 
In this note, a simplified procedure based on energy consideration, has been 
developed, for the solution of steady-state vibration of a system with combined 
viscous and Coulomb friction damping, subjected to frequency in dependent — 
- frequency dependent excitation, which yields results essentially same as 
- the exact solution . The Proposed method uses equivalent damping (3) which 
. assumes that if the damping in a system is small, the total damping effect 
can be represented by that of an equivalentdamper, 
i, Mass-Spring-Dashpot System with Coulomb Friction—Simplified Procedure.— 
_ The equation for steady-state vibration of a mass- eS ~dashpot : ees with 
Coulomb friction n (Fig. I) is given by isi 


MX +CxX+ F= 


in which M - = = mass; = viscous | K = 
F = Coulomb friction, and Q(t) = Q,. Sin (wt + ) = applied dynamic force, 

Q. = force amplitude, w = circular frequency, @ = phase difference of the 

dynamic force, X = X, sin wt = at any time a = Gieplacement 
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- translational movements, radiation damping is very large when “compared with 4 


4 
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Sisipate AE, due to viscous ous damping per cycle is 


by viscous damping per cycle of vibration is given by Baas, ve ihe 


is given by 


4s; =4FX,. 


= onsite 


1 Spring | Model with Combined Dry and Viscous Damping 

mass- motel can be written as ye 


. (8) 


— 
| 
&g cycle due to equivalent viscous damping C,, (3) is given [im 
| 
q 
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TECHNICAL NOTES 


te factors M,, = X,M/me, in m = the eccentric mass; 


and e = the eccentricity, for nc dependent excitation and M. =X, K/Q. a 
frequency independent excitation are given by 


Ty 


Mew 


ne qe one 


+ 


" 
a5 | 
J 


‘in which damping factor D = steady. -state vibration, the 
~ 


governing conditions are 


for frequency d dependent and and frequency independent exci excitations respectively. 7 


4 Elabor co 
ait 
- 
4 
|_| 
q 
independent and dependent excitation with frequency ratio, w/w,, for D= 


OCTOBER 


0.2 and for different frictional force ratios—(F/Q,) or F 

the simplified method presented in this note. They are com the results 
obtained from Den Hartog (2) and Anandakrishnan and Krishnaswamy q af 
respectively. It can be seen that the results obtained by the proposed method _ 

compare very well with the exact solution, 

Figs. 4 and 5 present the relation between the magnification factor at resonance, 


and the frictional force ratio ‘obtained from the proposed method, and from rom 


— 


FIG. 2.—Variation 
Factor M_ with w/w, for Various Values “7 


of F/Q _—Frequency Ex- 
Citation 


a 


Ss 4. —Variation of Magnification FIG. 5. —Variation 
Factor at Resonance M,, with F/Q, for _—* Factor at Resonance M,,, with F/mew? 
_ Various Values of Damping Factor D— isk for Various Values of Damping Factor 


agnification 
rious Values 
pendent Ex- 


woda 


FIG. 6. = n of 0 /e at Resonance » with F/Q, for Various Values of Damping 


Factor Independent Excitation | 


FIG. 2 —Variation of w/w, at Py with F/ mew? fer’ Various | Values of of Damping 


| 
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those obtained from exact solutions for different The proposed — 
: method gives essentially the same results as obtained from the exact solution. — 
The strong dependency of the magnification factor with the frictional force — 
ratio is evident and the relationship with the frictional sh pape ratio up toa il 


of about 0. is almost linear. Peres 


excitations for different damping factors. For lower values of frictional oe 
ratio, the difference between the pooposed method and the exact solution is 


_ The new method could be used for solving problems involving vibrations | . 


of embedded foundations subjected to other modes of vibration. * 
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Discussions 
4 Discussions may be submitted on any Proceedings paper or technical note published i in any : 
Journal or on any paper presented at any Specialty Conference or other meeting, the ecmangy 
‘of which have been published by ASCE. Discussion of a paper/technical note is open to 
anyone who has significant comments or questions regarding the content of the paper / technical 
- note. Discussions are accepted for a period of 4 months following the date of publication < 
of a paper/technical note and they should be sent to the Manager of Technical and Professional q 
Publications, ASCE, 345 East 47th Street, New York, N.Y. 10017. The discussion period may __ 
extended by a written request fromadiscusser, 
The original and three copies of the Discussion should be submitted on 8-1/ 2-in. (220-mm) 
by I1-in. (280-mm) white bond paper, typed double-spaced with wide margins. The length of ‘ 


a Discussion is restricted to two Journal pages (about four typewritten double-spaced pages 
of manuscript including figures and tables); the editors will delete matter extraneous to the 
Sam under discussion. If a Discussion is over two pages long it will be returned for shortening. 
‘All Discussions will be reviewed by the editors and the Division’s or Council’s Publications 7 


- Committees. In some cases, Discussions will be returned to discussers for rewriting, or they | 
may be encouraged to submit a paper or technical note rather than a Discussion. = = ~—- a 


- Standards for Discussions are the same as those for Proceedings Papers. A Discussion is 
_ subject to rejection if it contains matter readily found elsewhere, advocates special interests, 
is carelessly prepared, controverts established fact, is purely speculative, introduces personalities, a 
or is foreign to the purposes of the Society. All Discussions should be written in the third — 

; Ss and the discusser should use the term ‘‘the writer’’ when referring to himself. The a 
author of the original paper/technical note is referred to as “‘the author.” 
- Discussions have a specific format. The title of the original paper /technical note appears 
at the top of the first page with a superscript that corresponds to a footnote indicating the 
"month, year, author(s), and number of the original paper/technical note. The discusser’s full - 
should be indicated below the title herein an with 


‘The discusser’s title, company affil and should appear on first 
4 page of the manuscript, along with the Proceedings paper number of the original paper / technical 
; note, the date and name of the Journal in which it appeared, and the original author’s name. i 
_ Note that the discusser’s identification footnote should follow consecutively from the original 
_ paper/technical note. If the paper/technical note under discussion contained footnote numbers — 


| and 2, the first Discussion would aa with footnote 3, and ‘Subsequent Discussions would — 
a Figures supplied by the discusser should be designated by letters, starting with A. This also 
applies separately to tables and references. In referring to a figure, table, or reference that 
_ appeared in the original paper /technical note use the same number used in the original. = 
4 It is suggested that potential discussers request a copy of the ASCE Authors’ Guide to 
the Publications of ASCE for more detailed information on poopenetion =e submission of 
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= experimental technique developed by the authors s seems to be extremely _ 
interesting. However it does not appear that the results obtained can only be | 

_ analysed qualitatively as indicated by the authors. In fact a granular material 
can be as isotropic i in so f far as its uniformity coefficient is superior 


4 joint of the | particles i is , stochastic, which 2 means that the mean number of contacts = 
_~per particle, N., is equal to six. This is due to the fact that the structure 
_ made up by the particles is isostatic if the presence of inactive particles is 

ignored (these being the particles which can be extracted from the mixture 
The writer thinks that there is no definite relation between N. and n. The 
jaa Ni, depends on the arrangement of the particles (stochastic, cubic, 
the pa . . .); m depends on the granulometry of the material, the shape 
_ of the particles, their mineralogic nature, the size of the mould containing the 
material, and essentially the energy exerted, including the specific weight. __ 

_A number of contacts : superior to six probably results from the fact that — 

in some places i in the material, geometrical arrangement have occured. Besides, 

a number of particles do exist in stochastic structures which can be very close 

to one another without being in contact. The experimental technique and the _ 

viscosity of paint might explain the fact that such cases are considered as real _ 

G Finally the writer wishes to congratulate the authors for their fine study 0 on | 
the relative » displacements between particles. The | experimental technique s seems 
co promising and should lead to better approach in the knowledge of the 


| itt problem of the relative movement of particles in granular material. 


by Ro ert D. Holtz,“ M. ASCE 


ORE. 


The author’ s observation that the product of the number of contacts per 


particle, N., and the ‘porosity, is approximately three is mos ‘most intriguing, 
and it prompted the writer to look at some other data. These data suggest © 
5 that the relationship may not be so simple. For example, Oda (4) in his studies 
of coordination number (average number of contacts per particle) and shear 
strength reported on some experiments with glass spheres. These data are 
- summarized in Table 1. Only in a few cases does the product, N. x n col 
May, 1980, by Dimitri Athanasiou-Grivas and Milton E. Harr. 
Asst. Dir., Dept. of Civ. Engrg., Univ. of Clermont II, Clermont-Fd, 
“Assoc. Prof., School of Civ. Engrg. Purdue West 
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er TABLE 1 —Date on Glass Spheres wom Oda (1977) 


rt 


 ters(5.2mmand | 


od 
Four different diame- 


- (5.2 mm, 12.3 | 
| 


‘mm, 16.4 mm, and 
247mm) | 


Date on Granules Soils fre from M: Marsal (1973) 


‘TABLE 2. 


ficient 


7 
Material 
Pinzardaran sand and gravel 29 0.22 36 rt 
Malpaso conglomerate 8 | 042 | 0.30 59 
Mica granitic gneiss gradation) 
Francisco basalt grade. 


Francisco basalt grada- | 
Granero slate gradation A 
El Infiernillo diorite | 
MM Granero slate gradation B 
‘Mice grantic gneiss gradation 


— 
i 
Two different diame- | 0.612 Pap 
| 
| us} o4 | 031 | 564 | 175 
0.96 | ILS | 0.34 0.25 | 671 | 1.70 
10 | 98] 045 | 031 | 5.94 | 1.84 
25 | 7.7] 048 | 0.32 
60 | 4. (0.64 |) 
00 | 24] O58 | 037 | 996 | 366 ff 


a 5), approach three most of the time, especially for iets packings, the product 
_ is significantly less than three. From a comprehensive study by Marsal (3) of 
ten granular materials primarily from rock-fill dams (Table 2), it appears that “a 
- the average value of the product, N. x n (Col. 7), may be two or less, i 
‘These data are plotted on the authors’ Fig. 9. Neither the linear nor the 
s proposed hyperbolic relationship fits these data very well. Can the authors — Pa 


materials that different from ‘the ‘materials previously e examined by the 
_ and others? Since the deviation seems to be greatest for very loose and a - 
igs arcing Yo. dire sirineg tsa 10) 


iol 


Linear relationship 


 Porosity,n % 


nN 


et contacts per partici, Me gan agate does 
FIG. 9.—Relation between Porosity n and Number of ( Contacts per Particle N.. , with 
additional data by Mi Marsal (1973) and Oda (1977) 
dense packings, should relative density or density index somehow be c considered 
3. Marsal, R. J., ‘‘Mechanical Properties of Rockfill,’’ Embankment-Dam Engineering, 
2igal Casagrande Volume, R. C. Hirschfeld and S. "J. Poulos, Eds., John Wiley & Sons, : 
New York, N.Y., 1973, pp. 109-200. 2990 
> 4 Oda, M., 7 “Co- ordination Number and Its Relation to Shear Strength of Granular 
17, ‘No. 2, June, 
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Closure by Dimitri A- Grivas,? A. M. ASCE and Milton E.. Harr,* M. ASCE | 


The writers thank both Holtz and Gourves for their interesting comments 
on the experimental technique employed and the results obtained i in this study. = 4 
Indeed, some of Oda’s data mentioned by Holtz are ‘different from those § 
examined as they correspond to packings of glass spheres made under special * 
| conditions. The writers were aware of these results at the time of their study, 
_ they chose however not to include them in their analysis. How could there 
x exist a real soil with ¢ an n average value ow the number of c contacts per times 
be Ba 
that for each particle with of contacts, N. =2, there exists 
particle, N. = 1, and that for many particles with N. = 2 there also exist 
l many particles with N. < 1. Such particles certainly do not abound in nature. — 
_ Even in the case of ideal packings of identical spheres, the oie el 
a Marsal’s data on the other hand represent only rough estimates on N, (Marsal, 
1973), obtained with rock boulders. The examined packings covered a wide 
range of values for the uniformity coefficient and the results reflect mainly - 
the effect of the shape of the boulders on the N.-n relationship (n = porosity). 
_ The comment made by Gourves that the viscosity of the paint may introduce | - 
“les enaiins is valid and extreme care must be taken when the 
_ experimental technique is applied. The pattern depicted on the surface of the 
_ particles after shearing may provide a verification for the true existence of — 
_ The particle arrangement (packing) and the associated average number of 
are important factors in the deformation and strength behavior 


in ‘such media i is a a process involving the destruction and generation of contacts 
- (with or without particle breakage). This process takes place in discrete stages, 
.- Stage occurring when the shear stress at one or more contacts exceeds 
a thrusthold value (i.e., the available frictional resistance). Although there exists 
no physical relationship between the number of contacts per particle and porosity, 
oy relationship between their statistical representatives, N. and n (n is also an 
average quantity) may be sought through a statistical analysis of measurements 
_ There are many factors on which N. depends, the uni uniformity coefficient 
mentioned by Gourves being one of them. To uncover the effect of each of — 
these factors on the N. -n relationship and, subsequently, on the deformation — 
and strength behavior of particulate media requires a considerable « amount of 
=a as well as theoretical study on the subject. 


{ | 
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“Prof., School of Civ. Engrg., Purdue Univ., W. Lafayette, Ind. 47907. 


TO TRAVELING 


= _ The following corrections should ginal paper: 
_ Page 185, Line 4th from bottom: Should read ‘‘no longer” instead of “not 
wa 


* 


188, Line 1: Should read ‘‘Eqs. 19” instead of “Eqs. 
189, Line above Eq. 29: Should read ‘‘Eq. instead of “‘Eq. 
190, Line 2: Should read ‘‘approaches zero” instead of ‘‘approaches’’ 
February, 1981, by 


| 
187, Eq. 15: Should read instead of 


of madly 
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